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PEEPACE 


The current acoustical literature provides us \vith numerous data 
and long tables concerning the sound absorption coefficient of the 
various materials that are on the market. These data together with 
the physical features of the material (the thickness of the layer, 
the porosity, the diameter of the pores, its stiffness) pro^dde the 
.skilled designer of absorbing materials \vith some idea as to how 
new materials can adequately be designed. In this book we have 
tried to give the design a more scientific basis. The principles 
underlying the wave propagation through media — porous or not — 
are described. The wave propagation through porous media has 
been treated at considerable length, being of great impoidanee for 
the majority of absorbing materials. 

Except for the last chapter (VIII) all theoretical considerations 
and measurements are confined to normal incidence of sound. 

With a view to the increasing application of absorbing materials 
behind perforated panels, due attention is paid to this .subject 
(chapter VII). 

The origin of the book has been a reque-st of the editor to the 
first author for a book on the .subject of acoastical materials. The 
latter wrote the te.xt in the grim war winter of 1944 — 1945 when 
all laboratorj’ work was utterly impossible. Due to a serious shortage 
of printing facilities the actual printing was delayed so much that 
it was eon.sidered necessary' to review the text on account of the 
publications that had become available .since the war. At the first 
author’s request the .second author accepted the ta^k to do this. 
It turned out that this meant hardly le.s.s than revmiting the whole 
text. IMoreover the chapter on resonators was not included in the 
original scheme and it due to the initiative of the second author. 

As far as we know a book with a similar aim docs not exi.st. 
which may seiwe as a jastification of adding a new book to the 
exten-sive acoustical literature. Furthermore, the existence of the 
latter may help to meet certain shortcomings in this book, of which 
the authors are convinced there are. Way it, nevertheless, turn 
out to be valuable for the scientific manufacturer. 


August 1949. 


C. ZWIKKEE 
C. W. KOSTEN 
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CHAPTER I 


SnffLE THEORY OF SOUND ABSORPTION BY 
HOMOGENHJOUS LAYERS 

5 1 PROPAGATIOJC GONSTAN-T AAT) -tt'AVE niPEDAN-CE 

If we consider a plane sound wa%’-e travelling in the direction 
of the positive x-axis in a homogeneous isotropic medium extending 
to infinity, the sound pressure depends upon the time t and on 
the distance a: as a damped sine; 

p(x) —Aexp^ ja>(t — ax\^ (1-01) 

in which j- — — 1, <u = 2 • frequency, c = velocity of propagation 
of sound, and exp (...) is the well-known sjunbolical notation for 
At the site a: = 0 we have p(0)=A exp jat. Putting 
a,/c=y3 and a + j/i = y, we obtain for the damped sine the 
shorter analytical form 

p(x) =p(0) exp i—yx). 

The constant y, which, apart from its dependence on u>, is 
determined by the nature of the medium is called the propagation 
constant of this medium; its real part a is called the attenuation 
constant and its imaginarj^ part p the phase constant. 

In the same way the velocity of specific volume displacement, 
V, can be expressed by a similar fonnula 

v(x) =n(0) exp ( — yx), 

in wliieh the same constant y enters as in the formula for p, 
because in a travelling M’ave the ratio p/v must be independent 
of X. "We shall use the symbol v for the amount of material volume 
passing through a unit surface in unit time, or the volume current 
density. Only in a homogeneous medium (free air, compact solid 
medium) is v identical with the material velocity. For example, 
in a porous medium with rigid solid skeleton, the velocity of 

Sound abgorbin" materials t 
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\olumc disphcement t is smaller thin the %eloeit% of the \ibratni" 
air, the ratio being equal to that of the soluiue of the accessible 
holes to the total volume of the medium, 1i The constant h is 
called the porosity or roiity factor and is one of the cle^lentar^ 
properties of the material whieh will ph^ an important role in 
the following paces 

In analogj to electrical practice the quotient 

z{i)=p(x)/i(i) (102) 

IS called the specific acoustic impedance at tlic site x For an 
unlimited medium, 2 must be independent of x This impedance 
IS a material constant and is called its iinic impcdoiicf, repre 
sented b\ ir As p and t are gtncrallj not in phrsc with each 
other, ir IS a complex qnantitj 

If a periodical pressure p(0) is applied at the site j-wO of an 
unlimited mc<lium, the dependence of p and i upon x and i is 
full} determined b} the two quantities y and IV Thus y and 
full} determine the acoustical bebaMour of the medium 

i S ^iI^ED\^CE OF A LA^EB OF FINITE THICKh'F^S 

Wc shall next consider a la}cr of uniform constitution, deter- 
mined b} the properties y and IV, haimg a tliiokne<is ^( 2=0 to 
js=I) Supposing it to be loaded at x = l b} the arbitral 
complex impedance r, , what is, then, the impedance z^ at ® = 0 
(Fig 1)! Part of the sound wa\c will be reflected at the end 
r = so that p is the superposition of an incoming wa\e and a 
reflected wa\e 

p(x) =p,c\p {y (I — x) } + Prcxpl — y (I— x) ) 
t (i) = (Pi/lV) O'^P { Y (f — Jr) ) — (Pr/ll ) e'^P { — y (^ — 1) } , 
in which Pi and p, are ob\ioiisl} the possibl\ complex; pressure of 
the incoming and the reflected wave nt x = /, just inside the haver 

The pressure adds as a scalar, the velocity as a \cctor, hence 
the minus sign in the formula for t ^\s a lioundarv condition 
we have p(I)/i(l)=r, Witli the aid of Has it mn\ easih be 
verified that 

Pr/p.= Ua— 11), (103) 

and introducing this into the preceding equatioas one fimls for 
the impedance at the site x=0 
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cosh yl + W sinh yl 
2 , sinh yl + ir cosh yl ’ 


(1.04) 


a formnla wellknown from the theory of electric cables and filters. 

This formula contains, as a special case, that of the medium extend- 
ing to infinity. Then z, = IT and the formula gives z^ = z, = IT'. 
Another special case is that in tThieh z, = co , in which ease 


2i = TT coth yl. (1.05) 

This latter ease is of great importance for the investigations, 
because this is the one which can easily be realized by loading 
the sheet of absorbent material with a completely hard back wall. 
Bacldng absorbent layers with a rigid wall is a normal way of 



Fig. 1 

Laver backed bv an impedance c. 


applying these materials in praetiee, so that this case is at the 
.same time of great practical importance. 

The other limiting case is obtained by putting z, — 0. Then 
fl.04) yields 

Zi = TrtanhyL (1.06) 

Because the impedance of a layer of air with a thickness of one 
quarter of a wave length, backed by a rigid wall, is zero (see 
(1.05) remembering that y = Jo/c, since there is no damping), 
this case can be realized by placing the absorbent layer at a 
di.stance of -V from the rigid wall. 

It is now possible to calculate in principle the acoustic im- 
pedance of any combination of sheets of different characters and 
lengths by the .sueee,ssive application of (1.04) (Fig. 2). One .starts 
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bj computing z^-i from then 2 ,_« from , and so on 
until fmaU^ 2 , is found The calcuhtion is complicated and tin. 
some and might better be replaced b\ a geometrical method * 

} 3 iruE niPEDWCF n*. of rmE ur 

The theoretical deduction of the quantities 11 and y for am 
medium is aluajs accomplished in cssentialh tiic same uaA In 
starting from the equation of motion for the aibrating medium 
and from the equation of conlinnitj (eonscnaition of mass) 'We 
shall now consider the simplest case, \iz, tliat of free air without 
taking into account the effects of damping 

The equation of motion is found bj appU mg X e w 1 0 n s 


2. 2, 


Ji. 


Fig 2 

SufCcwiTO computation of lie impedance of d multilayer mstem 



equation (force = mass X acceleration) to a thm lajcr of mr 
of thickness dx 


Sp 

3 ( ’ 

in which p^ = densitj The equation of continuita is 


(107) 


in which 


Eliminating t 


^ I «?/y 

aTTF’ 


(lOS) 


A'o = - 


dp 
dp/p« 

differentiating (1 07) 


with respect to x 


Poo ep Feldtkellcr a lerpoltheonc and Chajtor M }" Fip C“ 
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and (1.08) with i-espect to t and equating the expressions for 
S -v/9 X dt thus obtained, gives a differential equation for p 


== P» 

dx^ 7t„ dt^ ' 


(1.09) 


Since we expect a solution of the form p — Aexp {jot) exp ( — y^x) 
we may put: — 3/02: — yo and d/dt — ja>, so that (1.09) reduces to 


from which we get 

y,=dzjo,v^;jK;. ( 1 . 10 ) 

By comparing (1.10) with (1.01) we see that the physical inter- 
pretation of the constant V' K^/Po velocity of propagation 

Co of sound waves in free air, the plus (minus) sign for yp having 
to be taken for waves travelling in (opposite to) the direction of 
the positive 2 :-axis. Substituting tliis value of — d/dx together with 
jo) for d/dt in either equation (1.07) or (1.08) yields for the wave 
impedance 

p/v = IFo = = PoCo • (1-11) 


Prom (1.11) we learn that the wave impedance is real. By in- 
serting the known values of p„ and Cp at room temperature it 
appears that IFp has the approximate value of 420 kgm~2 sec”^ 
(42 cgs). 

If plane waves from the air impinge upon a wall of specific 
impedance z, reflection udll take place; the ratio of the pressure 
of the reflected wave to that of the incoming wave just before 
the reflecting boundai’y is given by (1.03), which runs in the 
present notation 


P’- ■ ■ 

Pi z + IVp • 


( 1 . 12 ) 


We shall call this ratio the complex I'eflection coefficient. By 
.squaring the absolute value one obtains the reflection coeffi- 
cient for the energy, which is eomplementai’j' to the absorption 
coefficient : 


! 


2-irp 

I p^ 

1 ^ i 

s + 11^0 


ap = l 


2 


9 


(1.13) 
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the index zero being added because this coefficient applies onl\ 
to normal incidence 

"Wliat IS said of free air m this section will be sliown to be of 
quite general aaliditj The equations of contmiiita and motion 
maj in manj eases be wTitten down in evneth the same form 
(107 and 108) although the constants p and K avill lia\e a dtf 
feront meaning thea will be dependent on seaeral more or less 
clementan properties of the medium and fmalh will lu\c to It 
considered as complex quantities 

In order to a\oid confusion all sjmbols relating to free air 
avill when ncccssarj be proaided with an index zero as was done 
in this section po Vo o <^o 

A verj important a\aj of application of absorbing materials is 
the direct fixing of a laaer of such material to a rigid bael wall 
In this cose the impedance at the front side was shown to be 
(cf § 2) 

2 = 1\ eothy? (1 Od) 

Substituting the general expressions for and y of (1 10) and 
(1 11) jnolds 

2 =l'^Ap coth Jail ]■' p/A (1 !■!) 

as a general expression for the impedance of an> laacr backed ba 
a ngid wall lie haae to make an exception for some poroiw 
flexible laaors The question of the computation ot A and p 
aaill be considered in § § 7 to II § 6 

j 4 GEOilETRX IX THE COilPLFA PL\NF 

The understanding of such formulae as (1 Oo) and (111) maa 
be facilitated b> plotting the function coth yl and other complex 
functions e of u in the complex plane If we allow « to assume 
all aalues from 0 to the function under consideration fofibws 
a certain contour in the complex plane This eurac as a rule 
appeals more directlj to our imagination tlian does an nnaljtieal 
formula The point is to separate the ical from thi imaginan part 
so that 2 takes the form r(M)=x(*i) + jjM and to plot 2 in 
the plane avith x and y (Iwlongiiig to the same aaliic of the pani 
meter «) as coordinates 

For example the function 2 = 1 — }>• is a straight lino stnrtin. 
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at the point a: = l of the real axis and going downwards (Fig. 3) ; 
z — exp (jia) =eos<a + j sin (a is the circle -K^th unit radius. 
Multiplying a complex number by exp (ja>) means turning it 



Comple.x representation of the circle involute 

round the origin through an angle <o. So by multipljung the point 
A (Fig. 3) by exp (ja) we come to the point B and it is obvious, 



Kg. 4 

Complex representation of a hj'perbola and a lemniseate 

that B lies on the involute of the circle; hence z = (1 — ji») exp (ja) 
is the circle involute. 

The function z = ]/ 1 + > is represented by a rectangular 
hyperbola (Pig. 4), which may be proved as follows. Put z — x + 
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jy = Kl + jo and squnrc — y* + 2 j x!f = l + ju> Separatin': 

the real from the imaginarj parts shows us tint x® — y* = l and 
this IS the common formula foi the rectangular hjperbola Further 
2xy=b>, now since u is alwajs positive x and y have the same 
sign The function 1 + is represented b> the upper branch 
of the hjperbola, the lower branch being represented bj Kl — jio 
The generalization z = vrherc /(«) is an arbitrary 

real function of w, leads to the same hjperbola but the o scale is 
changed This remark is one of a general ^ahdltJ, we niaj al 
wajs pass on to a noir «>-scale on the same curve bj substituting 
fM for 

Inverting a compIc\ numlior z=M exp (j^) gives 
( — if), so the modulus is inverted, and the sign of the argument 
IS the opposite of the original one As the Icmniscatc is the inverse 
curve of the hjperbola the former is represented bj the functions 
(1±;„)-1 

Bv inverting a straight line, one obtains, ns is well known from 
olementarj geometi^, a circle In Fig 3 the inversion of the half 
line 1 + jw IS drawn as a half circle running from the point x = l 
on the real axis towards the origin 

■More examples illusircatmg (he simple geometrical representation 
of functions In contours and vice versa maj casilj he obtaineil 


f 5 atjyMFTRIOXL R^I»R^S^-^T\TIO^ OF COTtl 


In Fig 5 the function coUi yi is plolte<l with the aid of* the 
formula 


coth y? = coth (a + jP)1 — 


sinh2al — jsm2^t 
cosh 2 «Z — cos 2 pi 


(115) 


Sex csxiftsni ftf ^ j>nd Jt JAwkv a sujjtjJ 

And indeed for large values of the argument, coth yl tv opproxt 
mated bj a logarithmic spiral We arrive at this result, becaiLse 
for large values of yf, coth yf ~ 1 + 2 exp ( — 2y0 

Leaving aside for the moment the shift over the dustance 1 
111 the direction of the real axis, the function 

> J Kvbncr Voniogram« of Compkx Iljpcrbolie Function*, p. 2o 
Copenhagen 1017 



§ 5 GEOiTETEIOAI. EEPBESENTATIO^’ OF COTH I 9 

2 esp (— 2 y?) =2 exp {—2al) - exp (—2 j^l) 

represents a spiral witli modalxis M — 2 exp ( — 2 al) and argument 
c = — 2^L Hence tiie relation bet-u-een II and c is 

and this is. for a constant value of a/p, tbe usual formula for 



Fig. o 

Compl-es representatioa of coth -,l 


the logarithmic spiral in polar coordinates. It is indicated bv a 
dashed curve in Fig. 5. In the same wav it can be shown that 
tanhyl converges towards the logarithmic spiral 1 — 2 exp ( — 2 yl). 
The geometrical interpretation of the constant a/p is this. 



Fir. S 

The slope of coth yl 


that it is the tangent of the “slope" of the spiral, which is 
constant along the logarithmic spiral (Fig. 6). The greater a. the 
greater the slope, and the more rapidlv the fnnction converges 
towards its apex. Figs. 7 and S represent two cases with slopes 1 
and OJ. respectivelv. 
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For an arbitrary mednim a and are both functions of o and 
so IS the slope a//3 Taking I as a constant, yl \anes because of 
the variations in a and /? with » In this general case we get a 



Fig 7 Fig S 


A eotli with slope 1 A coth with slope 0 1, 

stroBg damping- TPsak damping 

modification of the logarithmic spiral with a slope \ar>ing along 
the curve In Fig 9 the case a//J *= const /u is illustrated, the 



Fig 9 

A coth witJi a slope mrersely proportional with frequency 

slope approaches the talue of zero for large values of tu the con 
tour approaches to an asymptotic circle 



§ 6 G'SCMETHICAL EEPRESENTATIOX OF FORMXTLA (1.12) 11 


Multiplication of coth yl b 3 ' a complex number IF = M exp j<p 
(cf. 1.05) means enlarging the contour bj" a factor M and rotating 
it through an angle <p around the origin, as in the example of the 
circle involute. Multiplication by a complex number Tl^, which in 
itself is a function of is accomplished b}' appljdng this process 
for each point of the spiral, taking together corresponding points 
(for the same value of <d) of the multiplier IF and cothyZ. If IF 
does not run fast through the diagram, something resembling a spiral 
may stiU be seen; if, on the other .hand, IF runs fast, the spiral 
appearance may be lost, and a kind of damped sine cun'e with 
cunmd axis maj’’ result, this curved axis being the contour of IF. 


5 6 GEOilETRICAL EEPRESEKTATIOIC OF FORMULA (1.12) 


We shall now give a geometrical interpretation of ( 1 . 12 ) 


r 


2— IFo 

3 + IFo ' 


( 1 . 12 ) 


We plot s in a complex plane and also plot the points Wo and 
— IF„ on the real axis, z — IFo is the vector going from IFo to z ; 
z + IFfl is the vector going from — IF, to z, and the complex 
reflection coefficient is the quotient of these frivo vectors. Now in 
general the quotient of two vectors 


M, exp (jVi) 
.l/„ exp (jpj 


^ exp { 3(91 — Po) } 


is a vector with an absolute value equal to the ratio of the 
absolute values of the two original vectors and an argument equal 
to the differences of the arguments. 

Hence the absolute value of Pr/p. is the ratio of the two vectors, 
shown in Fig. 10, its argument, i. e., the phase angle b.v which Pr 
is ahead of p. , is the angle A in fig. 10 

The absolute value of r and therefore, the value of the ab- 
soiption coefficient ag — l — 1 r j- remains comstant along contours 
for which il/i/il/, — constant and from elementarj’ geometry we 
know that these contours are circles. For 100 % absorption il/j = 0, 
and the circle reduces to the point W<, (Fig. 11). On the other hand 


1 In accordance \vith usage in electrotcchnics wc denote the real part of s 
by R, its imaginaiy part by X, F being reserved for tlie admittance F = l/c. 
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■\ve also know from elementary geometrj that the locus of point z 
with the same value of A is also a circle, now going through the 
point 17o — 1^0 (Fig 11) Moreover the two sets of circles 



Modulus and arg^imcnt of tli® complex reflection coefficient r 
(eqoabon 1 12) m the z plane 



Fig 11 

The complex reflection coefficreat r in the impedance plane, 
phase Jump a and absorption coefficient <i«, the circle diagram 
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§ 7 

for constant a„ and constant A intersect orthogonallj-. This may 
be proved from the general property of transformations in the 
complex plane 'which are knoivn to be always conformal, i. e,, that 
angles of intersection are not altered. Now since in the r-plane 
the contonrs of constant A, which are the radii, intersect ortho- 
gonally the contonrs of constant [ r |, which are circles around the 
origin, the same must be true for the tn'o sets of contours of A 
constant and f r [ constant in the z-plane which is obtained by 
the transformation (1.12). The same result may be obtained in a 
more direct but more tiresome way by calculating the quantities 

and (~\ from (1.12), and observing that thev differ 

' OOt, 'a \ cA 'ag 

only by an imaginary factor, 'svhich means that the two directions 
«o-eon.stant and A-constant intersect orthogonallj’. 

The graphical representation of Pig. 11 is of verj- great value 
in two respects. If a theoretical expression for a as a function 
of frequenej* is obtained and the corresponding contour is plotted 
in Fig. 11. one sees at a glance the dependence of the absorption 
coefficient a„ and the phase jump A upon frequenej-. Converselj-, 
when, as is the case with most stationarj- wave methods, and A 
are measured, we have at our disposal a practical means of finding 
the corresponding value of a by a graphical method and without 
anj' further computation. 

« 7 WAVE niPEDAtICE OF IJIPER\*IOES MEDLV 'WITH 
INTERXAE FRICTION" 

We .shall now extend our considerations of § 3 to the ease in 
which internal damping occurs. As soon as this occurs the densitj- 
variation is no longer in phase with the pres.sui’e variation and 
the quantitj' dp/dp is no longer a pure real number’. Now bj- its 
definition the “corapre-ssion modulus” K is p dp/dp and bj- Kr 
we shall denote its real part. The dj-namieal value of K may 
then be mitten 

K = IU (1 -b jtan 5). (1.16) 

in which S, the loss angle, denotes the phase difference between p 
and the variation of p. This lass angle is alwaj's small and as 
our O'n'n mea.surements and those of other investigators have 
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slio\Mi, this 6 IS to Ji reasonable degree independent of u ,thns 
greatlj simplifying the problem Jir, hoireier, has proied to be 
dependent on u, at least at \erj Ion frequencies 
As before (§ 3) the equation of motion is 
8p dv 


restricting ourselves to a single travelling nave m the direction 
of the positi^ e x axis The equation of continuity is 
0U 1 0p lu> 

« sT 


from nhich we immediately. deri\e 

r = . ^y = }^W (J17) 

As S IS small ue may put approximately, 

m which case 

which means that the velocity of propagation amounts to 

c=i'^A,/f, 

so that 

r = and W=pc{l + ,±) (US) 

and the impedance of a layer of thickness I backed by a perfectly 
hard wall is 

2 = 17 eoth coth j j (^1 — J j (1 1^) 


(1 19) is a special form of the general expression (1 14) 

We see that the acoustical beluXMOur of impervious materials 
can be predicted from the xalues of their stiffness 7f, , their 
density p and their loss angle 8 
As a/j8 = S/2 IS not dependent on a, the coth comerges to a 
logarithmic spiral and as S is small the slope is also small Jlore 
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§ 7 

over its apex lies at the point pc (1 + j S/2), that is just a little 
above the real axis. Comparing this ease \nth that of undamped 
propagation ive .see that there are two corrections due to the 
occurrence of 5. The correction of Tl' is only of slight import- 
ance. the correction of y is more important as it gives rise to the 
creation of a spiral. The impedance of a layer of finite thickness 
is, therefore, pretty well illustrated by the curve of Fig. 8 and 



Impedance contour of a eample of cellnlar rubber, 
frequencies added along the contour in 100 Hz ' 


in our subsequent discussion we shall ignore the deviation of IF 
from the “undamped” value pc. 

An example is given in Fig. 12 where the impedance is shown 
for a sample of cellular rubber with completely closed cellular 
holes -. 


1 In order to avoid misunderstandiiig tvc tvant to make clear that e. g., in 
the point designated by 7 the frequency is 700 Hz (crcles/sec), but that 
a = Or • 700 radians/sec. 

- C. tv. Kosten and C. Zwikker, Pliysica, 8 (1941) 933. 
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"Without entering into the description of the method of measure 
ment something must be said about how Pig 12 and the following 
figures of impedance contours have been obtained Most of the 
experimental work has been carried out with an older type of 
interferometer with a rather small accuracy, especially as to the 
measurement of the phase jump A For better understanding the 
results have been smoothed without changmg the character of 
the contour The true measuring points may be found m the 
origmal publications The numbers added on the impedance curve 
are frequencies in 100 IIz 



Fig 13 

Absorptioa by cellular ruWier strong selectire absorption 

After each turn of the z contour it comes into the neighbourhood 
of the point Wj, , gmng rise to a high, absorption, coefficient ov er 
a relatively small range of frequencies, these absorption peaks 
being separated from eadh other by large ranges with onlj poor 
absorption (Fig 13, true values) 

Such resonance peaks occur in the neighbourhood of those fre 
quencies for which cothyZ is real Prom (115) we conclude that 
this is the case when 

Sm2;SZe=0, le 2pl = mr 

(resonance for n odd, antiresonance for n even) The resonance 
frequencies can be computed with the aid of jS = w/c, giving 
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S ‘ 

fls = ?iTc/2Z. The first resonance {ii = l) lies at v,—-cl2l, the 
second {n = Z) at a = Z-rcl2l. Expressed in the -wave length 
<\. = 2trc/(j, the resonances are found at: A. = 41, 41/3, 41/5, etc. 
The distribution of the velocity amplitude is illustrated in Fig. 14 
for the three lower resonances. The ratio between the frequencies 
of the first and the second resonance was well established in our 
measurements on cellular rubber (Fig. 13, curve for a layer of 
31.8 mm thickness). 

It should be mentioned here that the constant Kr cannot be 
deduced from its .static value, the dynamic value being several 
times greater. If no mechanical measurements are available Kr 



Velocity amplitude at the first three resonance frequencies 
may be computed from the first resonance frequency, which is: 



An easy way of determining the loss angle from the acoustical 
measurements is found by applying (1.15) to the first anti- 
rasonanee. Then sin 2 pi — 0; cos 2 pi =1 and 

coth yl = coth ttl =: 1 /al. 

Xow from (1.18) a = oS/2c and in the first anti-resonance 
1 = ttc/o, hence coth yl = 2/n-S. 

The corrasponding impedance is from (1.16) 

2.r.nr„.=pc2/:r8 (1.20) 


Socsd afc^ofbfo^ nati^rialj 
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and S can no^\ be computed from the measured \alue of 2 m 
the first anti resonance The \alues of S so determined agreed 
•well TMth those obtained in a purelj mechanical uaj 

“We maj also deduce a simple formula for 2 m the case of 
resonance Then sin2/31=0, cos2j8l= — 1, and from (115) 

coth 7!= tanhaf « al, 
and ■mth Z=7rc/2«, a = <fl8/2c we obtain 

2„, =/.ci-S/4 (121) 

\%hich formula maj also seme for the calculation of S An accurate 
^alue for Zn, may readilj be found from (1 13), since Oo is known 
wnth a high precision 

i S TflE GEVLBAL EQUATIONS GO\EIi\I'>,G THE TTAIE 
PROPAGATION IN POROUS MATERIAL ITITH 
RIGID FRAME 

"We shall discuss the beha\iour of porous materials for the 
time being with a iigid frame and direct our attention to the 
motion of the air 

The first difficulty is the definition of the velocity of the air 
m the material "NYhereas in an mfinetisimal volume only one 
pressure p exists the ^eloclt3 changes from point to point m the 
pores owing to their irregulantj Let us define the \elocity v 
in a certain direction therefore by the volume of air passing 
through a unit of surface perpendicular to that direction 1 e by 
the velocitj of specific % olume displacement Denoting the porositj 
(or ca^utj factor) bj h it inll be clear then, that for a gnen 
\elocit> gradient dv/dx a 1/A tunes larger 0p/3t will be found 
tb^n in free air Hence the equation of continuity maj be written 
m close analogy to the preceding cases as 

du A dp dp 
0x po dp 3f 

This time the quantitj p^dp/dp need not be real Leading vis 
cositj effects aside for the moment there are two limiting eases 
m which p and p iibrate in phase with each other In the first 
limit the transmission of heat to and from the solid frame is so 
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rapid that the enclosed air is kept at constant temperature (iso- 
thermal case, as was supposed in the old theory of X e w t o n) . 
In air at constant temperature a fixed functional relation between 
p and p exists (Boyle’s law). In the other limit the transmission 
of heat is so slow that the air vibrates adiabatically (Laplace’s 
■supposition). Then again, a fixed functional relation exi-sts be- 
tween p and p (Poisson’s law). In both cases dp/ dp is real. 
In general, however, there will be heat exchange with finite speed, 
so that, during the compression phase, the air is hotter than the 
solid frame, and during the expansion phase it is colder. Hence 
for the same value of p the density is smaller in the compression 
phase than in the expansion phase (Fig. 15) ; p and p are not in 



15 

Losses in tte- stiffness, hysteresis 

phase, dp! dp and also the modulus p^dp/dp are complex quantities. 

It seems imlikely at first sight that viscosity effects might 
influence the modulus (stiffness) of the air. since the modulus 
describes the behaviour of air during compression and expansion, 
which air may be supposed to be motionless as a whole. An 
important group of absorbent materials is formed by the acoustic 
porous plasters. Xow suppose that part of the holes in the 
materials can only be entered through one narrow pore. A time 
lag between p and the corresponding variation in p will now be 
found due to the ‘'RC-time” of these side holes, which can only 
be pumped up with difficulty. Therefore, phase difference betwc-en 
p and p cannot be put down with certainty either to defective 
heat exchange or to the latter side holes. Such holes may occur 
%'erv- ea-dly in the interior of the macroscopic globules of which a 
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plaster consists Of course this effect will be present to some 
degree m all practical matenafe ~ a\i 11 be denoted bj 7i, 

bj which, for a travelling wave, the equation of continiiitj again 
attains its universal form 


P, m which 

A h dp h 

The equation of motion for the enclosed air is 
0p L , 


(122) 


(123) 


As compared to the case of free air we observe three changes viz 
the appearance of the struciure constant K, the porosity ft, and 
the resisfance constant a The appearance of the latter will be 
clear from the fact that we have to tahe account of the visoositj 
For the steady flow state the term m dv/Bt cancels out so tJiat 
then, <r IS defined as the ratio of pressure gradient and velocitv 
of volume displacement 

There are several effects by which the densitj of the air seems 
to be enlarged, giving nse to the factor k/h > 1 This factor has 
been split up into two parts k and 1/ft tlie latter describing the 
rather elementary influence of the porositj, whereas k comes 
mainly from two effects, the structural properties of the material 
playing the moat important r61c, this is the reason for calling A 
the structure factor Supposing the porositj ft to be decreased 
without changing the structure, the pressure gradient, needed for 
the acceleration of the air, should obviouslj be taken 1/ft times 
greater in order to obtain the same velocity of volume displace 
ment v k contains a factor lying between 1 and Va because, as will 
be shown, in chapter II, owing to the internal fnction of the air 
a term between 0 and pg/3 is apparently added to the density p„ 
The influence of the structure may be elucidated by the following 
examples 

Imagine a structure like that m Fig 16a with pores placed 
mebned with respect to the maeroscopical pressure gradient 
Along these pores the pressure gradient is only cos 6 times that 
in pores having the same direction as the maeroscopical pressure 
gradient The acceleration of the enclosed air has m its turn 
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a component in the direction of the macroscopieal pressure- 
gradient, which again contains the factor cos 0. For this model, 
therefore, the structure factor k contains the factor l/cos^0. Should 
the orientation of the pores be distributed at random, then the 
resulting value of this factor amounts to 3. 

Again, consider a medium after the model of Pig. 16h, provided 
with lateral cavities located beside the air current in the main 
pores. Notwithstanding the application of a pressure-gradient the 
air in the lateral cavities remains largely at rest, thus acting as 
if it were heavier than the air in the stream lines and giving 
rise to a factor k greater than unity. In other words as far as 
this effect is concerned; since the pressure-gradient only succeeds 
in accelerating the air in the main pores, leaving the air in the 
side holes at rest, the porosity in the equation of motion should 



radius«=2r radius°>2r radius=tr 

a be 


Fig. 16 

To the concept of the structure factor. Three moilels haring the same 
porosity li and the same resistance r, but different structure factor fc 

be taken as the air in the main pores per unit of volume, i.e. as 
h/k. k, therefore, is eqital to the ratio of the total air contents 
to that of the main pores and is essentiall}" greater than unity. 

Yet another reason for the occurrence of k is the possible vibration 
of small parts of the solid skeleton, which also tends to enlarge 
the apparent inertia of the air. Although it cannot be denied 
that in principle this effect may occur, there are several reasons 
for warning the reader against over-estimating of this effect. In 
most cases (the more rigid absorbent materials) the effect is small 
indeed and the main contribution to k is due to the structure 
(e.g., k—3). In the case of very flexible materials the motion of 
the solid material maj* be considerable; in this case, however, a 
more adequate mathematical treatment can be given in which the 
effect of the internal forces between the air and the .stnicturc is 
coiTeetl}' taken into account (see Chapter III). As we restrict 
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ourselves m this section to ngid frames this effect mav be left 
out of consideratiorL 

The factor mentioned above between 1 and */, is a function of o 
Denoting it hr we may spbt up I into the factors and A,, 
where A, now is fully and eselusivclv determined bv the internal 
structure of the material 

It IS hardh pcesihle to calculate A beforehand from a description 
of the sample Nor is it possible to find it esperimentalh in am 
other war than by acoustic measurements This means that we 
cannot predict the acoustic behaviour of a porous material quanti 
tatively from such simple quantities as its porositv and air resLs 
tance According to experiments the structure factor A nonnalh 
lies between 3 and 7, but it may attain anv value exceeding 1 

Introducing the complex quantitv kpjh + afjo as the apparent 
densitv, p the equation of motion for a travelling wave is also 
brought into the universal form 

yp = ;<i)pi m which p = kpjh + (124) 

and as before we find for y aud TT 

y = Tr = ]'^Sp 

Putting P=y Po (1 — J tan S ) m which 

tanS; =<Th/«Apo and = •+• ;tanSir), 

we find for the impedance of a laver of thickness I 

} tan 5^) (1 + j tan St) coth juf 

I / ^ I / 1— jtanS^ 1 ^ 2a) 

hh. V 1 + ;tanSi J 

$ 9 THE BEJ5lSTO*CE C»^STA^T 9 

According to Poise uilles law the lelociti of lolume dis 
placement through a conduit with radius r in the direction of the 
positive X axis for steadv flow is 

~r* 3p 
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in which t] denotes the viscosity of air. If there are n parallel 
passages distributed over a cross section of unit area the amount 
is n times greater, and, as the porosity h equals riirr^, 


V 



dp 

dx 


In materials with pores inclined to the a;-axis or with pores with 
side holes (Pig. 16), v is decreased by the factor k, so that finally 


from which follows 


_ h r® dp 
k 8 1 ) dx’ 

k 8 7] 


(1.26) 


This is the value of a for a stationary current with a parabolic 
distribution of the velocity in the conduits. 

The expression (1.26) will also be used for materials not having 
identical cylindrical pores throughout the material. In such cases 
r has to be looked upon as a kind of mean value, defined by 
(1.26). 

In the case of rapid oscillations the character of the current 
changes. Only the outer layers of air are affected by the rigid 
walls; the inner parts of the pores, however, are swinging to 
and fro as a whole, the velocity being only controlled by its mass. 
In the limiting case of only a thin layer of air along the walls 
being influenced by friction, Helmholtz’ value for tr is to 
be applied ^ 

cr—^ ^ 2tju>p„. (1-27) 

As pointed out by Crandall (1. c.) the thickness of the layer 
which is hindered in its motion by the walls is of the order 
rj/oipa, and by a thin layer we mean a layer for which this 
value is small compared nith r. So the conditions depend upon 
whether the dimensionle.ss quantity p.— is great or 

small compared with unity {p » 1 Helmholtz, p « 1 Poiseiiille). 


II. B. Crandall, Theory of vibrating systems and sound, App. A., 
London 1927. 
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IS comparable with Reynolds’ number in hydrodynamics 
Both quantities are a measure for the relative importance of 
inertia forces as compared with viscons forces Their form is 
largely analogous (Re &= pvd/iy) 

Now, which of the two formulae (1 26) or (1 27) will have 
to be applied in the case of absorbent porous materials? Or will 
it be necessary to deal with an intermediate case giving rise to 
more complicated mathematics* The answer to these questions 
IS that all cases may occur on occasion 
The limiting eases will be dealt with in a correct way m II 
§ 4 An approximate (?) expression for <r in the intermediate 
frequency range will be given m II § 2 

We are fortimate in that the conditions for the two limiting 
cases, considered here, approximately eomcide with those for the 
two limiting cases of the compression modulus K mentioned in 
§ 8, and such that Poiseuille’s law applies to the nearlj 
isothermal case, while Helmholtz’ law holds in the nearly 
adiabatic case These coincidences are due to the close connection 
existing between the viscosity and the thermal conductivity of a 
gas, which, according to the kinetic theory, arc proportional to 
each other, and, m fact, are both linked up mth tlie mean free 
path of the molecules At the same moment, therefore, that the 
walls cease to influence the motion of the gas current in the 
middle of the conduit, they cease to exchange heat with it 



CHAPTEE II 


MOEE DETAILED THEOEETICAL CONSIDEEATIONS 
CONCEENING VISCOUS AND THEEMAL EFFECTS 
IN POEOUS MEDIA 

5 1 SOUND PHOPAUATION IN CYDINDBICAL TUBES AND PORES 

In §§ 8 and 9 of Chapter I several problems Avere touched upon 
that arise- Avhen studying the sound propagation through porous 
rigid media. The behaviour largely depends upon whether Ioav 
or high frequencies are involved. At extremely low frequencies 
the air resistance and the compression modulus have their static 
value, the real -density is increased by 33 % ; even this increased 
density is, however, negligible in comparison ivith the imaginary 
part of the complex density, ct/jo. At extremely high frequencies 
it is the density Avhich has its normal value; now the air resist- 
ance is "abnormal” (Helmholtz’ value) and the compression 
is adiabatic. 

We shall noAV proceed to the case of intermediate frequencies. 
A correct approach to the solution Avill be made taking Kirch- 
h off’s theory as a basis {§ 4). We then no longer consider the 
case of extremely low (high) frequencies, but that of simply low 
(high) frequencies, confining ourselves furthermore to the simple 
ca.se of the propagation in cylinders. In that case argument and 
magnitude of the complex quantities K and p, and therefore of 
IF and y, are slightly changed as compared Avith the corresponding 
expressions at extremely Ioav (high) frequencies, Avhile the mathe- 
matics are still relatively simple. In K i r c h h o f f ’s theory- 
A-iseous forces and thermal conduetiAdtj' are taken into accoiuit at 
the same time. It turns out that the correction factor of K only 
contains the thermal effect, that of p -only the auscous one. At 
first sight this is very plausible although AA-e do not see Iioav 
this can be proA-ed beforehand from the differential equations 
underlying Kirehh off’s theorj-. If, hoAA-ever, aa-c a.s.sume direetlj- 
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at the beginning that we maj omit the viscous effect ivhen the 
complex quantitj is to be calculated and on the other hand 
the thermal conductivity may be disregarded when calculating p, 
the derivation of K and p can be given comparatively simplj 
Moreover we are then in a position to derive expressions for K 
and p for all intermeduite frequencies and not only for low and 
high frequencies, from which of course, corresponding expressions 
for IV and y are easily found, olso for any frequency 

For convenience of the reader it seems most practical to gi\e 
the latter deri\ation first, accepting the necessary assumptions for 
the time being The derivation of p {§ 2) is essentially that of 
Crandall ^ repeated here for convenience The dernation 
of K will be given taking Kirchhoff's correct theory as a 
basis (§ 3) 

In order to avoid the assumption that K and p can be com 
puted separately , we shall then give the complete derivation m 
which viscous and thermal effects arc taken into account at the 
same tune (§ 4) We have to confine ourselves then, however, to 
low and high frequency resp This derivation is given for com 
pleteness but can be omitted by the reader, nho is not interested 
in all details 

$ 2 THE AIR DENSITY IN CYLINDERS DISREGARDING 
THERMAL EFFECTS 

The denv ation is analogous to that of Poise uilles 
law except for the introduction of an inertia term The 
total driving force on an annular nng of volume 2ifrdrdx 
IS — (dp/dx) dx2 TT r dr The actual velocity is a function of the 
radius r, giving nse to a viscous force on the annular nng 

Finally an inertia force has to be taken into account 
ju) p„ 2 TT r dr dx t 


II B Crandall Theory of t« 6 ra(tR <7 ssylems and swnd App A, 
London 1927 
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Equating the driving force to the viscous and inertia forces one 
obtains 


Cpfcx = 


( . 5? 2 

' ■'“po — r-ir 



I 


V. 


This raav be 'svritten 


— ^ — (JL 

Tj cx ^ cr- 



— ( 2 . 01 ) 


the solution being 

in svhieh -7^ is B e s s c I's function of zero order. At the cvlin- 
drical surface ( r = 1?) v must vanish, from which A can he deter- 
mined. Substituting this value of A in (2.02) gives 


1 cp ( .7„(7 t) ) 

3x ' J,(m^ 


(2.03) 


The mean velocity v over the cross-section is now found by 
replacing J^(kr) by its mean value over the cross-section 


B 

— ] J„ikr) 2^rdr. 
o 

This integral may be evaluated with the aid of the well-known 
property of Bessel’s function^ 

n <i 

J xJo{x)dz^xJi{x) I =aJ^{a) 

0 0 

in which -7. .stands for Bessel’s function of order unity. "U’e 
obtain 

= -hikR), (2.01) 

which substituted in (2.03) yields 


I E. Jaanfce acd F. Erode. FurktioMrtGfeln, p. 2ir,, Berlin ir‘33. 
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L i 1 ?_ J.(t-B) I 

jupo dx I A,i? ^ 

Defining the eomplex mean density by tlie equation 
_ dv 


2 05) 


we obtain the expression for p by comparison of this equation 
with (2 05) Replacing AJ? by its equualent fi — j p takes the 
torm 


P = 


po/! 


2 J, (y V—J) I 

\/- ]) I 


(2 OS) 


Approximations for low (high) frequencies may readilj be de 
Twed If ^ < 1 (i e , at low frequencies or in narrow pores) we 
may take three terms of the power senes for Jg and J, 

Jo — j) = 1 + V# J/*® — V et P-* 

Ji Cm V — j) = Vj II (1 + V» 3ix^ — V 10, p.*) 

leading to 

P = Vo Po 

or p = VaPo(l — 6j)j/wpo^®) 

1 St) 

or p = VoPo + — 

4/ J. 1 

or p = Va po + -— • 

3a> 

which IS in agreement with equation (126) 

If on the other hand ;4>10 (i e at high frequencies or m 
mde tubes), ne may put JJJ^ equal to j when the imaginary 
part of y — j IS taken positive, giving 

P Po (1 ^y J/m) I high frequencies 

or p = Po[l+ (1 — j) ) (208) 

which IS in agreement with equation (127) 

The contour of p as a function of p is given in Fig 17a At 
Ion frequencies (p small) the real part is equal to 4 pJ3 i e the 


M<1 

low frequencies 

(2 07) 
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density is increased by 33 % as a consequence of the parabolic 
velocity distribution over the cross-section. The imaginary part is 




Fig. 17 

Density and stiffness of the air in cylinders 


equal to the stationary value o-stai./jc). At high frequencies 
(/i great) the density deviates from the normal densitj'" po by a 
complex correction vrith an argument of minus 45° (see (2.08) ). 


J 3 THE COMPEESSIOH HODULHS OF AIB IN CYLINDEES 
DISB.BOAB.DING VTSCOHS EFFECTS 

If sound waves are travelling in the cylinder under consideration 
we can distinguish at any point of the tube the sound pressure, 
the three components of the vector velocity, and the excess 
temperature. All are functions of the position in the tube and 
the time. The dependence on time of all quantities will be 
supposed to be as exp jatt, so that d/dt can be replaced by ja. 

The experiment that i^dll virtually be carried out is the periodic 
compression in the axial direction of a cylindrical volume of air. 
Compression gives rise to generation of heat which, owing to the 
thermal conductivity of the air, tends to flow to the cylindrical 
boundary. Dissipation in the axial direction must of course not 
be taken into account, since such dissipation is of no importance 
in travelling waves, and this is the ease for which the compression 
modulus K is to be computed. 
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The excess temperature will of course depend upon the radius r 
and must be zero at the boundary r=Ji The same mil applj 
to the radial velooitj of the air The variation in densitj 11111 
also depend upon the radius, but will not vanish at the boundarj 
The pressure mil be assumed to be constant in the whole \olurae 
under consideration This is justified bj the fact that the radius 
of the cylinders under consideration is extremelj small compared 
to the wave length 

As to the \ariation of p m the axial direction the situation 
is still simpler dp/dx can be made zero by compressing the layer 
of air under consideration from either side exactly sjunmetncallj , 
iiz, by equal and opposite displacements of the flat terminal sur 
faces of the layer Speaking in terms of waves we can sa> that 
a thin layer of air is taken out of a standing wa\e at the position 
where p os a maximum Then dp/dx obviouslj is zero 
In the next 5 the correct derivation will be given in close 
connection with Rayleigh’s survey of Kirchh off’s theory*, 
making use of a notation that deviates as little as possible from 
Rayleigh’s The same notation will be used m this § 

Let 

p,Pjs= excess and equibbnum pressure resp , 
ip, pg — excess and equilibrium density resp, 

8 r, J’(,e= excess and equilibrium absolute temperature resp, 
v=; particle velocity (vector), 
s = condensation = Sp/pg , 

K =: Cp/Cp = ratio of specific heats, 

5' = a kind of relative excess temperature = ST/Ck — - llTo* 
Aft = thermal conductivity, 

Of = s Jond of reduced thermal conductivity 
A = Laplacian operator — + cjcl 

We then have the following four equations, which can be called 
the equation of continuity (2 09), the compression equation (210), 
the equation of motion (211) and the thermal equation (212) 


1 Rayleigh, Theory of sound II, 2na ed, p 310 a f Cambridge 1896 
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O 

£) 


div v = 

= — ds/dt 

(2.09) 

P = 

= Po { ^ + (k — 1) } 

(2.10) 

gradp = 

= Po 3 v/dt 

(2.11) 

de'/dt = 

- ds/dt -}- vAB' 

(2.12) 


(2.09) needs no explanation. It is the normal equation of conti- 
nuity if one bears in mind that s — tp/pa- (2.10) is a form of 
the law of Boyle-Gay-Lussae for small variations in Po, 
po and Tfj . (2.11) again is the normal equation and needs no expla- 
nation. Equation (2.12) expresses the fact that the temperature 
can rise by compression as well as by heat conduction. The com- 
pression term again follows from Boyle-Gay-Lussac, since 
for adiabatic variations 


and therefore 


T — p’'- ^ , constant 


ST/T,= iK — l) 8p/po. 


The conduction term is the normal one from the Fourier 
equation. 

Writing jta for d/dt and eliminating s between (2.10) and (2.12) 
we obtain 


U^i:^]0r + JL 

\ V ) Po 


Since p is independent of the space coordinates, the solution can 
be given at once 



Po 


+ AJ,r 




the second term being the solution of the reduced equation, the 
first a particular solution of the whole equation. 

A is found from the boundarj- condition 6' — 0 for r = R, gimng 

— jaK'/v ) 

/o (EK — jWv~)_ 

Substitution of this expression for B' in (2.10) gives 


K Po 
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P 

*Po 



JAr k-Wv) ~| 


The normal definition of the compression modulus K is 



_P 

s 


The ratio p/s is, howe\er, a function of r We, therefore, are 
need of a definition of the mean value of K over the erosssecho 
Denoting again mean values over the cross-section hr a b 
the required mean value of E must be computed according 
the definition 


p^Ks, 

for, then, the equation of contlmllt^ i> again brought into i 
usual form (lOS), since 

ds ___ 

“■ Tr“TF 


Notv p is constant, so we onlv have to compute s Averaging 
means averaging of the Bessel function in r This is nceon 
pbshed as in § 2 equation (2W) bj increasing the order of tl 
Bessel function br 1 , dmdiug it hy half the argument an 
finally replacing r by R 

If tve furthermore replace R I — j«c/v by its equivalent I — 
in which S=s)^ = 0 86 for air, we obtain 


E=kP^ 



Bp. k — 


-(k~1) 


Jo {BpV~j)J 


(2 13 


Comparison with the expression for p of § 2 (2 06) reveal* •'C^era 
analogous features Afost important is that the arguments of th 
Bessel functions are almost identical, since B « 1 This mean 
that the significance of the expressions “low”, “high’ and inter 
mediate frequencj ” is the same for the densitv p and the com 
pression modulus K 

Approximations for low (high) frequenciea for K are easih 
derived from (2 13) At low frequencies the argument of th 
Bessel functions is small, we then have 
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and 


2 


JABp.V-3) 
Ja (Bn V — j) 


1 + Vs {Bi^V-jy, 



iix < 1) (2.14) 


The modulus of K therefore does not change to the degree of 
approximation involved (see Fig. 17b). 

At high frequencies the argument of the Bessel functions is 
large, giving 

2 J, {B,xV=r^ 2 

if v~j is assumed to have a positive imaginary part, and 

A = kP„ [1 + 2 VVZTj (,, _ i)/bA . (2.15) 

y > 10 ) 


Therefore, both argument and modulus of E now change, K at- 
taining the limiting adiabatic value xpo in a direction that is 
inclined 45° to the real axis (see Fig. 17h). 

Comparing the results of p and K there turn out to exist very 
simple relations between the loss angles (the arguments) of both 
quantities, viz., (for air) 

tan Sp • tan S^- = 0.16 for /i < 1 

tan Sk = 0.46 tan Sp for p. > 10. 

"Without giving the derivation we shall write down the expression 
for E arrived at if p is not assumed to he independent of r, but 
still independent of x. The latter assumption is irrelevant, since E 
\vill be equal whether the volume under consideration is part of 
a travelling or of a standing wave, therefore, whether 0p/9x = O 
or 7 = 0 , 

E=^ Po X 


2 2 /,(/?]/ — A,) 


rV- 

— Aj JgiR y ^l) 

Ry'- 

-A, Ja{Ry—\^J 

ii 

2 /^ (/?]/_ A,) 

(l-^) 

\ Ja, ' 

2 J,(FK 

\ iw ' 

aK-A, J„(i?l/-A,) 

Ry^. Ja{RV 
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in which Xi and are roots of a quadratic equation from 
Kirchhoff’s theoi^y ^ and E is defined as p/7 

i 4 EXTEVSIOK OF KIRCHHJOTT THEOBY OF SOUND 
PROPAGATION IN (CYLINDRICAL TUBES AND PORES 

Taking into account at the same time viscosity and heat 
conduction, Kirehhoff developed hts bnllant theory describing 
the sound propagation, among other things, in tubes and pores 
of cylindrical shape 

The way m which we shall tackle the problem is that expres 
sions are derived for IV and y at low (high) frequencies, from 
which K and p then can be computed An expression for y at 
high frequencies is already given by Kirchhoff Also an 
expression for y at extremely low frequencies is available The 
corresponding expression for y for simplj low frequencies ^nll be 
denved m this § by giving the derivation to a higher approximation 
Expressions for TV will be denved from the definition TV«=p/w, 
m which p and u mean the average pressure and axial xelocity 
over the cross section Properly speaking p ought to be the mean 
force per unit surface, which is not exactly equal to the mean 
hydrostatic pressure The difference is the xx term of the fnction 
tensor (see Rayleigh*) which can be shown, however, to be 
negligible as compared with the hydrostatic pressure 
In order to avoid a repetition of Kirch hoff’s theory we 
shall take Rayleigh’s version of it as our starting point 
References to Rayleigh arc given with §, page and equation 
from the second edition (Cambridge 1896) Rayleigh’s 
notation will be slightly changed 

The four fundamental differential equations of the previous 
section remain valid except for the equation of motion (2 11) 
in which frictional terms enter 'The new equations of mofrmi 
are Rayleigh’s § 348 p 319 (2) 

du 1 3p 7 ^ _7_ 

dt po dx Po 3 po 


1 Rayleigh Theory of sowid H, 2nd ®d, p 321, equation 19, Cam 
bridge 1896 

2 Eaylcigh, Theory of sound II, 2ad ed , } 345, Cambridge 1896 
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and cyclic equations for v, to, y, z, when u, v, w denote the three 
components of the vector velocity. 

Elimination of all other variables than 6' giv^ an equation in 6', 
from which 


is a solution, if 


6' = A,Q, + A,Q, 
AQ^ = X,Q,, AQ,=:X,Q^, 


(2.16) 


■where X, and A, are the roots of a quadratic equation, -viz. 
Rayleigh’s § 348, p. 320 (16). Good approximations are 
^i=(i“/Co)^> X^ — jaK/v. The -viscosity turns out to be of much 
less importance for A, and A, than the heat conductivitj-, 

XoTV our object is supposed to ibe cylindrical in shape, further- 
more all quantities are supposed to vary as exp ( — yx), i. e., a 
travelling wave is supposed to exist. Finally a rotationally sym- 
metrical solution is the onlj' one that interests us, therefore the 
velocity at any point is given by the axial component ii and the 
radial component q. 

A, Aj, A^ being constants to be determined next, correct ex- 
pressions for u, q and 6' as functions of r are given (§ 350, 
p, 324, (8) (9) (10)), running in our notation: 


u — AQ + A^y (jo/Aj — v) Qi + A^y (/w/A^ — r) (2.17) 


jW rf - 


in which 


dQ 

dr 



e'^A,Q, + A,Q, 


(2.19) 


rf—y/pai 


Q^Jo{rV y^ — Hr}') 

Q^=^JArVy^ — X,), 

Q^.=JAr Vr-K)- 

These exact solutions have to fulfil the boundary* conditions at 
the cylindrical surface r — R: 


u = 0 \ 

q—O j for r — R. 

6' = 0 j 


( 2 . 20 ) 
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In these equations A, and arc unknoivn quantities still 
to be determined To peiroit a solution the determinant of the 
coefficients should be zero, whence at r^R 


r 


> (1 1 \ dinQ , \ 

xJ dr ‘^Va, '’f 


din Qi 




din Qf 


(2 21 ) 


This equation is Rayleigh’s § 350 p 324 (11) 

To a thxrd approximation we may write for small values of the 
argument of the Bessel functions (which applies to Q, Qi and 
Qi at low frequencies (§ 350, p 326)) 

dlfit/o {z)/dz= — (1 + Vs 2* + Vis 2*) (2 22) 


Introducing this into (2 21) taking for Q and Qj the expres 
sions gi\en above, an equation m y®, y* and y* is obtained from 
which y* must bo calculated In doing this only the leading 
terms should be retained This means that B/wave length 

and are considered as negligibly small quantities in com 

panson with unity All terms m y* cancel out among other terms 
and one finally retains 



(for low frequencies) 


This equation is in complete agreement with the equations for p 
(2 07) and K (2 14) calculated separately The relation should he 
y*=(j<u)*p/K and turns out to be complied with 
For high frequencies we may at once borrow a suitable expression 
from Rayleigh § 350, p 325 (15) 



which again complies wtb p and K, calculated separately and gnen 
by (2 08) and (2 15) rcspectivelv 

Once y has been taken as to fulfil the condition (2 21) the 
ratios A/Ai and A/ A, can be demed from equations (2 17) 
(2 18) and (2 19) Giving a suitable value to A one finds 
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• 4-1 = — VyQi 
^2 = l/r <?2 {aJB,)- 


(2.25) 


a,Ji stands for the complex argument of the quantities Q, at r — B. 

The argument of Q at a variable point r, is, therefore, a,r. 
Introducing (2.25) in (2.17) gives 

_ • f 1 1 ) ( j(a \ Qxia^r) fjo, \ ^,(02^-) 

^“Va, Q{aR) Va^ "J QxiccxR) \x, QAccJR)- 

The average value u is obtained by averaging Q, over the cross- 
section, which may be performed as in § 2, see (2.04). The in- 
tegrated Bessel function of order zero leads to a Bessel 
function of order unity. 

A corresponding equation for p is obtained by successively intro- 
ducing (2.25) in (2.19), calculating the condensation s from 

s = e'—iv/M AS', (2.12) 

where A is put equal to A^ and A^ respectively, substituting s 
and 6' in the expression for p (2.10) 

P = Po [s + (k — 1) O'], 

and averaging over the cross-section with the aid of (2.04). 

Now this value of p is not exactly equal to the mean force per 
unit surface, since it is only the hydrostatic part of it. Property 
speaking a frictional term corresponding to § 345, p. 313 (5) 
should be added. It can be showm, however, that this additional 
term can be neglected for the approximation required. "We then 
obtain 

PoC'f — /") Qi/ QA°iR) — Poi*^ — ’'-^2/ QJ 

-J«y(l/Ai — l/A2)9/Q(alS)+7(ifl>/^i — v)Q,/Qi(«ii?) — yU<»/X„ — v)QJQ 2 {ad 

(2.26) 

In order to derive approximate expressions for low and high 
frequencies suitable expressions for Qi/QxicciB), A^ , A, and y arc 
introduced. Retaining the leading terms, we obtain 
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W=,.c. [l + 1^) 1/j^] 

for high frequencies (2 27) 
11 = I^Va Po + y for W frequencies (2 28) 

From these equations and the corresponding equations for y 
(2 24) and (2 23) A and p are easily obtained with the aid of 
the wellknown equations (11?) The result is gi\en in Table 1 


TABLE 1 

y ir S A2CD p IS Gi LIKDERS DURn ED FKOif 
EIBCHHOFF S TIlEOBr 


1 

low frequencies ^ 

high frequencies 

H=]/ wpoAVi? < 1 

n=y mpoBVv > 10 

a _ 8 f- , __ 

pflfi* L 6 ij 



-(’"'- 1 ^) 1 /'^] 


II =p .^<.[3 + (/^ — 

1 


=p.[i+VsU-i)'^] ; 

lL = ,p.[l— 2(1/— -j^) 


[/ 

4/ . 1 877 

0— Apo+ 

p = p.[l + 2[/^^^] 
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and espenment are Jn good agreement. Experiments hy the 
aiitJiors and Mr Tan den Ejk*, aEthowgh less acenrate than 
those of Kaje and Shernatt* ga\e the same result. 

} 5 TIIE CONSEQUENCES OF KIRCHHOFF S THEOEl FOR THE 
WA\E PROPAGATION IN POROUS MEDIA 

'We maj raise tie qiwstion, whether the results of Table 1 
can be applied directlj or after small alteration to porous media 
"We shall see that this is onli the case for media inth j>ores 
of ejlindrieal shape of the same diameter throughout the i^bole 
material, a rather hj-pothetieal case Ordinary materials, especially 
those containing pores ^ aiding greatly m cross-section along their 
length, gne rise to great difficulties For materials \\hich only 
contain pores of one diameter, the equations of Table 1 arc ^alld, 
proMded. K is divided by h, and p is multiplied bj K/h The loss 
angles 5; and will obviously remain unchanged, because in 
ever} pore the propagation is exactly as given m Table 1 there 
fore the ratio of resistance forces to inertia or stiffness forces, 

1 e, tan S, and tan S/c rcspcctnel}, will be unchanged for the 
whole as well 

■NVhetlier the propagation mechanism is of the PoiseuiIIo 
Kirchhoff tjiie or of tic Helmholtz K i r e h h o f f 
t}T>c, ma} readilj be found bj calculating p from 

p s= — ySakpJhotxki , 

see (126) The quantities K, h and have to be inserted, 

h and <r*u, are taken from measurements, whereas J is found from 
the resonance frequency of an absorbing laaer or put equal to 3 
as a rougli estimate If/i <1, Poise mile Kirchhoff 
holds, whereas /i>10 certainly indicates that the Helmholtz 
Kirchhoff mechanism is valid 

The theory fails in the case of media containing pores the 
cross sections of which \ary greatly along their length In order 
to make this clear, consider the extreme case of a material in 


iC Zwiikcr, J tna dea E^ji, awl C Kosten Fhr/s^ca 10 
(1943) 239 

2 and 6 her rat t, Proc Eoy iSoc A, 141 (1*^33) 123 
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whicli the air content is present in rather large holes, mutually 
connected by very narrow pores. The air resistance is mainly 
due to these latter pores and, owing to their small cross-section, 
may essentially be of the Poiseuille character. The behaviour 
of the compression modulus is governed by the velocity of heat 
exchange of the air, and, because the greater part of the total 
air content is found in the rather large holes, this heat exchange 
may occur so slowly that Helmholtz — Kirchhoff applies. 
Therefore, it may happen that 1 < p < 10, and that nevertheless 
p complies with the first column of Table 1, whereas K is found 
from the second column of Taible 1; in other words: the motion 
might take place almost isothermally, the compression almost adia- 
baticalljn The viscosity losses are governed by the conditions in 
the narrowest parts of tlie pores, the heat losses by those in the 
■widest parts. 

The following general trend, therefore, may be expected: 

a. If /t = 1, as calculated from vstat. , Poiseuille 's law for 
tan 5f vnU hold good, but tan may be considerably greater than 
0.028 p?. This 'n'ill be shown to be of the utmost importance for 
the absorption at low frequencies. 

b. If p=:2, -w^hich may be about the limit for the applicability 
of Poiseuille ’s law, tan may already have readied values in 
the neighbourhood of its maximum (0.1 — 0.2), 

c. If p>20, Helmholtz-Kirchhoff’s theory might be 
applicable both as to tan Sk and tan Sf , i. e., column 2 of Table 1 
might hold good (after the changes K-^liK and p-^^hp/k. 

A rough estimate of p. at 500 Hz for different materials is 

acoustic plasters ; <r = 5 • 10® — 10® k/h 1 = 4 p = 4.5 — 1 

hair felt : cr = 5-10‘ — 10® k/Ji = 3 p= 1 — 0.3 

wood fibre plates; cr = 5-10® — 10’’ k/h — S p = 0.3 — 0.1. 

This means that eases a and h will indeed occur, that is that the 
motion -nnll frequently be isothermal, while tan Sx be greater 
than according to the first column of Table 1. 

As we shall see, the complex quantities p and K, and therefore, 
also their loss angles, may be calculated from impedance measure- 
ments, so a check of our theoretical considerations is possible. A 
question of more or less academic importance may be answered 
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too in tlie mean time In the eqoation of motion the structure 
factor k was necessarily introduced, it seemed impossible, however, 
to check whether k originated m oblique pores or in side holes 
Now an idea about the relative importance of both effects may he 
expected to be gamed from measurements of A' If oblique pores 
are preponderant. Table 1 should apply, if side holes are more 
important, the required /i value for K will be greater than for p 

The same results as deserrhed tn thts section would haie been 
found if the heat conductivity \h were supposed to be several 
times smaller than the actual one This causes a shift of 7f along 
its complex contour (Fig 171>) m the direction of the adiabatical 
value, giving nse to greater Si values at low frequencies, smaller 
at high frequencies The ratio of the apparent heat conduetivitj 
Xapo and the normal value must be looked upon as a new 
material "constant" It will depend upon the frequency and the 
structure and must be obtained from experiment 

$ 6 DISCUSSION OP THE THLORCTIOAL RESULTS OF 
KOERLSGA, KBONIO, AND SMIT 

The behaviour of porous materials was described in the beginning 
by the material constants h, k and a For matenals with ideally 
eyhndncal pores these constants seemed to be sufficient Consider 
ing practical matenals, however, a further “constant” had to be 
introduced, the ratio A^pp/A* (< 1) Although the need for the 
introduction of the latter constant is obvious and the constant has 
a physical meaning nhidi is easily grasped, it is clear that this 
way of attack ceases to be elegant, the more so since it is rather 
difficult to give an even rough estimate of the ratio Aapp/Aj 
A mathematically more sound way of attack was proposed bi 
K'ornnga, ATronig, and Smit', aiYAoagA fAfs method hs& 
the disadvantage of only being applicable m geometnea]]} well 
defined and simple cases, even then being quite labonous Thej 
considered the propagation of sound waves in a medium consist 
mg of spheres, situated on a cubic lattice, taking into account 
viscosity and heat effects, thus applying Kirchhoff’s method 


1 J Eorringa, R Kronig, and A Smt, Phystca, 11 (1945) 209 
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(who considered viseositv and heat effects in tnbes) to a problem 
treated by Rayleigh, viz., the propagation in a medium of 
spheres, neglecting damping effects h 

Their results are given in a few simple graphs, gmng y (dis- 
regarding a factor ja/c^) and TV/p^c^ as a function of frequency 
for different radii of the spheres and porosities. Their results 
es-sentially corroborate our views and give us an idea of the magni- 
tude of the ratio -Wpp.A^ for the media under consideration. 


TABLE .3 

I- Ayi) Aapp/;ft FEOir KOERrN'GA, KEOXIG, AjSD SUIT’S 
THEOET FOE PACKFE SPHERES 


7i 

diam. 

di.st. 

k 

^prj,/Xft 

0.8.5 

1 

0,66 

\ 

1 1.1 

0.45 

0 75 

0.78 

1.1 

0.45 

0.6.5 

0 87 

1.2 

0.45 

0.55 ! 

0.95 ' 

1.2 

0.45 

0.48 i 

1.00 ! 

1.3 

0.45 


Table 3 gives a survey of our interpretation of their results. 
The arguments of p and K, in the frequency range for which 
K o r r i n g a et ah gave numerical results, are rather small in- 
dicating that we are working in the Helmholtz’ region. Ac- 
cording to our scheme the wave impedance at high frequencies 
should be \-^ k/h times , and, since h is known, k may be com- 
puted. One really finds figures greater than unity of the expected 
magnitude. From y and TT the complex K and p can be computed 
and. of course, the ratio tan SAtan 5p interests us areatly; Kirch- 
hoff predicts in the ease of tubes 0.46 for this ratio, which 
according to our conception should be 0.46 

The results, obtained in the paper under discussion lead to 


I Rayleigh, PTilI. iXan^ 3i (IS92) iSl-, Scv^nfic papers, Vol. i 
fI903) p, 15. 
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ratio’s X»pp/Xft of about 04 to 0 5, very accei)table values for this 
case Because ive are, now, in the Helmholtz’ region, this 
means that Kirchhoff’s correction is two times smaller than 
in tubes On the other hand it completely corroborates our Mews 
and justifies the suggestion that the contribution of heat exchange 
losses in the more important ease of Poiseuille flow is in- 
creased by the irregularity in cross-section of the pores 

} 7 POROUS MATERIALS WITH HIGH AIR RESISTANCE (/^<2) 

We shall now proceed with the consideration of the absorption 
characteristics of media m the extreme cases of small and large /», 
and, in order to give some idea of the influence of the above 
mentioned ratio of eondnctivities, this ratio mil be inserted into 
the equations We must bear in mind, then, that the quantity v 
IS proportional to the thermal conductivity 
The general expression (114) for the impedance of a lajer of 
thickness I 

z = coth y p/K 

may be written for small \alues of the aigument of the coth (i e , 
far below resonance) as 

z^K/]>ot + (229) 

For, if 25 « 1, 

coth 25=-^ 4- 4- (2 30) 

X 6 

(2 29) is the impedance of a mass Zp/3 mounted on a stiffness A/Z 
The effective mass of the lajer per unit surface, therefore, turns 
out to be ^/a of the total mass/m® 

If /t < 1, the values of the first column of table 1 for K and /> 
may he /nserfed h am? trii’iafe jidds m ths ess& J'J nor 

_£ ^ \ 1 4- ( ^ ^ i , Po I 

PijCg 3 pgCg t \ al / ktff 20 h* I jtuZppCoA 

^ jaZ jpQ _ R~3^ (2 31) 

3 P<J®0 

in which o denotes the steady flow air resistance It should be 
borne in mind that Cg stands for the sound velocity in free an, 
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c^ — \^xpa/pj. Far beIo~ resonance the mass reaetanee term in 
X may be dropped. The layer behaves like a spring -srith losses. 
We Xen have 

X —pa/mlll (2.32) 

X is neither dependent on k. nor on c. and indeed all esperimentaJ 
contours of a show the same trend at low frequencies (Fig. IS). 



SitibOtLt-c cvrvBS of measured covtovrs. freqcT-EciFS ia ICO Hz 


The absorption coefneient in this region is approximately pro- 
portional to R. 

G,=4p,c,B/Xh (2.33) 

Therefore, F is of great importance for the absorption coeffi- 
cient. E passes through a minimnm if v varies. As a rule 
(since X, 20 7>' ss 1) = 2 p,c, . corresponding with a 
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\eilical line for 2 at a dislaj^ce 2poCo iTOia the imaginarj axis 
Under extreme conditions (A,pp =Xit,l = l,?i = l) ^min w V 2 po^o 

Inserting X from (2 32) and taking J?„ia =poCo eqiiation (2 33) 
gues 

Ojaln ■ (2 icut 

We should hke to emphasize that (2 29) and (2 31) onlj hold 

1) if P ois eui He’s law holds (1 e , about w < 7icr/2 i.po), and 

2) if the expansion of coth into tno terras is alloned (1 e, 
about tu < pJahP) 

The lowest of these limiting frequencies forms, of course, the 
limit of applicability of the equation For materials with low 
a (plasters) the first condition is the more severe and traces a 
limit at 50 — 1000 Hz On the other hand, for materials ^vlth 
high o, the first condition is satisfied at all audiofrequencies 
but the second is passed already at rather low frequencies (10 Hz 
for very high air resistance) An easj control is obtained m 
the complex plane According to the second condition, all 
impedances to which equation (229) applies should have an 
imaginary part about 15 times greater than the real one, le, the 
argument should be between — 55® and — 90® 

The correction m U betnecn the brackets of (231) due to 
Kirchhoff’s effect in is only of importance if 
To give some idea of the magnitude of this factor m practice, it 
IS given m Table 4 for a few materials It is of great importance 
uhen the air resistance is low' 

TABXjE 4 

BOUGH L 6 TIILVTE OF MROIIHOFF S CORRECTION 
IN THE POISEUILLE RANGE <1 — 0 02) 


matensl j 

corr factor 

plasters I 

25—100 

felt 1 

1—2 

ivood fibre | 

1 


C 11 IC 0 s 


Appl M Sfseoreb, B1 (19^9) 241 
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For large values of a we know' that coth converges tow'ards 1, 
and, p being almost purely imaginary, a frequency contour for 
coth as in Fig. 7 is to be expected. Putting cothyl=l, 2 equals 
the wave impedance, approximately 

z = I^Tpr V'' 1 — i tan Sp . 

Now — u) is a rectangular hyperbola (Fig. 4 low'er branch), 
which explains the form of the upper part of the contour shown 
in Fig. 18 curve c. This is the case “with materials wdth extremely 
high cr. At low frequencies their absorption is relatively .high, but 
they do not reach really high values of , even at high frequencies, 
because their wave impedance remains appreciable. Materials with 
smaller a behave like Fig. 18 cun-e a and b. At high frequencies 
Poiseuille no longer holds, coth shows maxima and minima 
(Pigs. 5, 8, 9), giving rise to a spiral-like impedance contour and 
an absorption characteristic with pronounced maxima and minima. 
Pig. 18 curve a holds approximatelj' for 20 mm hair felt, c for 
wood fibre plates. 

§ 8 POROUS MATERIALS WITH LOW AIR KESISTANOE (/^>20) 

As in § 7 equation (1.14) holds for the impedance of a layer I 
on rigid background, and, again (2.29) is its approximation for 
low frequencies. Furthermore, because the modulus of p is much 
smaller than in the case of p- <2, tire approximation is valid up 
to higher frequencies this time. A further approximation is allorv- 
ed, since both K and p have arguments, 8^- and Sf, of the same 
order, whereas the modulus of the mass reactance at low' frequencies 
is negligible as compared w'ith that of K. Dropping, therefore, the 
term in p gives for low frequencies 



leading to an absorption coefficient proportional to w. Of 
course this dependence cannot remain valid for extremelj' low 
frequencj', since Helmholtz’s value must then be replaced by 
Poiseuille ’s. This time the real part of z varies as 
giving rise to the euiwed impedance contour often found for 
materials w'ith verj' low' resistance at low frequencies. 
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For the higher frequencies we obtain, when only the leading 
terms are retained, and, moreover the irrelevant small loss angle 
m the wave impedance 13 dropped 

{}^lV^7e^{l—3 (1 + 046 Kwws,/2}1, 

which IS represented in the complex plane by a spiral with a small 
slope of less than 073 Sf. Now since 8f in this case is proportional 
to <i>" the slope steadily decreases, although there is no asj mptotie 
circle os in fig 9 
For very high w z approaches 

Z=VKrPr'= ~VkpoCaX{2 to 3) p^Cg, 
n 

therefore, the absorption characteristic will have very pronounced 
maxima and minima “around” a mean value, which amounts to 
about 80 % almost independently of the material under con 
sideration 

Ar can be determined from resonance frequencies, for which 
I iKF/<!.=»/4, 3Tr/4, etc 

or from anti resonance frequencies, for which 
«»onr«» fc/C(, = ir/2, IT, etc 


j g A METHOD OP VERIFICATION OF THEORETICAL RESULTS 

The correctness of theoretical results can be verified by evperi 
mentally measuring the impedance of a layer of thickness I If 
the wave impedance and propagation exponent, used for the theo 
retical expression of z, were correct, good agreement vith experi 
ment should be found Discrepancies will generally be found , it, 
then, turns out to be rather difficult to find the reason for these 
discrepancies The best wa> of getting a good insight into the 
material and its behaviour is, however, to compute the complex 
quantities K and p from the experiments The reason for dis 
crepancies between the actually found K and p on the one hand 
and any couple of theoretical expressions for K and p on the other, 
can be found much more easilj than when comparing experimental 
and theoretical results of the impedance z of a layer of thickness I, 
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§ 9 

becaiise the mathematical distance between the fundamental as- 
sumptions and E and p is much shorter than that between 
these assumptions and z. 

This way of attack, which is common practice in the electrical 
field, was perhaps first applied to the propagation of sound waves 
in absorbing media by W ii s t A possible scheme for the com- 
putation of K and p from adequate measurements is the follounng 

3 , = W eoth yl W = V'^ K = jo^W/y 

Sj = iT tanh yl coth yl—V zjz^ p = Wy/ja 

Zj and Zj are the measured impedances for the layer backed by 
a rigid wall and a colirmn of air of one quarter of a wave length 


- Re fe/So) 



Fig. 19 

Jleasuremcnts of the complex stiffness and density 


respeetivelj' (see equations (1.05) and (1.06) ). 17 is now calcul- 
ated from the product ZiZ,, 7 from the ratio zjz„. The ne.xt step 
is the computation of E and p fi’om the well-Itnoum relations 
between 17, y and E, p ( (1.17) ). 

Although this way of attack is doubtless the most promising one, 
only a few measurements of this kind have been carried out. 

In Pig. 19 a geneia.! survey of results is given comprising 
those of IV ii s t ( W) , P e r r e r o and Saecrdote^ (P, S) 


^ II. tVu.st, Hochfrcqncncicchn. n. FlcVirooVusW:, 44 (1034) 73. 

- 11. Ferrcro and 6. Sacerdote, ATioro Cimenio, 4 (1947) .3. 

Sound absorbing materials 
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and of C AV K os ten, M L Kasteleyn, and J V’an 
den Eyk^ (K, K, E) Wlist, Ferrer o, and Saccr 
dote described results for felt only, whereas the results of 
the last group of investigators (K o s t e n et al ) also include 
measurements made on acoustic plasters These latter results 
especially lack accuracy probably due to the fact that the samples 
of plaster were insufficiently homogeneous 
All p curves shou the upward trend ivith frequency u hich was 
to be expected (see Fig 17a) Except for the measurements of 
F e r r e r o and Sacerdote, the real parts of p are nearly 
constant indicating that the structure factor is constant With 
the results of Ferrero and Sacerdote, the structure factor 
obviously decreases with increasing frequency The structure 
factor is of the order 2 for all felt samples 
The compression modulus K of felt lies somewhere in between 
the isothermal and the adiabatic value In the results of W u s t 
an unmistaliable progression towards the adiabatic value is to be 
noticed The loss angle of Ferrero and Sacerdote has 
about the maximum value m agreement with their 5f 
The general tendency of W u s t s results is the same as that 
of the results of Kosten et al for felt His resistance seems 
to be a little smaller giving rise to the theoretically expected 
frequency dependence of K The results of Ferrero and 
Sacerdote, although extremely accurate deviate much more 
in general trend Moreover it seems rather difficult to give a 
theoretical explanation of them Since tan 8 ^ is of the order 
unity K might be expected to be rather strongly dependent on 
the frequency whach is not the case 

The behaviour of the plasters 1 and 2 is much more complicated 
The real part of /> is hpjh, yielding k/h values of 12 5 and 21 
resp The real port of X is proportional to 1/h Assuming the 
variation to take place about midway isothermal to adiabatic this 
jields li- values of 0 34 and 03 rtsp The structure factor there 
fore, IS found to be approximately 42 and 6 3 resp, which are 
reasonable values The low ft-values are not in agreement with the 
statically measured ones, which are about 0 7 This will be due 


» C W Kosten, M. L. Kasteleyn and J Van den Eijk, un 
published. 
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to the fact that in the statical measurements all cavities are taken 
into account, also those, which, in the dynamical situation, do not 
take any part in the mechanism. This fact also explains other 
features of the curves. In the ease of plaster 2 8 a- is much greater 
than its theoretical maximum. Moreover the variation of the real 
part of K is greater than could he explained from the mere 
transition from iTieoih. to jTadiab. (1:1.4). The effective porosity 
obviously decreases more or less with frequency, which decrease 
is accompanied by extra losses in K. 

The question arises whether the measured loss angles are in 
agreement with those calculated from measurements on the steady 
flow resistance o-stat. • This turns out to be the case for 8,, ; the 
calculated values for 8 a are, however, smaller than the measured 
ones. Since 8p for the plasters is rather small, the relations be- 
tween 8p and o-etat. for the Helmholtz-Kirchhoff case 
wll approximately apply 

tan 8p ■ — ^pr 

(see (1.26) and (2.08), bearing in mind that pr = hpjh) 
and furthermore Sa can be calculated from 

tan 8 a = 0.46 tan 8 p . 

The results are given in Table 5, showing good agreement for 8p. 


TABUD 5 

CALCULATED AIsVD ilEASURED LOSS AKGLES AT 1500 Ha 
FOB PLASTERS 





tan 


tan 8 A 


O^stat. 

Pr 

calc. 

meas. 

calc. 

meas. 

1 

5500 

12 

0.11 

0.1 

0.05 

0.15 

2 

34000 

24 

0.19 

0.15 

0.1 

0.5 














CHAPTER III 


THEORY OF SOUND ABSORPTION BY IIO’MOGENEOLS 
POROUS AND ELASTIC LAYERS 

f 1 POBOCS SIATEBIALS WITH ELiASTIC FRAME* 

If the possibility of nbralion of the solid material, dae to 
Its finite stiffness, is taken into aceount, the mathematical problem 
of the coupled vibrations of frame and air arises 

This problem will be dealt with in this and the following 
sections In order to distinguish clearly between the quantities 
of frame and air respectively all quantities relating to the 
frame are wntten with an index 1, those relating to the oir with 
the index 2 Let 

Pi s= force actmg on the frame per unit ares of cros!> section 
of the sample, 

ti = mean relocitv of the solid material at x, 

Pi = densitv of the solid material given by the mass of solid 
material per unit volume of porous material, 

Ki = specific stiffness of the frame to be taken complex if ne- 
cessary, 

Pj = excess force acting on the air per unit area of cross-section 
of the sample, the equibbnnm force hai> been subtracted 
therefore it is called the excess force , is eloselv related 
to the mean sound pressure at x, which is 1/h times greater, 
if 7i denotes the jHirositv, 

V, = mean velocity of the air at x, this quantitv is 1/ft times 
greater than the velocitv of volume displacement, used up 
to now, 

p, = density of the air as given bv the mass of the air content 
per unit volume of porous matenaJ, it is 7i tiroes the den 
sity of free air. 


I C W Ko«ten and C Zwife&er, Phy^ea, S (Wl) 06s, C W 
Kooteu, Thesis Delft I'll’ 
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Kz — normal modulus of air, if nec&ssaiy to be taken complex 
in the ease of losses being present (thermal or 'viscous), 
Po = total equiKbrimn pressxire of the air, 

P ~ total pressure of the air during the •vibration, 
s — coupling coefficient, to be defined next. 

The equations of motion are 


3x 

3P2 

0X 


0l/'l 

3f 

-f s (Ui- 

-f,) 

(3.01) 

'W 

- + s (rq- 

-rq). 

(3.02) 


The first term of the right hand side of the fiist equation is the 
ordinal^,’ mass reaction force; the second term stands for the force 
per unit volume of porous material exerted by the air on the frame 
if Uj ^ fj . The definition of s is obviously the force per unit 
volume imr unit difference between the velocities of frame and 
air. Oving to the law of action and reaction the same term also 
occurs in the equation of motion of the air, however, of course, 
\vith the opposite sign. 

The coupling coefficient s can easily be connected ■^^•ith well- 
kno’wn quantities. Supposing i\ to become zero 'V'vithout changing 
the structure, (3.02) must as.sume the u.sual from (1.23) from 
which follows 


s — jap^ih — 1) 4- 7iV (p2 = ^'Po) (3.03) 

This means that coupling even occurs in the case where <7 = 0. 
The equations of continuity will be shown to be 


dx h~ dt 


dt 



+ (1 — h) {K._ — P„) 



(3.01) 

(3.05) 


The first equation expresses the fact that the force in the frame 
will increase by a concentration of solid matter as well as of air. 
The term containing /ij is the normal contribution which would be 
present in vacuum, arising from the stiffness of the frame. The 
nece.s.sity of the second term can easily be demonstrated as follows. 

Suppose dvjcx to be zero in a special point of the material. 
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If, furthermore, it is supposed that the frame does not deform if 
submitted to a hydrostatic pressure (e g, porous sample in vacuo 
as compared to porous sample m free air), then it is clear that pi 
must increase together with p,, if the frame does not deform 
(3Vi/3a;c=0) Indeed Pi/(1 — h) and p^/h are the pressures in 
solid material and air resp, which should \ary bj equal amounts 
if dvjdx be zero 

Equation (3 05) states that the air pressure will increase bj a 
concentration of air as well as of solid matter, the first term of 
the right hand side relating to the influence of air concentration, 
the second to that of solid material (3 05) maj be deru ed as 
follows 

Pi — hP — \Po=aetual force — equilibrium force, (306) 
Pj = ;ip, (3 07) 

if p denotes the variable density of air 

P and p are connected through the compression modulus of air 
(complex if necessarj) 

dP dp 

/fi ■” Po ' 

which equation, with the aid of (3 06) and (3 07) jields 

h dt k dt p, \ dt ha#' 

Now 3h/3# can bo replaced by (1 — h) Zv^/dx and Sps/S# bj 
— Pi ^Vi/dx, which after substitution j lelds the equation of con 
tinuity of the air (3 05) 

The physical meaning of (3 05) can be elucidated by the fol 
lowing qualitative consideration If in an infinitesimal volume 
solid material is introduced owing to 0t»,/3x not being zero, the 
air pressure wll increase The oir force p* will not be increased 
if tbe compression takes place isothermally, for, then, the pressure 
increases m the same ratio as the porosity decreases Indeed the 
second term of the right hand side of (3 05) cancels out in this 
case, owing to if tbe compression is isothermal This ivill 

applj at low frequencies, thus materially simplifying the considera 
tions As we do not confine oiirsdaes to this case the complete 
expression will be used 
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§ 2 TEAVELLIXG WAVES IX ELASTIC POEOES IIEDLV 


In order to find expressions for a travelling sine-shaped vave 
we substitute — y for d/dx and ja for d/dt in equations (3.01), 
(3.02), (3.01) and (3.05). ja is a given number, y has to 
be solved, giving us damping and velocity of propagation of 
the supposed wave. 

Indeed the quoted equations of motion and continuity now form 
a set of four homogeneous linear equations in , p^, , and m , 

which only permits a solution, if the determinant of the coeffi- 
cients is zero, thus giving an equation for determining y. If 
furthermore we write 

r = y/i<j (—l/complex sound velocity), ) 

S — s/ja (=generah'zed air density). ) (3.08) 

this determinant is 



P1P2 ■ (P. ^ 

^ JiK^K, 

Since only even powers of P enter into (3.10). the four roots 
of the equation may be written ± T' and ±: T". For the moment 
the negative roots may be omitted as they are related to waves 
in the negative z-direction which do not interest us for the general 
considerations. 

VTe see that at a given frequency two waves are possible in 
the positive x-direetion. The amplitude of each wave will be 
dependent on the boundary conditions, We should like to 
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emphasjze that these two na^es Tnay not be considered as the 
\va\e in air and that m the frame reap The Tvase tjT)e relating 
to r' takes place both m the air and in the frame, and so does 
that of r" 

A fanly good analogy is to be found m optics, nbere double 
refracting media e\ist One single irave impmgmg on such a 
medium gives ti\o refracted wa\es In the acoustical case it seems 
ho\\ever rather inadequate to distmgiush between an ordinary 
and an e'^traordinary beam 

In special eases P' or r'' maj assume a numerical lalue uhich 
reminds one strongly of the propagation constant of a wa\c in 
the air or frame separately This will be the ease if the couplmg 
IS i^eak Hor\e\er such a ware m air or frame is alwajs 
accompanied by a naie m frame or air of the same r 

For r equal to either P' or r" the equations m Pi , p* , and t , , 
of Vihich the coefficients maj be taken from (3 09), maj be used 
for computing all ratio’s between p,, p, , t?j and e g, pjvi 
and Pj/i’a, which bare the character of partial wave impedances 
Since r may take one of the two \alues r' and all ratio's 
ha^e also t^vo values We can distinguish, therefore, four wave 
impedances 

lV,"=Cpi/r,)" i 

w\ = (pjv^y ir "= (p,M)'^ ) 

The ratio’s (tj/ti)' and {vjv,)" decide whether the role of the 
air or the frame is preponderant in the wav e type under considera 
tion. The wave to which the smaller ratio belongs, maj be called 
the disturbed frame wave, the other the disturbed air wave, since 
tg is the velocitj of the air 

The ratio of aii} two quantities may denied from the four 
fundamental equations The following general expressions for 
11^2 und t,/i\ may easilj be derived 

^ _ ftPg + (ft + pi ) S SKtT^h / (1 — h) _ 

^ t'. {fh-hS/(l~k)}T 

(l^Ji)SKg+hK,KX^—PgKg~SKt—(l—mp2-^S) (Aa— Po)/A 

7ih.j r* — ps — s 
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ry P 2 P 1 P 2 + (ft + ft) ^ P^K^T'—SEJ:' j 

" {—K,T^^Sni + p,)T ~l 

(p^-hS) (1 — h) (K, — P„) r + iiK,sr I 
(1 — 7t) (/£,~PJ r- + S' ’ / 

— hK,T^ + hp^ + S (1 — h) {K, — P„)V- + S 
'Uj (1 — h) p, + S — hKj!^ + p2 + S 


(3.13) 


(3.14) 


All special eases are comprised in the general one. E. g., put- 
ting TCj ~ CO will yield the ease of acoustic plasters, S—vj applies 
to materials with extremely high air resistance or to normal 
materials at extremely low frequencies (see equations (3.03) and 
(3.08) ). At low frequencies assume its isothermal value Pq . 

At high frequencies »S' is equal to {h — l)p 2 and a reasonable 
assumption for is 1.4 P^ (adiabatic value). 

A few results for special cases may be taken from Table 6. 


Thro interesting eases deserve special attention. If (3 = 0 one 
of the wave tjqjes is completely coupled, (?; 2 /ai)' = l. F' and 
(TFj' 4- TTj') have the values which are to be expected. The second 
wave type is infinitely slov' (r"=co); the corresponding wave 
impedances are infinite, so that this wave type seems to be of 
little importance. It is, however, the onl}' possible tj'pe in rigid 
materials like acoustic plasters. This case is especially worth while 
studying for more or less flexible wood fibre materials. These are 
generally treated as rigid ones (F''), whereas it is verjr likely that 
the other wave type (IT") is often of much more importance. 

The other extreme case is that of complete decoupling, which 
turns out to occur only if at the same time (u= cc (or <r = 0), 
7i = 1 and h = 1. That decoupling does indeed take place is clear 
from the values of Uj/Ui (0 andco ). Moreover F and IF have the 
normal values. 

If h ^ 1 (last column but one of Table 6) frame and air remain 
coupled systems, even at high frequencies, where the air friction 
is of no importance. 


§ 3 GRAPHICAL BEPKESEXTATIOX OF THE BOOTS OF THE 

f-EQUATIOX 

The roots of the F-equation (3.10) can be computed without 
difficulty. Their dependence on the frequency and the material 



VALUES OF r', r TT/, AITO (vt/tO’* IN SPECIAL CASES (SEE ALfiO TABLE 7) 

ABBREVIATIONS, USFD IN THIS TABLE ONLY, ABF INDICATED UNDFR THE COLUMNS 
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constants is very complex and hardly surveyable. The following 
graphical representation of the dependence on frequency under 
various circumstances may be of great help. 

Writing 

K. =r- 

P__ hp^ + S p, + S p,{l — 1i)+S — — 

hK^ hK^ hK^K^ ■ h 

_p_ P 1 P 2 + (pi + Pi) S 

h 

the r-equation takes the form 

— Eu -h F — 0. 

Its solutions are u' and u" 

2u' ■=E+ VW, 

2u" = E—Vl^ 

if D denotes the discriminant of the P-equation 

B = FJ — iF. 

In order to allow for a graphical interpretation something must 
be assumed about hTi and K„ , and to simplify things we shall 
suppose both Jfj and /fj to be real constants, i. e. it is a.ssumed 
that internal friction in the frame is of negligible importance and 
that the Kirehhoff effect (II § 3) and other effects, w’hich 
might cause to be complex, are absent. As we saw before, the 
loss angle for K is often small indeed and the contribution to 
sound absorption is for far the greater part due to viscosity. 
Pi and P 2 are constants by definition, and AT, are assumed 
to be constants, so the frequency now only enters implicitly 
through S, since 

S'=(fc— 1) p^+W<j/3<n. 

Assuming and to be real and constant means, of course, a 
restriction to the general theorjv 

One sees at a glance that E is of the general fonn a + h/jw, 
furthermore it is easilj' verified that I) can be written 

D = A — B/^= + C/j^ 


(3.15) 
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m which A and B are essentiallj positive constants, C can be 
negative or positive 



hK^ ) 

{K, + X,)* 


(Jf p, 4. p,) 


f'Pt f 

hp, ) 


(ifp. Cp,) (316) 


The character of the freqnencj dependence of D, and therefore 
of u' and it", changes diseontinuously if C changes it sign Prom 
(3 16) one easily deduces that this tak^ place when approvimateh 

^ = and 1 resp (i£ p ,4 p,) 


If the ratio KJK^ has a value in between these limits C is 
negative, in the other cases it is positive 
Fig 20 shows the way m which r' (= ///«) and r" (= 



Fig so 

The fODstraction of r* »nd r ' 


mai be constructed F/2 is represented bj an upward straight 
line parallel to the imaginary asos D is obviously a parabola (see 
(3 15) ) , D/-i has been drawn m Fig 20 for a positne ^alue of C 
For negatne values of C the other half of the parabola (above 
the real asis) would have been found 

The square root of J)/A vields which added to or sub- 
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Tlie consequence of this is to be seen m the values near the 
en d points of the T-curves (the interchanging of V pJRi and 
V kpJhE^ in passing from the first to the second case) 

If Ki/K^ IS varied contijjuonsly the r graph also taries conti 
nuously from the first to the third graph of Table 7 The ‘ hyper 
bola" and 'semicircle” approach each other, at a certain value 
of KJKi the curves come into contact This requires D — 0, i e 
0 = 0 and A — 5/<i)* = 0 from which we see that r' may coincide 
noth T" if KJK^ has a definite value and the frequencj is chosen 
adequately The same coincidence again occurs if D becomes zero 
for the second time 

A high value of a (or low frequency) means that the roots are 
situated at the extreme beginning of the respective curves One of 
the roots is then equal to V {p^-¥ pj)/(Kj -f K^) A low value 
of a (or high frequency) means the opposite, i e , roots at the 
extreme endpoints 

"Which of the two curves of the graph is called the r'-ourve 
or the r”-curve respectively is arbitra ry It seems a good system 
to denote all cun es ending m U'pjK^ by r' and to call the corres 
po nding w ave the disturbed frame wave The curves ending in 
y kpi/hki may then be labeled with r" and considered to belong 
to the disturbed air wave This nomenclature can be justified bj 
considering numerical values of the ratio vjv^ It might be 
expected that the ratio vJv^ for the disturbed airwave has a 
greater value than unity or at least greater than that for the dis 
turbed frame wave This indeed turns out to be true This maj 
be verified with the aid of the expressions for Vx/Uj at b» = 0 and 
o) = oo in Table 6, if is supposed to be very small and very 

large respectively 

Prom a practical point of view it is important to know whether 
frame and air are to be considered as closely coupled or not This 
can be estimated as follcrws For the frame- wave Vi > Vj The 
transition between coupled and rather decoupled vibrations of this 
wave type may be a&sumed to take place where the inertia term 
equals the friction term in magnitude (see equations 

(3 01) and (3 03) for "P Vg), from which the decoupling frequency 
IS found 
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Partial decouplings of the air-irave may be expected for similar 
reasons at the higher frequency 

<!)" = ll^a/llp2 . 

The practical importance of partial decoupling may be seen from 
the following reasoning. At low frequencies (complete coupling) 
one wave type is undamped (r = real), the other is indeed 
strongly damped but absent owing to the high wave impedance. 
As soon as the frame-wave becomes free it wiU begin to contribute 
to the absorption. The air-wave will start its contribution to the 
absorption at a higher frequency. It should be home in mind that 
Ax and A, have been assumed to be real; as this is not strictly 
true, damping takes place owing to the loss angles in /£. and , 
even at very low frequencies. 

To give some idea of the decoupling frequencies in ordinary 
materials a rough estimate is given in Table 8; the steady flow 
value of o- is taken for simplicity. From Table S it is clear that 
the study of coupled vibrations is of the utmost importance for 
many usual materials. 


TABLE S 


BECOHPLLN'G FEEQUEyCIES FOE PEACTICAL lIATEELiLS 


material 

<T 

F = o,V2- 

1 ,/'=:o"/2v 

ac. plasters 


1 — 20 Hz 

100 —2000 Hz 

felt 

5-10‘ — 10'^ 

100 — 3000 

1000 —20000 

wood-fibre 

; 5-10= — 10- 

500 — 15000 

j 10^ — 2-10 ■ 


Figs, 21 and 22 give two numerical examples for sponge rubber 
layers. They may be deemed to be rather flexible. Other materials 
like soft wood-fibre plates, rockwool, etc., may have, however, a 
compliance of the same order. 

Expressions for the decoupling frequencies have been derived 
from the equations of motion, We .should bear in mind, however, 
that also the equations of continuity are coupled and even remain 
so at high frequencies. What a and ti" really mean can perhaps 
best be elucidated with the aid of Fig. 21 and 22, where i-' and 
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v" are indicated on the com^onding curves One sees that ther 
are easy means for giving a rough idea of the frequency scale 
that belongs to each curve The transition points on the r"-cunes, 
1 e , on the curves corresponding to the airwave, maj be thought 
to be too much m the neighbourhood of the region of free 
vabrations to he considered as belonging to the decoupling fre 
quenc:} range It can be shown that at frequencies which are about 
three times smaller than the indicated frequencies the argument 
of r", I e , the damping of the air wave, is about maximum, whereas 
the maximum of the imaginary part of r" ls alreadj reached 



Fig 21 Fig 22 

Th«reti<al enrre* for ft si, 1 = 7, Tbwrctical carves for h = l=s7, 

K, = 0^4 10* ^/m’, 2,= 1 4 10», E,=s04S 10», 2, = 14 10*, = 200 

f , = 120 kg/m’, 1 2, 'T = 2 10* f,ss7S,f=22 10< mks-umts 

C2l»-CC]lS 

at frequencies which are again three times smaller Although 
we did not check this in other cases we might perhaps conclude 
that the decoupling frequencj rattffe ls 0 1 to 1 times the decoupling 
frequencies given above 

( 4 SIMPLIFIED THEOEF, A = 1 

If the porosity is put equal to unitj, as was done in the authors’ 
original publication, the mathematical treatment can be given in 
a much more sjTnmetncal, and therefore, elegant way This is 
due to the coupling between the equations of continuitj cancel 
Img out 
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The most important results of the general theorj' will be 
summarized for this case. 


L K, 


+ 


K, J 




K,K. 


= 0 


w,=K^r 

ir 2 = 

n, — /fiF- + Pi + S S KjT- — p; 

Uj >S' — /fjF” + p, -r jS* K^T~ — p. 


The latter expression for v^/v^ makes it easy to verify that 


(v,M)' (u^ni)"=: — ETi/iT,. 


It should be emphasized that this simple I'elation exists indepen- 
dently of the frequency. For the derivation one needs r' - • F" 
and ^ + T" wliieh quantities are taken from the F-equation. 
If a function x defined according to 


vJv.^X 

(3.17) 

one obtains 


11 

1 

(3.18) 

furthermore from the equations of continuity 


p„ V 2 j hr. 

(3.19) 


which simple relation holds for both wave UTes (insert x' or x" 
resp.). 

Finally one easily verifies with the aid of (3.17), (3,18) and 
(3.19) the reciprocity relations 

v/ v/ ’ Vj" v" 

From (3.18) it follows that one of the x’s is, in absolute value, 
smaller than unity and the other greater, the former belonging 
to the frame-wave, the latter to the air-wave. 


5 5 ELASTIC LAYER WITH CLOSED FRONT SURFACE 

YTe no%v consider the particular case of an elastic layer to 
which the general theory applies, backed by a rigid wall and 

Sound ab^orbicjf matcriali ^ 
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at the front surface closed by a thin non porous covering layer 
of negbgible mass (e g, a tbm lajer of paint) 

The way in which this problem, that of the next section and 
other problems concerning elastic layers may be dealt with is 
in all cases essentially the same I’our waves are supposed to 
be present in the layer, viz, an incident and reflected y'uave, 
and similar waves of the /'-type The amplitude of the first 
wave — the incident /wave — is completely fixed if the velocity 
of the frame at a; *= 0 is given Let us denote this quantity by v'^ 
Once v'li IS given all other quantities v',, , p'li , and p'ji at x == 0 
can be computed since ad ratio’s (Vai/vii)', etc , are fixed for the 
wave under consideration The dependence on x is also knoivn 
through / The same reasoning applies to the other three wave 
types so that the wa\e pattern is completely determined if the 
four velocities at * = 0, /,<, t/,,, v"„ and t/'jr are given Finally 
the problem is to compute the ratio’s between these four amph 
tudes from the boundary conditions and from these ratio’s again 
the impedance at the front surface 
We shall demonstrate this procedure for the case of a layer 
with closed surface, this being perhaps the easiest case 
For brevity we shall introduce single letters for ratio’s of 
pressures and velocities, viz 

(not to be confused with the number 314), 

Vj/Vi c= f (thus distinguishing clearly from x> u’hich ivas used 
in the case 7t = l), 

Pj/vj‘=Wg, 

Pj/Vj t= TFj 

These four ratio’s occur of course with a single or a double 
dash according to the wave type involved 
The boundary conditions are (see Fig 23) 

2Vj = 0 for x — l 
2 Vj = 0 for x = l 
2 Uj = 2 Oj for X = 0 

These conditions give three equations m the four amplitudes 
from which their mutual ratio’s can be computed Written in 
full we obtain 
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f'li exp (— y'Z) v'^r exp /I + exp (— Y'l) + exp Y'l — 0 
f'ji exp (— y'Z) + l/ar exp y'l + v" 2 i exp (— y'l) + V^’^r exp y”l = 0 
u'li + I^ir + t>"ii + v'\r= v'^i + v'zr + t/^i + v''^^ . 

Expressing all velocities in the four frame velocities at x = 0 
yields for the matrix of coefficients 




'V'\i 

r"ir 

exp (— /Z) 

exp y'l 

exp {—Y’l) 

exp y"Z 

e'exp(— yZ) 

Y exp y'l 

Y' exp (— y"Z) 

Y’ exp y"Z 

1 — c' 

i~Y 

1 — Y' 

1-5," 



Fig. 23 

Bonndary coaditions for a layer vith closed front surface 


with the aid of which one easily verifies that 

n'lr/ v\i = — exp (— 2 y'Z) ; v'\r/ v'\i = — exp (— 2 y"l) ; 

v'\i _ 1 — / , coshy'Z _ exp Y'l 
v\i 1 — <p" cosh Y'l ' exp y'l ' 

The impedance at the front surface is 

- / S Pi d* 2 Pa 

\ S Uj /jp _ Q 

V\i + -h V\r + p’.r + P’\i + p".i + p'\r + P'',r 
V'li d- V'jr + V"ii + V"ir 

The double dash velocities of the denominator can be eliminated 
with the aid of the third boundary condition (see the last line 
of the matrix of coefficients) ; the numerator can be written in 
^li 1 t/ir , u"ii , v"jr by first writing as -rpx and then changing 
all pressures into velocities with the aid of p = ±: ITu. 

After introduction of (3.20) one finally obtains 
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Z = 1?', (1 + 7^) ? coft y'l + 

+ W", (1 + ,7") coth y’’l (3 21) 

r — T 

The impedince may be considered as a series connection of two 
impedances We are not allowed to interpret this result thus 
that the impedance is com'posed of that of the frame and that 
of the air In the extreme case of \ery loose coupling (u»= 

^ = 1) and h = l (321) maj be greativ simplified ouing to 
(see Table 6) 

11 'i =\^ AlPl = 03 , 


^ / from (3 04) ond\ 

\(3 05) fo^/l = l^ 

= 0 


11 =: 11 y'Vir", 

11 -,= 


We then get 

z=VKj^w\hjalV'^jK^ + V'KJ,<x«aiJVJJK', (322) 

Now there is no objection against speaking of the frame 
impedance and the air impedance separately 

In the other extreme case of tight coupling the second coth 
in (3 21) cancels out leaving 

z= K(A, + if,) (p, + P,) coth iJ, V (p, + (.,)/(K. + K,) (323) 
as was to be expected 

Let us now discuss the results with a view to the absorption 
coefficient of the layer To begin with we shall consider in more 
detail the case of slight coupling in which the two terms of the 
impedance may be looked upon more or less as giving the 
impedance of frame and air separately Fofh terms are of sprnr^ 
form in the complex plane the spiral of the second term (air) 
being essentially the same as discussed in Chapters I and II 
shelving resonances and anti resonances in the audio frequencj 
range the lowest resonance frequency will occur somewdiere 
between 800 — 2000 Hz 

The first spiral (frame impedance) must be superposed upon 
the air spiral The general appearance of the resultmg impedance 
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contour is materially influenced by the ratio of the lowest 
resonance frequencies of frame and air. If the frame resonance 
frequencies are much lower than the corresponding resonance 
frequencies of the air (p^ » p^-, ca K^; soft materials like sponge 
rubber, rockwool, loose rvood-fibre plates of low air resistance) the 
character is as shorvn in Fig. 24. 

The general upward trend of the air contour is retained. The 
frame resonances cause fluctuations thereon. The most important 



Theoretical impedance contour for a layer ■with closed front surface 


consequence, however, is the shifting of the air contour to the 
right along the distance of the order of the wave impedance 
of the frame (K F,pi). The air impedance for materials of low 
air resistance at low and moderate frequencies normally has a real 
part of the order TTo- Especially at low frequencies an increase 
of this real part materially increases the absorption, which maj' 
be deemed a great advantage. 

In the reverse case of comparatively stiff frame (p, '/> p,; 
If, » » Jf, so that If ,/pi A'j/p- ; acoustic plasters, stiff wood- 



70 ABSORPTION BY ELASTIC LAYERS III 

fibre plates) the character <rf Fig 24 changes in so far that now 
the general upward trend is doe to the frame spiral, while the 
fluctuations are to be imputed to the air, these fluctuations now 
occur at much higher frequ».ncies (e g, 2500 and 5000 Hz instead 
of 440 and 8S0 Hz) while the shifting to the right O '' is 
so great that the contour conies in. the low absorption part of 
the complex plane 

The latter case is generally known, plasters and stiff wood 
fibre plates no longer absorb when superfienllr closed The 



Fig 

Improrwnrat of the sbsoxpticoi bj cJcwng the 'Urface 
of a very clastic layer 1 open 2 clo^ 

former case (superfieiallv closed rather compliant materials) is 
the most interesting one Much better results than with the same 
lajers in non closed condition maj be obtained especially m the 
important low frequenct range In Fig 25 an example of 
measurements is given for a sponge rubber layer of ren Jou 
air resistance The impedance curves (Fig 26, smoothed) indeed 
show the expected form 

The intermediate case, where the resonance frequencies of frame 
and air occur m the same frequency range gives rise to ^erv 
complex results 

More results and detailed discussions of the agreement between 
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experiment and theory may be fonnd in the original publication 
If a material has a small coupling between frame and air we 
are in a position to increase the thickness of the layer of air 
without changing the thickness of the frame by placing the layer 
at a certain distance before a rigid walk 



Impedance contours to fig. 2a 

Fig. 27 shows results of such measurements, carried out on 
the same layer to which Pig. 25 applies. Indeed we succeeded 
in lowering the resonance frequencies of the air, thus obtaining 
a compact impedance curre in the high absorption part of the 
complex plane. 


» C. W. Kosten and C. Zwikker, Plii/sica, 8 (1011) 933. 
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Up to now the coating layer was assumed to be massless, a 
hypothetical case as a matter of fact This mass may ^e^y easily 
be taken into account, if, namely, the mass per unit surface 
be denoted by m we only have to add to the impedance This 
means an increase of the upward trend of the impedance curve 
At high frequencies jiow maj be appreciable The coating layer 
thus gives nse to a rather vague cut off frequency For the low 
frequency absorption jim, is advantageous Another rather im 
portant advantage of resilient, superficially closed, layers is doubt 
less that they can be washed and painted 



Fig 27 

Impedan'ce contour < f a superficially closed layer backed 
with a layer of air of 78 mm frequencies indicated in 100 Hs 

$ 6 ELASTIC LAIEB WITH OPE^ FRONT SURFACE 


This case is perhaps of more importance in practice 
mathematics are a little more complicated 
As boundary conditions, we now have 


SOi = 0 

2 Vj = 0 

SPi _ 2p*^ 

1—h h 


for xe=l 
for X‘=l 

for x~0 


The 
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The latter condition is obvious since both frame and air at a: = 0 
are acted upon by the outer air pressure, the surfaces per unit 
surface of porous material being 1 — h and h respectively, the 
forces Pi and p„ by definition. 

If again four vaves are supposed to be present, given by 
v'lr, V\r (velocities at x = 0, i. e., at the front surface), 
the third boundaiy condition may be expressed in these variables 
as was done in the preceding § 5. 

AH pressures p, are expressed in pj by — then all re- 
maining pressures p^ are changed into the corresponding velocities 
with the aid of p = ±'Wv. We obtain 


— (i/i. — V\r) w 

+ W" { h ~ (1 —h) — v’\r) = 0. 

From this equation and the first two boundaiy conditions, which 
are equal to those of the preceding section, all ratio’s between 
^-''n f v'lr , -F'l. , and v’\r can be computed. 

After insertion in the expression for the impedance 


ut 

' VS(1 — 


S Pi + 2 Pj 


2 (1 — h) Vi + 2 liv^ 


L. 


we finally find after tiresome but elementarj* calculations 
1 (l~7!.-bV) {h— (1— 


Tr'i(7rW') 

{l—h + h^") { h—{l—h)7r' } 


tanh -/f-r 


tanh Y'l. 


(3.24) 


17" (^W) 

For materials for whicli Ji — 1 holds to a good approximation 
equation (3.24) simply runs 


ir'i(rr'-^') 


tanh y'l + 


17"i(Tr"— tt') 


tanhy'T (3.25) 


The impedance turns out to be the result of a shunt-connection 
of two impedances. In the case of slight coupling the impedance 
may be considered to be the air-impedance, more or less short- 
circuited by the frame-impedance. This short-circuiting is most 
effective for small frame-impedances, i. e ., when the frame is in 
the neighbourhood of resonance. 

Theoretical results of shunt-connections like (3.24) may best 
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be represented in an aeonstical admittance diagram From the 
general equation, giving the absorption coefficient as a function 
of the impedance (equation (112)) one sees at once that the 
absolution coefi}C 2 ent does not ebange, when instead o£ the 
impedance the admittance is msertcfd, provided impedances are 
measured in IVo as unit This means that the hues of constant 
absorption in the complex admittance plane are the same circles 
as used up to now m the impedance plane 



In Fig 28 an example is given The assumptions underljmg 
this figure are slight coupling, Pi ^ p 2 . Ki » The second 
term {air admittance) is a normal spiral with its first resonance 
at about 1000 Hz A complex tanh is represented by a spiral like 
a coth, only it begins in the origin instead of at infinity It 
should be noticed that the apex of the spiral shifts towards the 
ongin if the \\a\e-impedaiiee increases In our case, therefore 
the frame spiral is small and turns rapidly owing to the loi\ 
resonance frequencies The superposition of the two spirals gi' es 
a more or less modified air spiral (dotted in Fig 28) The m 
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fluence may be of some importance in the low frequency range. 

If the coupling is not small the terms air-spiral and frame- 
spiral gradually loose their significance. In the limit of very 
tight coupling one of the spirals cancels out, the remaining being 
related to the normal coupled wave (see equation (3.23)). 


i 7 EiEASTIC LAYER BACKED OTHERWISE THAN BY A RIGID WALL 


If an elastic layer is backed by a membrane of impedance Zb , 
the boundary conditions (S Uj — 0 and S Uj = 0 for a rigid wall) 
should be replaced by 



S Pi ~h S Pa \ 

— /i)Ui + 2 hv^ 'x=i 


Uj = Uj at X — 1. 


The condition Uj — expresses the fact that the receiving material 
of impedance Zb is essentially not double-refracting. If this material 
is a rigid porous one, we have to substitute the boundarj^ con- 
ditions by 


Oi = 0 at x — l, 



at x — l, 


if now Zb is the impedance of the receiving material as for instance 
measured in an interferometer. 



CHAPTER IV 


EXPERISIENTAL DETERJIINATIOX OP THE ELEMENTARV 
ilATERIAL CONSTANTS OF POROUS ABSORBENTS 
GOVERNING THE ABSORPTION COEFFICIENT 

} 1 INTRODUCTION 

A few elementary material constants entered into the theoretical 
considerations of Chapters I to IH, viz , the porosity h, the structure 
factor k, the air resistance <r and, for elastic lajers, the “compres 
Sion modulus” of the frame As we saw (I f 8) the structure 
factor can hardly be expected to follow from anj statical or steadj 
flow measurement It is rerj compbcately linked up with the 
structure of the material. The porosity can easily be determined 
from steadj -state or low frequency measurements, although there 
IS some evidence (see II § 9) that the value thus found is too 
great to meet the acoustical situation Anjhow, the measurement 
of the steady state or low frequency value will give us a rough 
estimate of the true value The measuring method for determining 
this value is given m IV § 2 

The dynamical value of the air resistance can be computed from 
the steady state value, provided we know whether Poiseuilles 
law or Helmholtz’s theory applies Therefore the steadj state 
value IS no doubt of importance (see IV § 3) 

The measurement of A, seems easj, but reallj it is not "What 
we really nant is a knowledge of the complex quantitj Ki, which 
can only be obtained dynamically It se«ns plausible that a statical 
measurement of the stiffness of a sample should give us the modulus 
of /f, However it is a rather well known fact that the djnamical 
modulus may be several times the statical one At low frequencies 
(say 5 — 50 Hz) Ki normally vanes so bttle with the frequencj 
that an extrapolation to acoustical frequencies seems allowed One 
disagreeable restriction should however be made We ought to 
measure the frame stiffness freed from any air stiffness So if a 
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is high (tight coupling) the measurement should be carried out 
in vacuo; if o- is low it may be sufficient to provide for ample 
passage to the surrounding air. Tliis difficulty may be avoided 
by computing from the measurement of the rigidity modulus G 
(against shearing). Properly speaking this only holds good for 
isotropic materials, but also for non-isotropic materials (e. g., wood- 
fibre plates) this procedure will probably provide us mth a reason- 
able accurate value of the complex quantity K■^ (see IV § 4). An 
idea of the degree of inhomogeneitj- might moreover be obtained 
by taking samples of different orientation. 

§ 2 THE POEOSITY 

Making use of Boyle’s law pV = constant, the porosity can 
be determined with the aid of the apparatus shown in Pig. 29^ 
The specimen of known external volume is placed in the ve-ssel 
of volume Fq this being filled up, when neeessarj’-, unth glass 
plates to reduce the dead volume. Inei’ease of volume is effected 
by opening the tap F, and by adding the volume of one or more 
of the spheres V'', V", etc., to F^ by lowering the vessel B. The 
change of pressure is measured Yuth the manometer E. Mean- 
while the level at C is kept at the same height with the aid of 
the taps at the top of tube D. 

The following details are worth mentioning. Both the set of 
bulbs V', V", etc., and the manometer part E — D have to be 
placed in a water thermostat. As manometer fluid, o-di-iso-amyl- 
phthalate was chosen, because of its combination of low density, 
low volatibility and low viscosity. Some results may be 
mentioned here. 

Porous plasters have porosities (/t) between the values 0.50 
and 0.75. 

For wood-fibre plates, values were found in the neighbourhood 
of 7^=r0.80. 

Pelt, sponge rubber, and cork samples all have high porosities, 
varjung between 0.85 and 0.95. 

Recently an attractive device was described by Leonard- 


1 .1. Van den Erjk and C. Zvikker, Fhysica, 8 (1941) 149. 

2 B. tv. Leonard, J. Acmtst. Soc. Avi-., 20 (194S) 39. 
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which permits the measurement of the porosity -mider low 
fiequency conditions At extremely low frequencies the wave 
character of the sound inside the absorbing material may be left 



Fig 20 
Porosity meter 


out of consideration The air content mainly behaves as a spring 
the compliance per m* of which is proportional to the total amount 
of air per xn* included in the material, le, proportional to the 
thichncss and the porosity The main features of Leonard’s 
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apparatus may be seen from the sehematieal figure 30. A resonant 
system is built up of rrhich the stiffness is for the greater part 
clue to the air content. Therefore the porosity can be computed 
from the resonance frequency. Leonard claims as advantages 



Fig. 30 

Leonard’s porositv meter 


that his method is not sensitive to changes of temperature, a 
great advantage as a matter of fact; moreover the volume of 
difficultly accessible pores and holes does not contribute to the 
porosity found. Indeed (see Chapter II § 9) the effective porosity 
tviU decrease cvith increasing frequency. 

« 3 THE AIR EESISTAXCE 

An apparatus for the determination of the stationary air- 
resistance a is given in outline in Pig. 31. On the specimen A 



Fig. 31 

Air resistance meter 
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a bo\ B IS placed m \>lueh a small subatmosphene pressure is 
maintained 11 ith tbe aid of an air pnmp This pressure is measured 
ivith a F u e s s micromanometer D The % elocity of the air current 
IS measured by a capilJarj’- flowmeter C In order to present 
any radial currents in the specimen, bo'^ B is enclosed tn a uider 
lox F in which the pressure is kept at etactlj the same values as 
in B (principle of the guard nng) This etjuabti is cheeked bi 
means of a drop of o-di iso amyl-phthalate H, ^ihich is to remain 
at rest inside the capillary J Bo-^ P is fastened along the edge 
on the specimen A bj means of plaster of Pans 

Practical \ allies of the air resistance vary between 10* and 
10^ mfe units In order to meet this wde range different capil 
larj tubes ha^e to be used, or what comes to the same thing 
hut IS much easier the nir resistance of the capillary tube may 
be changed by inserting one or more chroinel wires of different 
diameter m the tube 

An apparatus was constructed along these lines m the laboratory 
of the authors The calibration characteristics of the flowmeter 
were not straight Imes, neither when the \olurae \elocity was 
plotted against the pressure difference, nor when it was plotted 
against the difference between the squares of the pressure 
Pressure differences of the order of 002 to 2 mm water height 
(le, 02 to 20 N/m*) were used, yielding air-currents of about 
3 to 180 emVsee The inner surface of the samples was always 
400 cm’ Specific air resistance between 10* and 3 10^, therefore, 
could be measured The relatne error amounted to about 15% 
(see for practical values of €f II § 5) 

Measurements without a guard rmg are also common practice 
The cylindrical samples, needed for the measurements in the 
acoustic interferometer, are then used For the air resistance 
measurement the side surfaces of the sample ha'v e to be closed 
with solid paraffin in order to procure a parallel flow of air 
through the sample 

Leonard’ developed an apparatus for measuring the air 
resistance in which an analytical balance is used at the same 
time as air pump, pressure meter and velocity meter (see Fig 32), 


R W Leonard, J Aeoust Soe Am, 17 (1946) 240 
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which can be used for values of air-resisfanee as low as 1000 mks, 
using air-pressure differences as low as 1 X/inv The accuracy 
is about 2 %. 

The uniformity of a material may be checked by taking 
samples of different thickness of one and the same material. 
As a rule the surface of every layer is less impervdous than the 
interior part. 

For high absorption erI/3 should be of the order of (420 mks) 
(see H § 7). This demands specific air-resistances of the order 
of 5 • Ifr* mks ; these are found in high-resistance plasters and low 
resistance felts. 



Leonard '5 air resistance meter 


§ 4 THE COltPEESSIOX 3IODULUS 

As pointed out in IT § 1 hC, cannot be measured statically. 
If measured at low frequencies, the results have to be extra- 
polated to audiofrequencies; moreover due account should be- 
taken of the influence of the enclosed air. The air may virtually 
increase both the stiffness and the losses. 

In order to avoid measurements in vacuo we can have recourse 
to shearing measurements. Since shearing of a material does not 
give rise to a change of volume the stiffness of a sample against 
shearing is independent of the air-content, and since displacements 
of air in the interior might be a.ssumed to be absent or .small, 
the air will not materially increase the losses. So the rigidity 

Sonnd ah5orbtn?t inatrriala ^ 
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modulus G ■will practically coincide ivitli the value in vacuo 
Ej IS, according to its definition, the stiffness of umt volume 
of porous material if sideway extensions or contractions are pro- 
hibited Xow this IS neither the definition of the compression 
modulus K nor that of the modulus of elasticiU (Young’s 
modulus) E If the material under consideration is isotropic, one 
can easily derive 

2 G (»i —1) (tn —2), 
where in is Poisson’s ratio (»n>2) 

Now for all practical porous materials in will be rather high 
(sav 10 or more) so that we maa approximatelv put in equal to 
infinity, yielding 

A, ss 2(? 

The measurements can be accomplished nith anj one of 
the apparatus for measurmg mechanical impedances Suitable 
apparatus h’i\ e been described by if e y e r and K e i d e 1 *, 
Costndoni*, Bohme* and the authors* In principle the 
normal acoustic mterferometer maj be used for this purpose too, 
though special precautions h*i\e to be taken 
If the measurements are carried out at audiofrequencies anth 
samples, the resonance frequencies of which be in the frequeuoj 
range used, the impedance uiU be spiral like as in Chapter I 
Marrelloiis results of this type have been obtained by Bohme* 
with solid materials (eg, rubber) The authors measured at low 
frequencies so that inertia effects are of no importance Then, 
of course, the results most be extrapolated very far 

5 5 KEbULTb or ilEASlTREMEATS OF THE COMPLEX STIFFNEisS 

Experimental data on the specific stiffness and the loss 
ancle S of porous matenat from measurements in vacuo or 
derived from shearing experiments are not available 
The static ^'alue of A, is readily measured Its magnitude maj 


I E Mever Mitl L. Eeidel, Z terA IS (1^37) S'>£* 

» B C Costadonj, Z IccA PhytU., 17 ( 1026 ) lOs 
» JL Botiae, Atvat Z, 2 (1937) 503 
* C tv Kostea aad C Zirikker, Phynfa, i ( 1937 ) 221 
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be as low as Vioo of tfiat of air for soft fibrous materials, 
for soft sponge rabber and up to practically infinite values for 
plasters. 

Tbe best estimate of the dj'namic complex value of is 
perhaps obtained by multiplying the static value by the ratio 
•Ififlin./ffistat. obtained for non-porous materials of the same kind 
and assigning a loss angle to the value thus obtained at least 
equal to that of the solid material. It is therefore worth while 
to summarize briefly a few results of the ratio siat. and 

of the loss angle S for solid materials. 

An instructive example of a soft non-porous material for which 
data are available ^ is rubber. If properly vulcanized and relatively 
“pure”, i.e., without unnecessary filling materials, dynamic and 
static stiffness are practically equal, the loss angle is only of 
the order of 1®. If the material is more or less luider-vulcanized 
the dynamic stiffness exceeds the static one by a factor of 2 
or even more. At the same time the loss angle increases to 
values of 5° — 10°. Moreover the material is also less ideally 
elastic under static load; after removing the static load the 
initial height turns out to be less (permanent set, flow). The 
same effects as caused by under-vulcanization are foxmd by using 
filling materials like barium sulphate, carbon black, etc. Sjmthetic 
rubbers generally show the same behaviour as imder-vulcanized 
natural rubber, i.e., greater dynamic than static stiffness, loss angles 
between 5° and 10°, although it greatly depends, of course, on 
the kind and manufacture. The behaviour of this class of rubber 
and rubber-like materials is governed, therefore, to a great extent 
by elastic after-effect during the deformation. 

It is not within the scope of this book to enter into the 
theoretical aspects of the djmamic behaviour and the internal 
friction in metals, wood, etc. Since the characteristic frequencies 
of bodies of such solid materials, as computed from the static 
elastic constants agree with those actually found, it may be 
inferred that static and dynamic stiffness for these matei'ials 
are practically equal. The loss angle is mostly small, say 0° — 1°. 

It is very doubtful whether porous structures of similar 


I C. \V. K; OS ten, Proc. Jlubh. Tcclin. Conf., p. 9S7, London 1938. 
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materials as wood fibre plates hair felt glass silk etc ha%e the 
same behaviour Sleasurements earned out in air ' yield dynamic 
stiffness values up to 10 times the static one and in extreme eases 
even higher Friction between the fibres may increase the losses 
the loss angle S vanes between 2® and 20° Values of B greater 
than 10° are almost undoubtedly due to the streaming of the air 
in and out of the sample In the case that the flow of air eon 
tributes much to the losses both K-i and 8 are strongly dependent 
on frequency ® The low frequency values of and 8 approximate 
to the values in \aeuo It should be emphasized that whereas for 
normal incidence of sound on an infinite layer of absorbing mate 
rial the coupling is highest at the louer frequencies the reverse 
IS true in measuring the compression of small samples where the 
air may escape sidewajs The sample behaves as an air-container 
with leakage requiring a minimum time to assume the atmospheric 
pressure inside after having been pumped up b> compression 


C Van G jn and C Zwtkkar Ingen e r (Dutch) 56 (1941) 0 11 
E Meyer and L KeideJ Z Ud Ptyn! 18 (1987) 299 E Mejer 
Z Ver deiif Ing 78 (1984) 19o7 

* 0 Koaten and C Zwikker 


Fhysica 4 (1937) 843 



CHAPTER V 


3IEASUREiIEXT OF XORiLAL lilPEDAXCE AXD 
ABSORPTIOX 

} I I^'TBOD^CTIO^" 

Practically all devices for measuring the impedance of 
materials (from which the absorption coefficient may be derived 
with the aid of Fig. 11) make use of a measuring tube in order 
to create plane travelling waves. At the one end the sound 
source is placed, whereas the tube is terminated at the other 
end by the sample tinder test. The number of methods for 
measTiring the impedance is practically unlimited. The sound 
pres.sure is measured as a function of the distance from the 
sample, of the frequency or of the length of the tube. Three 
adequate quantities are neeessarj' in general (e.g., two pressures 
and one distance), provided damping of the travelling waves in 
the tube owing to the Helmholtz — Kirchhoff effect is 
neglected. A rather complete survey of the various possibilities 
has been given by B e r a n e k It is impossible and unnecessary 
to mention all methods here. ITe shall confine ourselves to the 
description of a few' of them, viz. 

a. constant length method, measxu-ing maximum and minimum 
pressure in the tube (see § 2). 

b. variable length method, measuring maximum and minimum 
pressure at the sound source (see § 3). 

e. variable length method. mea.snring curve width of pressure 
at the source (see § 4-), 

d. variable length method, measuring the electrical impedance 
of the sound source (see §§ 5 and 6), 

e. short length tube, measurement of pressure and velocity 
(see § 7). 


* L. L. Beran'clc, J. Jccmst, Soe. Aki.. 12 (1010) 3, taWe I, p. 1. 
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5 2 CON'STKNT LENGTH INTERFEROMETER 


If the sound source operates at sudi a !o^r frequenc.^ that the 
\\aAe length is greater than 1707 times the diameter D of the 
tube, onlj plane imcs can tra\el in the tube'' In the neigh 
bourhood of tlie source these plane waves arc disturbed if the 
sound source generates non plane waves The same holds for the 
neighbourhood of tlie sample A non homogeneous sample mii 
reflect the incoming plane nave as a distorted reflected one At 
some distance from tlie sample the reflected nave mil be plane 
again if the wave length is long enough Therefore nil interfero- 
meters can only be used np to frequencies loss than cJllD This 
restriction can not be neglected if the source provides us nith 
ideallj^ plane waves since all samples arc more or less inhomo- 
geneous and will reflect travelling non plane naves if the fiequcncj 
IS chosen too high Assuming tlint this restriction is fulfilled and 
denoting tlie impedance of the sample bj e we maj desenbe tlie 
sound field m the tube as the superpcisition of tuo plane waves 
an incoming and a reflected one On mg to absorption the reflected 
one will bo of smaller intensity than the incoming one We mav 
borrow from equation (1 12) tlio relation between the reflection 
coefficient r=p,/pi directly before the sample and the impedance z 


z 

IrT 


1 + 1 
1 — r 


(5 01) 


Remembering that the propagation constant y of nr is jm/Ca ne 
can nTite down the pressure at a distance x m front of the sample 

p(x) =pj exp ( — Jvix/Cf,) + rp, exp (+ jux/Co) (5 02) 

The exponentials have modulus umtj and onij clnngc tlio phase 
of the travelling waves (damping is neglected I) 

Writing as previously 

r=|r|evpja, 

(5 02) takes tJie form 

p(x) =p{exp ( — jux/Co) + (r (pj evp j(ci>x/Co + A) (5 03) 


RayloigU, Thtory of Sound 11, 2nd ed, p ICl, Csmbiidge 1S«<5 
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§ 2 


The pressure ‘p{x) vrill be a minimum for those values x„ of x 
for which both terms are in opposite phase, i. e., for 


{axjcg + A) — ( — laxjcg) =z — (2 n. + 1) 
where n is an integer. The latter equation yields 

A = — (2u + 1)^ — 4 WX. (5.04) 

This equation offers us a means of measuring A. A may be com- 
puted. from the position of the places of minimum sound pressure 
in the tube. It turns out to be too inaccurate to compute A. 
from frequency and sound velocity. More reliable values of A are 
obtained bj* measuring the positions of Uvo minima (e. g., Xo and 
x-y) and eliminating A. One sees at a glance from (5.03) that the 
ratio of maximum and minimum pressure is 


Pmai. ^ 1 ^ i 

Pmln. 1 1 ^ ] 


(5.05) 


This ratio can easily be measured, from which | r [ is found. So 
(5.04) and (5.05) are a basis for measuring modulus and argu- 
ment of r. 

Numerous instruments have been designed along these lines, 
mutually differing only in the way the pressure is measured. 
Turn a ^ used the ear, and attained only inaccurate results. 
Taylor® used the Rayleigh disc, placed in a side tube. Paris® 
made use of anemometer wires. So did Goldbaum and 
AVaetzmann^, Eckhardt and Chrisler® used an 
exploring tube, a tube of very small diameter in comparison with 
that of the interferometer, forming an acoustical connection 
between the microphone outside the interferometer and the 
measuring point inside. This technique was, then, adopted by 
numerous other investigators. Although the use of an exploring 
tube is convenient and rather reliable it cannot be denied that 
it is a source of error. Its impedance is not infinitely high, so 


1 J. Turn a, Sits.ber. Kaisert. Akad, JTiis., 3 (1902) 402. 

2 H. O. T a y 1 o r, P/iys. Bev., 2 (1913) 270. 

“ E. T. P a r i s, Froc. Phys. Soe,, 39 (1927) 269. 

•* G. Goldbaum and B, Wact zman n, Z. PhysiJ:, 54 (1929) 179. 

5 E. A. Eckhardt and V. L. Chrisler, Bur. Standards Sci. Paper 
.Vo. oSe (1926). 
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It disturbs the sound field Moreover even it its impedance uere 
infinite it -would disturb the field by its geometrical shape 
Rough calculations show that errois of the order of the ratio 
of the cross sections of the exploring tube and the interferometer 
tube m the impedance are to be expected A method that has 
the great adiantage of not disturbing the sound field is shown 
in Pigs 33 and 34 The microphone is built into the tube wall 



Fg 33 

Constant length interferometer wilhont explonng tube 



(see also Fig 35) behind a stiff membrane and can be moted 
together with the whole tube m -which loudspeaker and sample 
holder fit m an airtight way like pistons Th s method seems 
to have been developed independently by -various im estigators 
It might seem a source of error to use a membrane of 1 cm* 


1 tv AI Hall J Acouit Soc Am 11 (l^S) t40 J Harmans 
Alust Z 5 (1940) 215 J t d Egk C tV Kosten and tV Kok Jppl 
set research B 1 (1947) 50 
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or more, since it will be impossible to measure in that way the 
sound pressure at a “single” point. This objection, however, is 
not correct. If, for example, r = l p„,n. is zero; the pressure of 
the standing wave varies as sin ax/Cg Avith the place in the tube, 
if a; is the distance from the point under consideration to the place 
of minimum pressure, so the pressure is not zero in the Aucinity of 
this point, but is in opposite phase at either side; no force is 
therefore exerted on a membrane, placed symmetrically AAoth respect 
to the measuring point. In a similar way it can be shoAvn that also 
in the general case | r | < 1 no error is made in measuring the ratio 
Pmaj /Pmin. (sec Van den E y k, K o s t e n, and K o k page 88 ^) . 



Details of the raicrophoire of the interferometer of Fig. 33 
f 3 WEXTE AED BEDELL’S METHOD’ 

If the length of the tube is changed by moAung the sample in 
the tube and provision is made to keep the strength (velocity) of 
the source constant, the sound pressure at the source Avill be pro- 
portional to the impedance of the tube at the source, looking into 
the tube terminated bj^ the sample. The ratio of Pmu. to Pmm is, 
therefore, equal to the ratio of Zmai and Zmm • If the impedance z 
of the sample is knoAATi, aa'c might compute this ratio with the 
aid of equation (1.04), a tiresome way, as a matter of fact. "With- 


’ E. C. AVente and E. H. Bedell, TScll Tech. J., 7 (192S) 1. 
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out any calculations the result may be read from the circle dia- 
gram, Fig 11 If the length of tiie tube vanes, the impedance at 
the source must run through a circle of constant absorption in 
Fig 11, for the absorption of the sample alone is equal to that 
of the sample with a layer of air in front of it The extreme 
values of the impedance obviously are the points of mtcrsection 
of the impedance circle with the real axis, for which one easily 
finds from (5 01) 


^(MtX 




1 T |rl- 


Therefore, the ratio of p®*, and pm,a is 


Pmla VT— 


(5 06) 


the same expression as (505) except for the second power (5 06) 
IS the basis for measuring |r| The argument A is readily obtained 
from the tube lengths at which extreme values of s are found, 
again witli the aid of Fig 11 AVjth increasing length of the tube 
the impedance runs through the absorption circle in a clockwise 
sense, covering 2 » in A for evorj half wave length So the first 
maximum will be found for 


the first minimum for 

(a + ir)/2ir, etc 

In order to provide a sound source with constant strength 
Wente and Bedell used a telephone acting upon a heaiy 
piston, the sound pressure at the source was measured with an 
exploring tube 


f 4 ODTECVK WTOTir 

This method resembles "NVente and Bedell’s method in so 
far that it is also a variable length method with a source of 
constant strength, what is measured is also the pressure at the 
source end of the tube There are, however essential differences 


t L L Be r 


k J Acovst Soe Avt, 12 (2W0) 3 
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To begin -vvith, B e r a n e k uses a “point” source at the centre of 
the source end by supplying the sound energy through a high 
impedance thin tube. The sound source, therefore, has a (high inter- 
nal impedance and gives a velocity' independent of the load. The 
pressure at the source end is a function of the radius, since the 
waves are not plane in the \'ieinity' of the source. B e r a n e k 
derives this function and shows that the most adequate place for 
measuring the pressure is at the source side quite near to the 
wall of the interferometer tube. Furthermore he does not measure 
maximum and minimum pressure, but derives the damping of the 
standing longitudinal wave from the curve width near maximum 
pressure if the tube length varies. This has the great advantage 
that no measurements are necessary at low pressure level. Damping, 
caused by energy losses at the interferometer wall itself is taken 
into account by determining the damping with a perfectly' reflect- 
ing wall instead of the sample. Prom these two values of the 
damping and the two corresponding resonant lengths the impedance 
can be computed. Various provisions have been made for improving 
the accuracy, e.g., a precision screw for accurate measurement of 
the displacement of the sample. 

Beranek’s method is original and in many respects admir- 
able. His proof of orthogonality of the eigenfunctions seems to 
be incorrect (loc. cit., p. 7), since the ratio of the normal 
component of the particle velocity to the pressure should be 
indei)endent of the eigenvalues at the tube wall, which certainly' 
is not the case. Jloreover it seems questionable whether all 
refinements of his apparatus are necessary in view' of other 
sources of error, e. g., mounting conditions of the sample, dif- 
ferences between various samples of the same material, the 
uncertainty of the definition of what is called the surface of 
the sample. 

§ 5 REACTIO^'-0^'-SOTJKCE lyTEEFEBOlIETEES 

The electrical impedance of all electrical sources of sotmd 
depends more or less upon the acoustical load. If the source 
is prevented from vibrating by fixing the cone, the membrane 
or the emitting surface in general, the electrical impedance, as 
measured in the normal way' at the input terminals of the source. 
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Will be called the blocked unpedanee If the emitting surface 
can vibrate, a different impedance Z is found The difference 
Z — Zi IS generally called the motional impedance, it will be 
labeled Zmot Now depends upon the acoustical impedance 
at the vibrating surface It is, therefore, plausible that the latter 
can be calculated from the measured value Zmot The method is 
extremely attractive in that the whole measuring equipment is an 
audiofrequency generator and an electrical impedance bridge 
Attenuators, amplifiers, microphones, etc , are no longer necessarj 
The method has the disadiantage that it has not yet been worked 
out in all details The results, obtained up to now, justify the 
opinion, however, that it can turn out to be of great value 
It can be shown that for all normal sources the following general 
equation holds 




+ z 


(5 071 


in which i, ilf and A are m general complex constants which 
should be measured as a function of frequency and z is the acous 
tical impedance at the source We shall derive this equation for 
the electrodynamic loudspeaker and the telephone 


a Reaction on an cUctrodynamte loudspeaker Let 
u = voltage on the coil, 
t = current through the coil, 

B = the constant magnetic induction in the slit 

I = the total length of the winding, 

II = velocity of the coil, 

then 

u=-Zit + Blv, (5 OS) 

in which Blv IS the e m f of induction due to the motion The 
equation of motion, in the case that the onlj external force acting 
on the system is the reaction force of the acoustical load, maj he 
written 

0 = — Bh + Z„eci. V, (5 09) 

where Bli stands for the Lorenta force on the winding Zmecti 
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is built up of two parts, the internal mechanical impedance of 
the cone (mass, stiffness, and resistance) and the mechanical im- 
pedance at the site of the coil due to the acoustical load. 

Elimination of v from (5.08) and (5.09) and putting Z for 
u/i gives 


Z^Zi+ . (5.10) 

If the cone were replaced bj^ a piston and if this piston were 
acoustically loaded vdth a specific acoustic impedance z, constant 
over the entire surface 8 of the piston, e. g., by placing the source 
in a suitable interferometer tube, we should have 






piston 


+ zS, 


which brings (5.10) into the general form (5.07). A normal loud- 
speaker does not in the least comply with this idealized situation. 
It was shown, however, by one of the authors ^ that (5.10) may 
be transformed into (5.07) in the general case that the velocity 
of the vibrating surface varies over the surface and the loud- 
speaker does not fit into the tube wthout leakage. Now M and N 
of (5.07) become complex functions of frequency which have to 
be determined, moreover L is not the blocked impedance, but the 
impedance that is found when the acoustical load is infinite. Owing 
to flexibilitj'^ of the cone and air leakage around the speaker the 
coil is not completely blocked if z — ca. 

h. Reaction on a telephone. K e n n e 1 y and Pierce' were 
probably the first to measure the mechanical load of a telephone 
by its influence on the electrical impedance. Pay and Hall® 
applied this method to a tube of %mriable length, closed by a 
sample. The validity of (5.07) in this case is simply proved by 
Pay and Hall from eneigy considerations. AVe shall follow 
their derivation, •which they also applied to the preceding case. 
AVe must, however, define verj' precisely the quantities used. If 


1 C. 'W. K 0 s t e n, Appl. sci. Hes. B, 1 (lO-iT) 35. 

2 A. E. Iv c n n c ly and G. tV. Pierce, J. Fran7:Un Inst., 200 (1912) 46". 
2 E. D. Fay and IV. JI. Hall, J. Acovst. Soc. Ant., 5 (1933) 46. 
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no current is supplied to the telephone the permanent magnet sup 
plies the equilibrium flux that exerts an attractue force on the 
iron membrane If we give a small displacement to the membrane 
the flux, and thereby the force, will increase This force has 
nothing to do -with any electrical input It is considered as an 
internal force of the magneto-mechamcal system The magnetic 
setup, therefore, behaves and is treated as a sjstem with a 
negative stiffness If at the same time a current is supplied, an 
other contribution to the flux is made, increasing the force on 
the membrane This force, which originates in the current will 
be called the force on the mechanical system, and is labeled F 
A similar distinction must be made m the e m f of induction due 
to variations in flux If the membrane is blocked and a current 
is supplied there is an e m f of self induction If at the same 
time the membrane is guen a velocity there oecuis an extra emf 
of induction caused by the velocity This emf wll be called the 
back em.f and be labeled e 

If furthermore i be the current and v the velocity of the point 
to which T IS applied the extra power supplied at any instant 
to maintain the motion ct, will be equal to the product Fv, i e , 
to the power transmitted at that moment by the system therefore 

et =: Fv 

Now if e =: Kv 

we ha%e t = F/t, 

and therefore c/* = A* v/F 

or (5 II) 

From this derivation it is clear that ^mech is the sum of the 
internal mehanical impedance of the sjstem negative magnetic 
spring included and the mechanical impedance of the external 
acoustical load, as felt at the driving point of F We have now 
arrived at an equation corresponding to (5 10) of the electro 
dynamical loudspeaker and follow the same method to transform 
the result into the general equation (5 07) in order to allow for 
non plane motion of the membrane and air leakage in the inter 
ferometer (see K os ten, page 93*) 
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5 6 THE DETERMINATION OF THE ABSORPTION COEFFICIENT 

MHTH REACTION INTERFEROMETERS. PRACTICAL ASPECTS 

Equation (5.07) is the general basis for measuring impedances. 
The constants L, M and N, however, have to be determined by 
calibration. This can be done by measuring the electrical impe- 
dances , Zo and Z, if the source is loaded bj- three knovm 
acoustical impedances z, , Zj and z^ . 

We can then express L, M and N in Z.^ , Z^ and Z, , by which 
(5.07) can be shoAvn to assume the form 

z _ . Z, — Z 3 Z — Z, 

P(fia ^ — ■2’2 ^ — Z, 

The quotient {Z-^ — Z^j{ZT^ — Z,) has to be looked upon as a 
new constant of the apparatus, being a function of the frequency 
only. 

If we choose for Zj , z^ and z^ arbitrarily the impedance of a 
column of air of “/g , V 4 and ’/g wave length respectively, we have 

2 i = ipoCoi Z 2 — O and Zg = x., 

Since only differences of impedances enter in this equation, which 
are, as a rule, one order smaller than the impedances themselves, 
some simple sort of compensation bridge is recommended bj” which 
the major part of the impedance can be compensated. 

Although this procedure seems to be correct, it cannot compete 
Avith the methods described in the foregoing sections as to 
accuracy and ease of handling. If a telephone is used, the 
impedance differences are only of reasonable magnitude AV'hen the 
telephone is tuned rather sharply to resonance. It, then, is, 
hoAvever, very sensitive to changes of frequency, either of the 
audiofrequency generator or of the resonance frequency itself, 
and therefore, to temperature. Using a loudspeaker, temperature 
sensitmty seems to be of much le.ss importance, but noAV the 
atmospheric conditions influence the behaAUOur of the paper cone. 
This means, that Z^ , Z^ and Z^ cannot be measured once for all 
as a function of frequency, but mnst be measured for each 
Z-measurement. 

If Ave content ourselves, hoAvever, Avith the mere measurement 
of the absorption coefficient, rapid and accurate measurements 
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are possible from which the result can be obtained from two 
graphs without any calculation For full details the reader is 
referred to the original publications ’ ^ The outlines are the 
following 

According to (5 07) every point in the complex z plane corres 
ponds to a point in the complex plane of the quantity + z), 

so every contour in the 2 plane corresponds to a contour in the 
other plane If the curves, corresponding to the two sets of circles 
of constant absorption coefficient and constant phase jump in the 
z plane (see Fig 11), are drawn in the itf/CA + s) plane one can 
read directly from the latter graph the value of r belonging to a 
given point in the 3// (If 1- 2 )-plane without calculation This 
transformation is elucidated in Pig 36 for arbitrarj values of M 
and M being chosen real for simplicity, however If M were 
taken complex the circle-configuration in the 3f/(A + s) plane 
(lower part of the figure) should be turned through the phase angle, 
assigned to U, round the origin O The two sets of circles are 
transformed into two identical sets, apart from a translation, 
rotation, and magnification The circle of complete reflection is 
m the 2 plane the imaginary axis, in the other plane the circle 
with diameter OA All z points with positive resistance are found 
within this circle OA Shifting this circle OA over the complex 
distance L gives the transformation of the acoustical a plane into 
the electneal 7, plane 

Now, if sample and source are placed at either end of an inter 
ferometer tube and the sample is moved m the tube, the impedance z 
at the source runs through a circle of constant absorption coeffi 
cient So does the electrical impedance Z in its own plane too 
It was shenvn by one of the authors^ that the absorption coeffi 
cient of the sample bears a simple relation to the diameter of 
the circle in the Z plane belonging to the sample, the diameter 
of the circle for complete reflection and the position of the point 
for complete absorption in the large circle It does not depend 
upon the quantity £ of (5 07), 1 e, upon the centre of inversion 
Furthermore, the diameter for cmnplete reflection and the position 


a W Koaten, JppJ sc*. Ses B, 1 (1947} 8a 

D H Bekkering and C W Kosten, B,1 (I94S) 203 
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of the point for complete absorption are not critically dependent 
on frequency, in the case of a normal loudspeaker still less than 
wdth a telephone, and may be measured once and for all. There- 
fore the absorption may be read from a graph like Pig. 37, giving 
as a function of d/D and D/i, where d,D — diameter in the 



Pig. S6 

The relation between the motional impc<laiicc and the acoustical 
impedance (reaction interferometers) 

electrical impedance plane of circle belonging to the sample and 
a sample with 100% reflection resp., I = distance of the point 
of 100% absorption to the 100 %-refIcction-cirele (see Fig. 38). 
d is measured and D and I arc known from previous experiments. 
A suitable electrical set-up, with which the diameters can be meas- 

Sound ab-iorbinj materials ^ 
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urcd, IS shown in Fig 39 The current through the low-impedance 
loudspeaker is kept constant with the aid of a high series resist- 
ance The voltage on the loudspeaker is measured with a tube 



■ - 

Fig: 37 

The absorpcioa coefiicieot as a /uo«tioo of tbo <Iiasieter of cioetneal 
impedance circlee and the unpedance for 100% absorption 
(reaction interferometer) 



Fig 38 

The rneaniag of the qaantiUes D, d and I 
(reaction mteiferoineter) 
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§ 7 

voltmeter. This voltage can be compensated mth an auxiliarj' regul- 
ated voltage. After comi)ensation the length of the interferometer 
is varied till the tube voltmeter indication is a maximum. The in- 
dication is proportional to the diameter of the impedance circle. 

The results are reproducible’ udthin one unit in the second 
decimal of tto . Therefore, the same accuracy seems to be obtained 
as vdth other methods, whereas the apparatus is essentially simpler. 

Although this has not yet been published it is doubtless that the 
method can be worked out for impedance mea.surements, too. Phase 
angles can be measured if due attention is paid to the position of 
the sound source. 



Fig. 39 

Electrical set-ap for reaction measurement 
§ 7 GELEK’S niFEDAEOE lEDIOATORa 


During the last war a very remarkable and ingenious apparatus 
was developed by G e 1 u k. He used the screen of a cathode ray 
oscilloscope as complex a-plane, on which the impedance under 
test was indicated automatically as to magnitude and phase, and, 
since his apparatus worked instantaneously, the complete a-contour, 
when varying the frequency, could be shown in a few seconds 
and photographed if desired. 

A normal loudspeaker is pressed against the sample or wall 
under test (Pig. 40). The sound pressure between the cone and 
the sample is measured by a small condenser microphone, thus 
trausforming pressures into capacity variations. The displacement 
of the cone is also transformed into a capacity variation. This 
was done by coating the cone with a thin conducting layer and 
placing a second conducting cone 6, which remained at rest, at 

1 D. H. Bekkcring aud C. W. Kostcn, Applied sci. Eesearch, B, 1 
(19iS) 205. 

2 J. J. Gelnk, Thesis Delft, 1946. 
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a small distance of the actual cone of the speaker This second 
cone IS called the anti cone and is amply perforated in order not 
to hinder the vibrating cone in its vibrations The eapacitj -van 
ations proportional to pressure and displacement re^eetivelj are 
changed into electric -voltages with the aid of t%vo similar high 
frequency deiices The di^lacranent -voltage is multiplied b> } 
(turning the phase by 90®) thus providing us with a voltage 
proportional to tlie velocity of the cone The velocitj was 
automatically kept constant by controling the output v oltage 
of the generator with the aid of a sound level recorder of tlie 
Neumann type The pressure voltage is then proportional to 
the ratio p/v and if we succeed m indicating that voltage on the 
screen of the oscilloscope as to magnitude and phase an impedance 
indication is obtained In order to solve tins problem an auxilnn 



Pig 40 

T1 « Bound source mtl anti-eone of G « 1 u k a impedance ind cator 

V oltage of equal magnitude as the pressure voltage differing 
90° in phase was derived by electric means The two voltages are 
supplied to the two pairs of plates of the oscilloscope leading 
to a circular image on the screen of radius proportional to tlie 
pressure i e to the impedance z Finally the electron beam is 
sappT'Sssi^ arvih erespS /at tba mnnient the 

velocity passes through zero in the positive direction So one 
single spot is indicated on the screen giving the impedance as 
to magnitude and phase 

The disadvantages are obvious The accuracj can hardlj be 
expected to be as high as with interferometers. Jloreover the 
apparatus is ver> complicated and needs calibration as to tlio 
magnitude In those cases however where rapid comparison of 
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impedance contours is wanted, e. g., in developing materials or 
when checking the properties of an absorbing layer in a completed 
building, the method may turn out to be valuable. 

Before entering into some of the details of the apparatus 
something should be said about what is actually indicated on the 
screen. The velocity indication is obviously proportional to the 
mean velocity of the cone and does not depend upon the velocity 
distribution along the cone. The condenser microphone measuring 
the pressure is placed quite near to the cone surface, so the in- 
dicated ratio p/n is approximately the impedance at the source. 
The volume between cone and sample may be treated as a short 
interferometer tube of effective length I equal to volume divided 
by surface. Since the influence on the impedance of a column of 
air is a clockwise turning of the impedance along a circle of con- 
stant absorption (Pig. 11), the change in phase jump A being 
2 - for everj" half wave length A, the correction due to the effec- 
tive length Z is a counter-clockwise shifting of z along a circle 
of constant absorption, the correction in A amounting to 2 ~l/^ A. 
One sees that this correction decreases in importance with decreas- 
ing frequency as was to be e.xpeeted. Moreover it does not affect 
the absorption coefficient. 

§ 8 DETAILS OF GELUK’S lAIFEDAXCE IKDIOATOK 

Both sound pressure and velocity are measured through the 
intermediary of capacity variations. The latter are transformed 
into current variations by a modification of a method described 
by Z a k a r i a s Since this method deviates considerably from 
the commonly used R i e g g e r scheme ^ it may be outlined here 
(Pig. 41). The varjing capacitj' Cm forms part of the capacity 
of a high frequency LC-circuit in a high frequency generator using 
the oscillator part of the Philips octode AK2 as generator tube. 
To the fourth grid of the AK2 a circuit is connected tuned to 
the frequency of the generator with Cm hamng its equilibrium 
value. The Philips tube AP 7 is used as amplifier. Calculations 
show that the plate current of the generator tube depends upon 
the capacity of the LC-circuit according to Pig. 42. The apparatus 

’ Zakarias, Tech. Mitt. Tungsram, 39 (1938) 103. 

2 H. Riegger, JTiss. Veroffentl. Siemens-TTerhen, 3 (1924) 67. 
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1 * Used m the •'trsught part of the curve thus, avo diuff distortion 
to a hurh deirree The slope can be shown to depend upon fns 
qnenev and constants of the apparatus that can easily he deter 
mined, thus providinn us of a means of measurms relative capacitv 
variations in an absolute way 



Fig 41 

THe of caD T%nano“» 


Let Us first consider the capacitv variations due to the iao*ion 
of the cone The^ capacitv camtioos are proportional to displace- 
ments of the movable part of the condenser The proportionalitv 
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Fi^ 42 

Cnire”* cbaractoii:*ic to Fig 41 


constant can be computed from measurements of the oquilibniim 
value of the vamnir capacity and that of the additional condeo-er 
in the LC-eircuit So the amplifier ?ive& a voItaire proportional 
to the displacement which is absolutely cahbrated This voltace 
Is transformed in the normal wav into a voltace ii, proportional 
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to the velocity with an BC-conneetion (Fig. 43) in which RC = 10“=^ 
giving an acceptable transformation up to 1000 Hz. 

The calibration of the condenser microphone in a relative way 
is usually performed by applying an alternating voltage of the 
same frequency as that of the sound to the condenser plates, 
superposed upon a constant voltage: u — Ug + UiCosiat, in which 
« 2 i.g. The attractive force between the plates varies vnth the 
square of the voltage. The equivalent sound pressure is therefore 

r 

fk 2 cos (lit. 

The term in cancels out since it is constant, whereas the term 
in -iq- may be neglected ; the Kelvin constant fk in this equation 
does not depend upon the frequenc 3 ’. The alternating plate current 
of the generator now takes the foim 

i — fmfk • 2 itfliq cos iof. 


Fig. 43 

Turning the phase by 90° 

f<a being a function of frequencj' depending upon the mechanical 
b6ha\'iour of the microphone and the electrical behaviour of the 
measuring apparatus. Measuring the amplitude of i as a function 
of G) with Uo and tq constant gives tis fa apart from a constant 
factor. For the absolute calibration measurements were made of 
the change in direct current caused bj^ a verj' small static excess 
pressure, which, again, was measured with a sensitive mici’o- 
manometer. 

The pressure voltage is applied to one pair of plates of the 
oscilloscope. The same voltage, shifted through 90° in phase, is 
applied to the other pair. This voltage is made with the aid of 
the arrangement of Fig. 44 (see also the vector diagram Fig. 45). 
At everj' frequencj' R{a) should have the appropriate value; it 
varies inverselj' proportional to the frequencj'. This was pcrfoimed 
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automaticallj b^ coupbng the shaft of the \ enable condenser of 
the audiofrequency generator to the shaft of the metal wire 
resistance that was constructed espeeiallv to fit the purpose 





Fig 44 

TurniBg the phase by 90* independent tjf frequency 
without changing the mognitude of the output eoltage 



Fig 45 

t ector dregram to Fig 44 



Pig *6 

Impulse generator 


In order to suppress the image on the screen, except for one 
single point, the sj-nchronous pulse generator of Fig 46 was made 
The two thjratrons EC 50 diride the dc Toltage of the rectifier 




Periodic poises; 
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The first (left) starts when the input ^oltage just becomes 
positive charging the condenser m its anode circuit in a leiy 
short time thus gumg a short pulse to the output This situation 
lasts until about V* of the period later the second thyratron 
discharges initiated by the increasing negative voltage on the 



Pig 48 

Ke nforcaag the image on the screen mth pulses 



Fig 49 

Functional vheine of Oelolcs impedance indicator 


grid of the EBC3 This brings the arrangement back to the 
initial situation Form and speed of the pulses may be seen from 
the photographs Fig 47 a h and c 
The pulse is applied to the "Webnelt ej Under of the oscillo 
scope thus making the image brighter at the moment of the pulse 
A functional diagram of the complete set up is given in Fig 49 



CHAPTER VI 


EXPERIMENTAL RESULTS; COMPARISON VHTH THEORY 
5 1 IXTEODIICTIO^' 

In this chapter a summarj' will be given of measurements carried 
out on characteristic samples. We should like to emphasize the 
great value of plotting the results in the complex z-plane, as %ve 
shall do in this chapter, as a rule. The behaviour is given by 
the z-contour, to which a frequency scale is added, generally in 
hundreds of Hz. One reads at a glance from the graph the fre- 
quencies of resonance and anti-resonance, the degree of damping, 
the absorption coefficient as a function of frequency; moreover 
irregularities in the sample manifest them.selves directly in irre- 
gular contours. Plotting real and imaginarj- part as a function of 
frequency, as frequently done in the literature, turns out to be 
less sun-eyable. 

A difficulty encountered in all impedance measurements is a 
certain vagueness as to the phase jump A. The resulting A depends 
upon what is defined as the front surface of the sample. With 
a rough front .surface there may be an uncertainty of 1 or 2 mm 
in this respect. If one chooses the front surface more outside the 
sample by an amount Ax, the measured phase jump decreases by 
the amount 2- (2 Ax/X) which is the angle over which the phase 
shifts in passing twice the distance Aa:. In the impedance diagram 
this means a cloek^vise shifting along a circle of constant absorption 
for all measured impedances. Especially for impedances high 
compared vdth Wo such a .shift in A may mean a considerable 
change in impedance and, therefore, in the apparent character of 
the contour. This should be borne in mind by studjdng contours 
in the z-plane. 

Since the concept of the structure factor J;, as introduced by 
the authors, is still not generally in use, it seems worth while to 
give firstly a few results demonstrating its raison d’etre, nature, 
and importance. We shall then give results of more practical 
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samples A\hich can serve for elucidating onr theoretical considera 
tions, in doing this v\e shall follow the same division as was used 
in the theoretical chapters Most of the results given have been 
obtained with an old interferometer, a prototype of Pig 33, whieli 
did not give very accurate results, especially as to the phase 
jump A As was said in Chapter I, we smoothed the curves when 
necessary without changing their character Beranek’s results 
(§ 4) are of course the original ones 

} 2 EXPERIMENTS ON ARTIFICIAL SAMPLES 

Experiments have been carried out on samples made of glass 
tubes stacked in the manner of a honeycomb, with the axis per 
pendicular to the front surface of the model In this case 
and the resistance can be expected to agree wnth Kirchhoff’s 
theory Placing the same tubes ineJmed with respect to the normal 
towards the front surface provides us of a sample with a structure 
factor k greater than unity If the inclination is 60®, we must 
expect A. = 4 Using the same tubes the thickness of the sample 
IS two tunes smaller The sound veloci^ in the direction of the 
normal is, however, proportional to Ka, i c , is decreased bj a 
factor 2, too Therefore the same resonance frequencies should be 
found This complicated reasoning may be omitted, of course The 
resonance frequency of a single tube does not depend upon its 
inclination However the result is that the resonance frequency 
for a given thickness of the sample depends upon the inclination, 
It IS explained vnfhont having recourse to the concept of an 
increased air density “by vibrating fibres” as has repeatedl} 
been done 

Fig 50 gives a view of the sample with slanting tubes (dimen 
sions of sample approx 20 X 20 cm*, length of tubes approx 
20 cm) The impedance contours arc given in Fig 51 The ab 
sorption as a function of frequency is shown in Fig 52 The 
resonance frequencies are about the same One may casilj verify 
that the impedance of both samples onlj differs bj a factor J/I 
1 e , 2 in our case, the slanting tubes giving rise to the greater 
impedance This may be seen directlj from the fact that a close 
packing of tubes with 60® inclination onlj contains half as manj 
tubes as a close normal packing In a formal way this maj 
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I IS K times smaller Pjg 51 is quite in agreement lierewith 
Quantitatively the value of the resistance was greater than accord 
mg to Kirchhoff's theory by about a factor 3 This turned 
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out to be -due to the contribution of the interstices befrvveen the 
tubes, as further experiments have shown. For further details we 
refer to the original publications^. Another interesting sample is 
shown in Fig. 53. It consists of five layers of perforated “non- 
porous” materials (w'ood), 4 mm thick, at mutual distances of 
4 mm. The holes are 3.8 mm in diameter; 232 of them are made 
in the sample of 20 X 20 cm^. Part of the interstices between the 
planks is occupied by the supports of the construction. Taking 
this into account the geometrical value of k (i. e., the total volume 
of air divided by the volume of the canals passing through the 
planks) was computed to be 7.84. The resonance-frequency was 
750 Hz. Taking into account the correction for the open end at 
the entrance of the canals, this points to a sound velocitj' of 
4 • 0.041 • 750 = 123 m/see, gi^dng a A-value of (348.5/123)' — 
8.0, agreeing very satisfactorily with the geometrically computed 
value 7.84. Again, vibrating fibres cannot be made responsible for 
this high value of k. The origin must be sought in the structure 
itself. 


i 3 IMPERVIOUS MATERIALS OF HIGH ELASTICITY 

The best material of this kind is cellular rubber. The results 
have already been given in Chapter I § 7 (Figs 12 and 13). The 
density of the material was 116 kg/m^, the static rigidity modulus 
G — 0.17-10^ newton/m^ (0.17 kg/cm'). From the resonance fre- 
quency one computes 3.3 • 10^ for the djmamie value of the total 
stiffness. Since the air content is circa 88 % (derived from the 
density 116) the air stiffness is about 1.4" 10V0.88 = 1.6 • 10'. 
Therefore the dynamic rubber stiffness is 1.7 • 10^, i. e., 10 times 
greater than G. Assuming arbitrarily Poisson’s ratio wi = 4 
gi%-es Aistat./G = 3, leading to Aj dyn./ITi «tat. = 10/3. 

^Materials of this type are obviously verj^ bad for normal use 
owing to the strong selective absorption. Only in cases w'here 
selective absorption is wanted should they be employed. 


1 C, Z^vikker, J. vandenEijk, and C. W. K 0 s t e n, Phi/sic<i, 8 (1941) 
1094; 10 (1943) 239. 
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} 4 POROUS IfiraBIALS ^VTTH RIGID FRAME 

"With a vie^\ to the theory of elastic layers it needs verification 
whether materials like felt, rockwool, wood fibre plates, etc 
may be regarded as rigid, smce their air resistance is iisualJj 
considerable Acoustic plasters are of course good examples of 
rigid materials The contours of the former flexible materials. 



e^en that of soft sponge rubber layers, with porous front surface 
(see § 5) are however, of the same kmd as that of plasters 
A typical contour for an acoustic plaster is given m Fig 54 
The theoretical cu^^es show clearly the necessity of a structure 
factor of about 4 One sees that a higher air resistance is ^ery 
likely to increase the absorption coefficient at least at the lower 
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§ 4 

frequencies. Further examples are given in Fig. 55a for hair 
felt (2 cm thick density 100 kg/m^; K os ten), Fig. 55b for 
J. M. Permaeoustie (one inch; Beranek ^), Fig. 55c for 
J. M. Aeoustex (0.88 inch; Beranek') and Pig. 55d for 
J. M. Sanaeoustie (Beranek'). Beranek ’s results can be 
descidbed theoretically, as was shown by Morse, B o 1 1, and 



ah c d 

Fig. 55 

Impedance contours of porous “rigid” materials: a hair felt, 

6 JT. M. Permaeoustie, c J. M. Aeoustex, d J. M. Sanaeoustie Pad; 
frequencies indicated in 100 Hz 

B r 0 w n The values of porosity, resistance, and structure 
factor, needed to meet the measurements are given in the legend 
of the figures. The agreement is veiy satisfactory except at low 

1 L. L. Beranek, J. Acousl. Soc. Avx., 12 (1910) 11. 

2 P. M. Morse, E. H. Bolt, and E. L. Brown, J. Acoii.st. Soc. Am., 
12 (1910) 217. 

Sound absorbio:: materials ^ 
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frequencies At low frequencies the structure factor should, as 
a rule, be taken higher This indicates that, at least m these 
materials part of the air contents can no longer take part in 
the djTiainie process at frequencies of, saj, 800 Hz and higher 
B e r a n e k also gives the impedance contour for Celotex C — 4 
(1 25 inch) Here the general theoretical expression for homo 
geneous materials fails which is not surpnsmg (see § 8) 

Fig 56 finally gives us an idea of the absorption characteristics 
of this class of materials Low absorption at low frequencies is a 
general feature 



{ j POROUS MATERIAL ELASTIC FRAME OPEN FRONT SURFACE 

Fig 57 shows the absorption characteristics for three specimens 
of sponge ruhher For Laticel tiie air resistance iS obviously so 
low that no coupling effect is noticeable (<r « 20000 mks) , this is 
corroborated by the s contour which resembles very much tliat of 
hair felt (Pig 55a) Quantitativelj the contour can be described 
with the static value of o and a structure factor of 7 The air 
resistance of the Dunlc^iIIo sample is about 10 times higher, 
giving rise to more or less coupled vibrations and the accompanjing 
selective absorption at the frequencies of resonance of the ‘ frame 
■wave” 
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It is a well-kno-ivn fact that the absorption coefficient as 
determined in the reverberation room is, as a rule, higher than 
that obtained at normal incidence in an interferometer. In Fig. 58 



Absorption by sponge rubber layers Trith open front surface 



Fig. 5S 

ilrasurcments in the interferometer a.« compared to those 
with the reverberation method 


a comparison is given. The reverberation results tvere obtained 
in a room of 165 m- with 10.3 m- Dunlopillo on the floor. The 
dotted line is computed from the interferometer results bj* 
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assuming the impedance to be real (which certainly is not the 
case as a matter of fact) and independent of the angle of 
incidence and averaging the absorption coefficient for random 
incidence 

} 6 POROUS MATERIAL ELASTIC FKAilL COATED SURFACE 

Figs 59 60 and 61 show absorption characteristics of similar 
lajers as described in § 5, this tune coated with a thin impenious 
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Fig 6d 

Absorption by sponge rubber (l&ticel) layers witli coated surface 



Absorption by eoaied sponge rabber layers 
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layer of 0.1 kg/m^. These materials have been developed on a 
laboratorj' scale only. The coating layer vas of the cellophane 
type. They may be considered, however, as typical of the class 
of materials owing their beha-viour to high elasticity and porosity. 
The influence of the coating can be seen from Fig. 62 (smoothed 
curves). 

In order to obtain a more practical material a laj'er of 
Dtmlopillo was coated with a thin impeiwious layer of cotton, 
latex, and paint. This layer is rather hea%-y and \vill impair the 
absorption at high frequencies. Careful reverberation measure- 

to 


ao 


1 


ments with 8.4 m~ of this material in a room of 165 m“ were 
carried out the result of which is shoivn in Pig. 63. A pure 
tone was used ; every measuring point is an average of 
100 measurements. For comparison a small sample was sub- 
sequently measured in the interferometer (Fig. 64). The general 
.^hape of both cuiwes is the same. In accordance with the 
general practice the reverberation method yields higher absorption 
coefficients. 

The results of this and the foregoing section could be used 




1 C. W. K os ten and C. Zwifcker, Fhynco, 8 (I9il) 947. 
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This impedance terminates the absorbing lajer of thickness d If 
wave impedance Tt and propagation constant y of the absorbing 
material are knovTi, the impedance at the front surface of the eon 
stmction IS readily computed from the general equation (104), 
ginng the impedance of a lajer of gi^en II and y backed bj a known 
impedance In general IV and y are complex, in which case eomputa 
tion IS tiresome It is however easj to find approximatelj the 
frequencies at vhich resonance occurs i e , where the absorption 
coefficient mil he maximum To this end vie neglect damping of 
the waves in the material, the impedance is then purelj imagmarj, 
resonance taking place when the impedance is zero The eomputa 
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Fig 6o 

Scheme lad eating the g^-neral Iio along whiA the resonance frequencies 
of a porous material on a layer of air can be computed 


tlou of the impedance maj be replaced by a construction along 
the imagmarj axis of the circle-diagram This will be elucidated 
with the aid of Fig 65 

The phase angle A of the complex reflection coefficient r is zero 
at the rigid wall it changes bj 2i- for everj half wave length 
through the air The phase angle in the air, just to the right of 
the material is therefore Z (4 t-/X) The impedance can be indicated 
at once in Fig 11 Now on passing from the air just into the 
material the impedance is unchanged the wave impedance makes 
a jump bj a factor K/h, however This means that the sets of 
circles of Fig 11 are multipbed bj ]^k/h The phase angle just 
inside the material therefore, may be read from the multipbed 
circle diagram from the known value of the imaginary impedance 
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or, what comes to the same thing bnt is much more convenient, 
we may divide the impedance by k/h and read the phase angle 
from the same circle-diagram. On passing further to the left in 
the material another contribution of the amount d (4 -A') is added 
to A, if A.' denotes the wave length in the material. We then find 
the phase angle A at the front surface which should be (2n + 1) - 
for resonance. Fig. 11 is merely used for a convenient explanation. 
The procedure may now be put in the form of simple formulae. 

The connection between an imaginaiy impedance z and the cor- 
responding A-ralue is read from Pig. 10 to be 

z = jlT cot 4 A. 

At the right boundary of material and air the impedance is 
continuous, therefore 

17 cot 4 Aniat. = 17o cot 4 Aa.r 

from which A^at. is readily computed since all other quantities are 
kno-ivn. Adding d(4-A') and putting the resulting A equal to 
(2 n + 1) gives as the resulting transcendental equation from 
which the resonance frequencies are obtained^ 

tan al/Co tan wd X^lc^ = h/ 

The influence of thin layers of air is sometimes veiy consider- 
able (Pig. 66) which explains the discrepancies frequently found. 

The same procedure may be used for finding the impedance in 
the ease that the damping is not neglected -. The result can then 
be obtained most conveniently by a combined graphical and cal- 
culation method. This 'vvill be elucidated for a two-Iaj’er system 
backed by a layer of air (Fig. 67). The impedance is constructed 
in the z/TT-plane, whereby it should be borne in mind that 7' 
changes on passing from one layer to another. The value of a/17o 
at the point 1 of the structure maj' be indicated at once in the 
2/17-plane, since A is known in this point (Z4'n-/A„). From points 
1 to 2 7^ is discontinuous corresponding to the 1 — 2-jiunp in the 


1 J. van den E ij k, C. W. Ko.‘=ten, and W. K o k. Applied id. nescarch, 
B, 1 (19i7) 50. 

2 C. W. Kosten, Xederlcnds Tijd.fchr. yatuiirk. 14 (194S) 309. 
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zlV/ plane Since W in. absorbing materials generally is complex, 
the point 2 obtains a certain real component Prom 2 — 3 the 
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s/lf-contour passes through a part of a spiral, the end-points of 
which part are related by 

from which must be calculated in order to find point 3. The 
jump 3 — i in the c/IF-plane again corresponds to the discontinuity 
in IT at the boundary 3 — i, etc. 

Experiments with the sample on a layer of air have also been 
carried out with a special Laticel t.sponge rubber) sample with 
coated surface. In order to provide free passage to the outflowing 
air the sample was ampl}" perforated. In Fig. 68 cuiwe 1 gives 



Fig. 68 

Sponge rubber sample with coated surface on a 78 ram layer of air: 

1 free hanging, 2 on a sheet iron perforated support (Fig. GVa) 

results of the free hanging sample. In order to meet more 
practical conditions it was then supported by the sheet iron gauze 
of Pig, 69a, yielding curve 2 of Fig. 68. The resonance frequencies 
are obviously shifted, explained by the radical change in the back 
side situation of the sponge rubber frame (free end — v fixed end). 
To eoimteraet the selective absorption the same sample was stuck 
to a support of Pig. 69b, yielding Pig. 70. Owing to the irrcgularity 
of the sample plus support no resonance bands can be distingui.shed 
in this ease; the measuring points suggest a very complicatwl 
absorption .spectrum to be present. The average absorption is, 
however, verj' satisfactory. 
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The behaviour of Celote'^ C 4 deserves a special treatment 
The impedance contoui for this material as given byBeranek‘ 
IS shown in Fig 71 Jlorse et al showed the impossibiUtj of 
e\plaining this curve with the theory of homogeneous porous 
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media. The .stracture factor (their “effective density”) should 
be taken rapidly decreasing with increasing frequency in order 
to meet the measured contour. Remembering the influence of side 
holes on the structure factor this is not surprising. The behaviour 
may be understood as follows. 

The main pores, in which the sound propagation takes place, 
are the macroscopic holes of 6 mm diameter distributed in a 
regular pattern all over the material. The rest of the material 



Fig. 71 

Impedance contour of Celote.x C — 4 (Beranek, frequencies in 100 Hz) 

is fibrous and porous; the air contents of this porous mediimi 
takes the place of side holes to the main pores, which are only 
aeces-sible through a certain resistance, viz., the resistance of air 
flow through the porous material. Now this resistance is veiy 
likely to be so high that at high frequencies there is no longer 
a possibility of air exchange between main pores and side holes. 
In Fig. 72 the approximate electrical analogue of a small part 
of a main pore is shoum. The whole material is to be considered 
as a long electric line with distributed self-inductance (mass in 
the main pores), distributed capacity (compliance in the main 
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pores) and distributed capacit% with series resistance (eomplnnce 
3f the side holes with entrance resistance) The ratio Cj/C, is 
qim ilont to the structure factor of the mam pore* i e , total 
iir contents duided air contents in the mam pores It is of 
tlie order 15 to 20 This agrees fairh well with the measurements 
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Electneal an-^logoti to CclotfT C— 4 


In order to meet the measured cune Alorsc ct al assume A=s3S 
at about 200 Hz 12 at about 2000 Hz Prom the icsonance fre 
quencies one computes 1C at about 700 Hz and 2 5 at about 3o00 Hz 
Calculations haae shorni that indeed an element ns shown m Fig 72 
satisfies the measurements reasonably 
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ABSORPTION BV RESONATORS 

$ 1 IXTRODUCTIOX 

In this chapter we shall deal -vHth the absorption by resonators 
of the Helmholta tj’pe. A great number of publications * on 
this subject have appeared. Especially in the case of the single 
resonator the confusion is astonishing — the more so since 
Rayleigh has already given aU fundamental knowledge that 
is wanted — so that a treatment of this problem at some length 
seems worth while. The case of perforated panels before a rigid 
wall can be treated with the theory of resonators. It is much more 
important and simpler than that of the single resonator and is given 
in a short cut by W. W i 1 1 m s ^ and V. L, Jordan^ {Akust. Z.). 
The development of adequate constructions for any particular case 
has remained, however, a question that belonged to the domain 
of the acoustical expert. In this chapter graphs will be given that 
enable one to make a good choice, even ^vithout understanding 
the principles underlying them. 

In principle a resonator consists of a vessel with volume V, the 
interior of which is in communication with the external medium 
through a channel, which may be narrow or rather wide, short 
or rather long, circular in cros.s-section or not. If an external 


’ Rayleigh, Theory of sound IJ, 2nii ed., p. 170, Cambridge 1896; 
E. C. IVente and E. H. Bedell, Bell System Teclm. J., 7 (192S) 1; 
E. tVi n te r g e rs t, Z. teeU. Pliys., 16 (1935) 569; Rschevkin, Tech. 
Phys. USSB, 3 (19,36) no. 6; E. Meyer, Elel:. Kachr.-Tech., 13 (1936) 95; 
tv. Zeller, Ahust. Z., 3 (1938) 32; W. ’VVillms, Ahnst. Z., 4 (1939) 29; 
V. L. Jordan, Ahust. Z., 5 (1910) 77; P. O. Pedersen, Ingenidrri- 
denskab, Skrifter, no. 5, Copenhagen 1910; P. V. Briiel, LydisolaUon og 
BumahusUJ:, p. 104, Goteborg 1916; F. Ingerslev and A. K. Xielsen, 
Ingenidrvidenskab, Skrifter, no. 5, Copenhagen 1944; P. V. Bruel, Summer 
Symposium of the Acoustics Group of the Physicol Society 1917, London 1919; 
A. Gigli, J. Acou.ft. Soc. Am., 20 (1918) 839; V, L. Jordan, J. Acoust. 
Sac. Am., 19 (1947) 972. 
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source of sound provides a sonnd pressure at the aperture of the 
channel, the air in the channel will he set into vibration, by uhic3i 
the interior is periodically compicssed and expanded and acts 
like a spring The kinetic energy is concentrated in the region of 
Mgourous motion, i e in the channel The potential eiicrg}, on 
the contrary , is to be sought m the intenor The sj stem is closch 
analogous to a mass m mounted on a spimg of compliance C (dis 
placement/force) 

The “mass” of the resonator is that of the air m the channel 
it vibrates approximately mth the same \cloeit> throughout the 
whole channel If the channel has the length 1 and the cross 
section S', the mass amounts to 

tti = />a IS 

The compliance is computed from the definition 
C = dx/dF, 

if dx js a small displacement of the air in the channel due to 
a small driving force dP Non 

di = — dlVS, dP^Sdp K^ = —Vdp/dV 

if dV denotes the volume of air added to the internal air, dp the 
excess pressure inside that is approximately constant throughout 
the vessel, and Kq the compression modulus of air Therefore 


The resonance frequency of the resonator is found from the 
fimiliar formula to be 

2 » =c.|/^ (7 01) 

if G = S/l and Co = VaVpo= the normal \clocitj of sound in 
free air 

The channel will seldom be of such a simple form that it can 
he treated as a cylinder with a certain length, and c%cn m this 
ease our treatment is only an approximate one, since there is also 
a contribution to the kinetic enci^, and, therefore, to the effee 
ti\c mass, at the beginning and the end of the clnnnel, nhcro 
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the flow of air •widens into the free air and the interior of the 
vessel. Therefore we have to add to Z a small length due to these 
end corrections. For channels mth circular cross section the total 
correction for both ends amounts to 0.8 times the diameter D. 
If I is reduced to zero by taking as “channel” a single aperture 
in the thin wall of the vessel the mass is even completely governed 
by these end corrections. In this case G is exactty equal to D. 
For full details we refer to Rayleigh^. 

The equation (7.01), giving vres. as a function of G and 1^ is 
illustrated in Fig. 73. Since for holes in a thin wall G xD, 



Fig. 73 

The resonance frequency of a resonator as a function of conductivity G 
and volume V. For a circular hole in a thin panel G is equal to the 
diameter of the hole (see Fig. 83 for thick panels). 

indications of ID are also given in Fig. 73. One should bear in 
mind, however, that this is only an approximation for thui walls. 
Walls are thin in this respect if the thickness is veiy small in 
comparison with the diameter of the hole. 

A most valuable picture is the following. G — S/l is obviously 
the electrical conductivity of a conducting body of the same shape 


1 Rayleigh, Theory of scntnd 11, 2nd ed., $ 5 303 a. f., Cambridge 1896. 


Sound absorbioc materials 
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as the channel, if the specific condiictiMtj of the conductor is 
taken umtj. This scheme is of quite general validitj and can 
help us to estimate G in complicated cases If the %essel is per 
forated at se^e^al places at eonaderable distances apart the con 
ductivities G can simply be added Taking 4 apertures of dia 
meter I? m a very thin Tessel gives t? = 4i? If the apertures 
are as near to each other as possible thej mutually affect their 
end corrections, since the passage to the free air is restricted, 
thus decreasing the conduetivitj of each aperture The minimum, 
however, is 6 = 2 D, the eondiictivitj of one aperture of a surface 
of the four apertures together 

For a c^hnd^lcal channel of nearlj circular cross-section we 
ha\e approximately, as mentioned abo^e 

G = S/(1+0SD) (702) 

Another rather simple case for which the cond^ct^^lt^ can be 
given IS that of one single slit of breadth b in an infinite panel 
of thickness d (Pedersen*) Per unit length of the slit we ha%e 

1/G = d/6 + 0 7 + (2/-) In 2 (7 03) 

Even when the end corrections in 6 are taken into due account 
our description of the phvsical behaviour of the resonator remains 
an approximate one The kinetic energ\ is not restricted fulh 
to the channel and its direct environment, the potential ene^g^ 
in the channel is not completelv negbgible In fact, the description 
is based on the assumption that the wave length is ven great 
in comparison with all dunensions of the resonator, a condition 
which IS fulfilled to a sufficient degree with normal resonators 
with a Tien to other approximations that have to be made if the 
damping is artificially increased Thcee who are interested m a 
TeSiBessteiit ihe nsimJ dismuptira) ^iven above are 
referred to a recent excellent article of Uno Ingfird®, deahng 
with cjlindrical resonators with one circular perforation 


1 Ravleigh, TUory of touitd U, tUrf } S03 a.f , Cambridge 

* P O PederseiijIydtrttiwlvFiiffr^lJ^etseT.Zaso'K'wiidfiwjlflb Stn/lrr, 
p 72, Copenhagen 1^40, *e al^ Ingersler and N i el 8 en (note 1 page 12<) 
and H Lamb, HydrodynanniW, p 531, 6th ed., Cambridge 19S2 

* tJ Ing&rd, J Acwi'St Soe 20 €65 
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It will turn out to be necessarj' to provide the resonator %vith 
a reasonable damping in order to avoid strong selective absorption 
at resonance only. If this damping is performed by filling the 
volume with some porous material we must bear in mind that this 
influences both the volume of the spring (since part of the volume 
is occupied by incompressible material) and the conductivity. The 
latter influence is linked up Avith the structure of the porous 
material. Displacement of the air through the porous medium 
(which takes place at the site of the inner end correction) is only 
possible through the main pores of the material, i. e. through a 
fraction h/k of the space, where h is the porosity and k is the 
structure factor (see Chapter I § 8). This means that the apparent 
air density in this region may be as much as 5 to 10 times greater 
than the normal value, which can be taken account of by sup- 
posing the inner end correction to be 5 or 10 times greater than 
according to the diameter of Hie hole. For thin panels this 
amounts to a 3 to 5 times smaller conductmty, leading to a 1.7 
to 2.2 times lower resonance frequencj^ than expected. In com- 
parison Avith this effect the increase of the resonance frequency 
due to the decreased volume Avill be small. 

A resonator has three functions acoustically: 

a. it absorbs energy by friction in or near the channel, 

b. it acts as a second source of sound due to radiation from 
the channel mouth and therefore is a means for diffusing 
the incident sound, 

c. it stores energy by its resonance, AA^hich energy is partly 
restored to the room after the incident sound has ceased, 
thus prolonging the reverberation time of the room under 
certain circumstances; this effect is governed by the rever- 
beration time of the resonator itseK, which, again, depends 
upon its damping. 

It vrill be clear that, if the resonator has little damping, it will 
only respond to frequencies to which it is exaetty tuned. Then, 
hoAvever, it aaIU be set into vigorous action and will absorb much 
energJ^ Its storing capacity is great in this case, its reA-erberation 
time long. Out of resonance it aaIU practically do nothing, so as 
an absorber it Avill be unsuitable for most purposes. 

In order to make it absorb in a reasonable frequency range tlie 
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damping must be matenally increased But then it is much less 
suitable for storing the energy and its reverberation effect is of 
no importance 

Since we direct our attention to the absorption properties we 
are in search of constructions that absorb to a reasonable degree 
at reasonable frequencies in a frequency range of reasonable 
breadth, and maj, therefore, neglect the reverberation effect for 
the room in which it v,nli be appbed 
We shall now proceed to the eompatation of the absorbing 
quabties of resonators and shall deal with the following cases 

a the single resonator in an infinite wall, 
b perforated panels before a ngid wall, 
c multiple resonators m panel form 

{ 2 THE SINGLE RDSONt.TOK D.* AX INTINTTE WALL 

Let US postulate a plane travelling wa%e, impinging nonnalh 
to a rigid wall of mfimte dimensions, m which one single resonator 
is incorporated of conductiviti G, volume T’ and resistance R 
(Fig 74) The resistance is a measure for the damping of the 
resonator and will be defined nest 


A 


F»g 74 

TSe * 11 ^ 1 ^ KsooAtfSr 

\s regards the absorption and the diffusing effect the resonator 
fully described by its impedance Z as a function of freqiiencv 
IS defined as the ratio of excess pressure at the site of the 
entrance p over the velocity of volume displacement t It can 
easily be expressed m terms of G, R and V, viz 
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This equation expresses the fact that part of' the external pres- 
sure is needed for overcoming the inertia of the air mass in the 
channel which can only move with friction, the remaining part 
of the pressure serving for the adiabatic compression of the 
interior. The system, therefore, resembles an electrical series 
connection. The inertia term for a cylindrical channel equals, 
apart from the end corrections, j<opJ/S and must according to 
Eayleigh he written in the general form jtapjG if the channel 
is not cylindrical and end corrections are taken into account. 
To be precise: if we write G for G/l in the expression for Z also 
the end correction at the outside of the channel is included; 
1 ) contains an inertia term owing to that external mass, i. e., p is 
not exactly the pressure in the mouth of the channel but the 
pressure in front of the external mass, say at a distance equal 
to the diameter of the channel before the plane of the mouth 
(point A in Pig. 74). The stiffness term of Z follows at once 
from the definition of the compression modulus that equals xpo, 
jo> ori^nating in a differentiation in order to derive the velocity 
from the displacement (electrical analogj' 1/jaC). 

Obviously R is the component of the sound pressure p at A 
which is in phase with the velocit}' of volume displacement 
divided by the velocity of volume displacement v. 

If the sound pressure p at A has an amplitude p, gimng a 
velocity of volume displacement of amplitude v, the absorbed power 
is 4 v^R, analogous to the electrical expression PR = 4 can 

be calculated from p by 

p = \Z\v, 

so that we only want the value of p in order to find the ab- 
sorbed power. 

Now this is the point that frequently gives rise to confusion. 
If a plane wave of amplitude p-, is impinging upon the wall, the 
pressure at A has neither an amplitude 2 p.- , as is found just in 
front of a rigid wall, nor p, , as assiuncd by several authors, but 
depends also upon the radiation impedance of the mouth. In order 
to find this sound pressure assxime a massless piston to be present 
in the mouth of the channel that can move without friction in 
the aperture. This piston is mrtually subjected to a veloeitj' of 
volume displacement v with the aid of an uuxiliar\j extcimal force. 
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Furthermore the incident plane irave Js supposed to be present 
The pressure at the outside of the piston consists of tiro parts, 
VIZ., 2p, due to the incoming wave and f due to the radi 
ation, if ZnA. IS the radiation impedance of the piston The pres- 
sure at the inside is simply Z\, where 2' obriouslj is Z minus 
the impedance due to the external mags (;upoA7„t /S) Therefore 
the condition for equivalence of pressure on either side is 

2p, — Zna.v = Z't, (7 05) 

the minus sign on the left hand side being due to the velocit\ 
being taken positive inwards into the resonator If t is given a 
■value fulfilling this equation there w no external force needed to 
maintain the motion of the piston and the piston may he removed 
•without, to a first approximation, changing the situation There 



Fig 75 

The electrical airalogpiie of the •logle monator 


fore (7 05) is the relation between the velocity t and the pressure 
of the incoming wave for a resonator without a piston. Ohviouslv 
Fig 75 is the corresponding electrical scheme The e m f 2p, 
acts upon the senes connection of Znt. and Z' This can be 
derived from the electrical theorem the current in any branch 
of any electrical network is equal to the current that is found 
through that branch if it were connected to the terminals of a 
source with an e m. f , equal to the voltage on the terminals oi 
the branch in the network if its impedance is increased to m 
fimty, and an internal source impedance equal to the impedance 
measured at the terminals of the branch in the netwoik, if the 
branch and all e m f ’s in the network are remov ed 

The acoustical “e m f ” m our case is the pressure p for Z' =» , 
1 e., 2p, The internal impedance of the source is ob'viouslv the 
radiation impedance, since this is the impedance that is found if 
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tie resonator is removed and the impinging wave is stopped. 
Prom Pig. 75 we see 


= + Z'), 


giving as an expression for the absorbed power 


abs. power = ^ 


2 Pi 

^r»d. + 2' 


|2 

1 Ke(^') 


if Re denotes “real part”. 

A better indication for the value of the resonator as an ab- 
sorber is found by dividing the absorbed power by the incident 
intensity ^ | Pi 1^/lPo , thus giving the total absorbing power A 
in m^ Sabine of the resonator 


,_4T7„Re(^') 

12„d. + • 


(7.06) 


^rad. consists of an imaginary part due to the external end cor- 
rection, which, added to Z' transforms Z' into Z, and a real part, 
which will be called the radiation resistance. The latter real part 
is, owing to the wave length being very great in comparison to the 
dimensions of the cross-section, equal to ^ 


Re (Z„d.) = 2 vlPoA" = Rrad. . (7.07) 


It is clear from (7.06) that we may replace Z' by Z, if the 
imaginary part of ^rad. is left out of consideration. 

A will be a maximum in the neighbourhood of resonance, i. e., 
when Z — B = vc&\. Varying R in order to obtain maximum ab- 
sorption at that frequency (i. e., at A. £= constant) we are obviously 
to take 

R Rrad. 

and find 

A„.a..=Ay„./2,7. (7.08) 

E. g., if the resonance frequencj’' be 345 Hz, i. e., A„s. = 1 m, about 
m= of incident energj^ is absorbed. 

Another feature of the resonator, its diffusing effect', may 
readily be estimated. Interchanging Z' and Zr^A. (7.06) immediately 


J Bayleigli, Theory of sound II, 2nd ed., } 319, Cambridge 1S9G. 
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gi\es the amount of power that is radiated from the month m all 
directions, measured in terms of the incident intensity (7 08) 
gnes the amount of diffused power if and the resonator 

IS in resonance The maximum of diffused power is howe\er 
4 times as great This can be shown by interchanging Z' and ^r.d 
in (7 06) and considering the case of resonance with R = 0 
The result of (7 08) looks very promising We should bear in 
mind however, that we generally want to absorb in a rather broad 
frequency range of, saj, two octaves at least If the resonator has 
so bttle damping as to give (7 0S) it absorbs strongly selectively 
and shows very little absorption outside resonance In order to get 
an insight into the effect of the resonator it is, therefore, necessary 
to study the absorption outside resonance 
The equation (7 06) also applies m this case If we introduce 

f =frequency/resonanee frequency = v/v,,, , 

fx = resistance of resonator/radiation resistance at resonance 
=:BXS„/2r IFo (soe 7 07), 
g = G/X„t , 


(7 06) can be shown to assume the form 


At resonance we have (/«=1) 


(7 09) 


2 7r (l+/i)*’ 


(710) 


which equals (7 08) for /*=! 

Outside resonance (/ 9 tl) A will be smaller than A„i A will 
be A„8/2 for 

cf + /t)* + (1 + (711) 

From this equation we can compute two frequencies /, md U> 
the number of octaves O between /i and /* being an adequate 
measure for the curve width 

There are five quantities that can be looked upon as measures 
for the efficiency of the resonator, viz , 
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-^ircs./A^res. = 3 . fulictioii of fi Only (sGe Pig. 76) that is maximum 
for fi — 1 (the part of the cumm for /i < 1 is not shown in 
Pig. 76), 

0 = the number of octaves between the frequencies where the 
absorption is Ar^J2 (see Pig. 77 and 78, dotted bnes) ; 0 is 
a function of /x and g, 



Fig. 76 

The number of m" Sabine (m- “open uindow”) absorbed at resonance 
by a single resonator, divided by the square of the ivave length at 
resonance, as a function of the ratio internal resistance/radiation 
resistance (P'1) 


Ares. 0/k^res. ~ 3 function of ft and g (see Pig. 77 lower part) ; 
it increases witli g, i. e., uith G; this quantity measures the 
effect of the resonator in an ahsohUe tvay. 

Ares. 0/0^ — a function of fi and 3 ; CP is a rough measure for 
the surface of the mouth (for thin vessels 6 — D); this quan- 
tity, therefore, is of some value aesthetically, since it gives 
the effect in comparison with the surface of the perforation 
of the wall; this quantitj* is not given in a graph, but it 
increases unth decreasing conductivity (diameter of the 
mouth). 
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^res. • 0 /Ares. 6^ = a functioii of [1 and g; at low values of jj. 
and/or g oi g. only (see the band of flat curves in the upper 
part of Fig. 77) ; this quantity is, like the foregoing two 
quantities, a measure for the absorbed area of “open -ivindow” 
times the number of octaves for which the “window is open”; 
the unit surface with which it is compared, is this time 
^res. G ; expressed in this unit its magnitude is approxima- 
tely 1 for g. adequate and g not too great. The latter fact 



Fig. rs 

Wall economy and selectivity of a single resonator as a function of 
internal resistance (relative measure /?) and conductivity G 


is illustrated in Fig. 79, a figure worth while remembering, 
since it tells us practically everj'thing. Ares. • 0/Ares. G is 
proportional to Ares. ' 0 ■ Kes./V, i. e., to the ratio of the 
effect Arcs. • 0 and the volume of the resonator and may for 
this reason be looked upon as measuring the effect from the 
economical aspect- since the quantity under consideration is 
nearly unity for p adequate and g not too great, we may 
eonehide that all reso7iators are equally econoynical {y ade- 
quate, g small). 
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O/Y^!^ = a hmehon ot n and g An esfimaljon of ffte effect 
of a resonator -with respect to the nail economy can perhaps 
best be gi%en bj expression the product open window X in- 
ter! al in units of the area of the wall occupied bj the reso- 
nator, IS a rough measure for this i-esonator surface 
(see Fig 78) 

Pig 78 shows that, in order to have the resonator ab'^orb with 
high wall eeononi! in a frcqncncv band of one octal e or more, 
p. should be 10 to 30, whereas G/k„, should be possibh 0 04 or 
more If the most economical resonator that is practicable is found 
from this figure, we can compute the eflect, or use Fig 77 for 
coincnience 



Pig <0 

The absorbed ar«i at resonance tim<* t)ie number of oclarc' 0 
for which tho absorbed area e\cted» fl is approximMely equal 
lo Oh , if i« Jmd ff are wasoonbje For a circular bole in a thin 
tancl G fa 0 


All curves m Fig 77 are straight for small value* of G/Xr„ 
Alathematical expressions for tlicm can easilv be derived from 
(711) assuming / lo be approximatelj 1 and g verv small m 
companson with 1 We obtain 

0=1442(1 +/i) <7 
0/A%„ =0 918w//(H- /i) 

A„. 0/G* =0918/t/ff (1 + p) 

A„. OA„, 0 = 0918,1/(1 + p) 

Ar« 0/7*'* = 10 64 ,1 l^'p/ (1 + /O 



1442 IS the decimal value of l/lo2, 0918 of 2/Tln2 
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It is surprising that these formulae hold -srith sufficient accuracy 
up to intervals of about an octax'e. 

5 3 THE HSE OF SDtGEE EESOXATOEE IN PEACTICE 

The data of § 2 wiU be elucidated by a practical example. Sup- 
pose a room of about 1000 m^ to have brick walls and a few 
windows. The room is to be used for meetings of about 20 people. 
The overall-boundary surface will be about 600 m-, the average 
absorption coefficient about 5 %. The optimum average absorption 
coefficient can be calculated from Sabine's law (saj" reverber- 
ation time wanted 1 sec) to be 28 %, so some 23 % has to he 
provided for. Eesonators are supposed to be prescribed for a&sthe- 
tieal reasons. 

In order to keep the costs as low as passible the most annoring 
frequencies, say 150 — 600 Hs, -vvill be most absorbed. Let us try 
to find suitable resonators of one single kind, resonance frequency 
300 Hz, interval 2 oetave.s. We calculate; 

= 340/300 = 1.1 m, 

G — S cm (high, but tolerable), 

= 0.027 (reasonable wall economy, Fig. 78), 
a;50 (Fig. 77 or 78 for 2 octaves wanted), 

J-rea. • 0/A%„. = 0.021-5 (Fig. 77 for jr = 0.027, /i = o0), 

= 0.0215 X 1.1-/2 = 0.013 m=, 
number of resonators wanted /m- = 23/1.3 = 17.7, 

V s; 1000 cm’ (eq. (7.01) or Fig. 73). 

There are several reasons for decreasing the number of resonators. 
The first one is that the theoiy of § 2 applies to normal incidence, 
whereas in practice the incidence is more or less at random. Xow 
the aeonstical “e.m.f.” is 2p, independent of the angle of inci- 
dence 6. The absorbed power is therefore also independent of 6, 
whilst the absorption coefficient of the wall is proportional to 
l/cosi&. At 0 = 60'^ the absorption coefficient is doubled, whereas 
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for 6 approaching 90° the absoiption coefficient increases with- 
out limit Even the effect averaged over the solid angle (see 
Chapter VIII § 2 equation (8 03) ) is mfmitelj high \Ye therefore 
think that m practice the effect will be at least a factor 2 greater 
than according to § 2 

A second advantageous eireamstanee is that the impedance eon 
tour of a resonator is not exactlj a straight vertical line in the 
complex impedance plane, as is assumed in equation (7 04) It 
IS sometimes curved, the resistance being approximate^ minimum 
at resonance This means that the resonance peaks are to be 
expected broader than according to the theorj’ of § 2 

It therefore seems justified to assume that in the case mentioned 
8 resonators per m* have a reasonable effect Of course several 
surfaces are unsuitable for mounting resonators (doors, windows, 
the floor), so that the number of resonators per m* of the remaining 
surface should be increased accordingly 

Although the solution of the practical problem of this § seems 
correct there is a reason to make another choice p. should be 50, 
which will be difficult to achieve without measurement of resis 
tances, e g bj impedance measurement An altemattve is to take 
about the same number of resonators but to use a certain spread 
m resonance frequencies /i can then be taken smallei This will 
even slightly augment the wall economy (Fig 78) "We then use 
resonators of ^=10 to 50 The effect will now be rather inde 
pendent of n 

Finallj we should bke to emphasize that m order to act as 
absorber the resonator should have a reverberation time much 
shorter than that of the room where the resonators are used In 
most cases this condition is easily fulfilled In case of doubt it 
should be checked Without giving the derivation we give the 
formula for the reverberation time 


Requiring t not to exceed 0 1 sec and putting n for reasons of 
economy at the minimum value 10, we can compute the minimum 
permissible conductivity (approx diameter of the mouth) as a 
function of the resonance frequency (see Table 9) At low fre 
queneies is seems advantageous to choose p. greater than 10 
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TABLE 9 

inELMXrM COAT>ECTIVIXT G AT TAEIOES EESOXAECE FEEQUEECTES 
FOE r < 0.1 SEC ANT) fi > 10 

. ^mia 

100 Hz 0.069 m 

200 0.017 

400 0.004 

800 0.001 

« i PEEFOEATEB PANELS BEFORE A RIGID IT ALL 

^Ye now proceed to the more important ease of perforated 
panels, mounted at some distance before a rigid wall and similar 
constructions that can be treated as a sj'stem consisting of an 
infinite number of resonators, aU channels ending in a rigid plane 
walL The resonators will be supposed to be distributed in some 
regular way orer the surface, the number per m- being n. This 
ease indeed is mbre important than that of the single resonator 
since it turns out feasible to make in a more economical way 
constructions that absorb verj- satisfactorily in a wide frequency 
range. As a matter of fact absorbers of this type are in widespread 
use. They usually contain some porous medium (e. g., glass silk) 
between panel and wall. 

Now from a theoretical point of view it is much easier to sup- 
pose the air resistance concentrated in the channels, i. e., near the 
perforations, which will, therefore, be assiuned presently. There 
is much evidence, however, that the theoretical residts wiU be 
applicable to practical constructions with a porous medium in the 
entire volume behind the perforated panels. Frictional losses will 
be maximum where the velocities are maximum, i. e., near the per- 
forations. The bulk of the porous medium is present in the regions 
where the air essentially is at rest and acts as a spring. This part 
of the porous medium only manifests its presence bj- reducing the 
air volume, thus increasing the resonance frequency a little. In- 
deed, the experimental results wiU be shown to agree fairly well 
with the theoretical formulae based on the assumption of concen- 
trated air resistance. 

The quantities governing the problem are (see Fig. 80) 
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V c= the volume per perforation, 
n = the number of perforations per unit surface, 

(? = thc conduetuity of one hole {S/l for a long cylindrical 
channel), 

B = the air resistance of one hole (pressure difference duided 
by \elocitv of \olum6 displacement through one hole dis 
regarding inertia effects) 

If the dimensions of V, and, therefore, the mutual distance of 
the holes, are very small m comparison to the ua\e length, a 
'specific acoustic impedance can be defined as z = prcssure/\elocit\ 
of \olume displacement per m*, from which the absorption coeffi 



t V ' 


Fig $0 

Quantities gownwng Uic behaTjour of perforated pMcls 


cicnt Oo of the construction can be computed wth the aid of the 
usual relation 


Oo=l — 


z-W, * 
2 + ir. 


(713) 


This expression can be obtained along the same lines as in the 
preceding section We then are to take the acoustical “emf" 
equal to 2pi, as before, the internal impedance or radiation im 
pedance being Wo It is clear that the external end correction 
has to be incorporated in z 

Since the vofume displacement per is % times (fiat of one 
hole we can put immediately (see (7 04) ) 


2 = -L [B + ,.p./c + (7 M) 

It is not advantageous to substitute (714) in (713) in order 
to find Cq ^Ve can better haxe recourse to the circle diagram 
(Fig 11), from which ean be read if z is known In resonance 
the impedance is R/n, a real quantity that can be xaned at will 
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for the same construction by adequately choosing the air resistance 
R of the perforations. If the frequency is varied around resonance, 
B can be considered as a constant and an impedance contour 
parallel to the imaginary axis is found. The frequency scale along 
this contour only depends upon n, G and V (see the imaginarj' 
part of z in (7.14) ) and is independent of R. One sees at a 
glance from Fig. 11 that, if R is increased for a given construction 
(?i, G and V constant), the frequency range of absorption is ex- 
tended. But, simultaneously, the absorption coefficient at resonance 
decreases %vith increasing B. Therefore, increased B means decreas- 
ed selectivity (in most cases a great advantage), but decreased 
maximum absorption coefficient (in most cases a disadvantage). 
The problem is, therefore, to find for every case the most adequate 
compromise between these two features. 

An adequate measure for the selectmty of the construction 
is the frequency range in which the absorption coefficient is 
more than 50 % of the maximum value (resonance). Putting 
z — r + jx in (7.13), where a: is a function of frequency, 
is obtained by inserting z — r (i. e., a; c= 0) . The critical value 
of X for which % is equal to is easil}^ found from (7.13) 

to be 

a:=± (r'-t- 1F„). (7.15) 

The critical a-value as a function of r can obviously be repre- 
sented in the complex impedance plane by the two lines at 45° 
slope through the point — 1F„ (see Pig. 81), from which it is 
clear that the selectivity is only small for not too small values 
of r (— B/n) . Taking 40 % and 75 % as reasonable limits for 
in practice, one sees from Pig. 11 and 81 that R/n should be 
somewhere in between 3 IFo a^d 8 IF^ , say between 1300 and 
3500 mlis-units. 

Inserting the expressions for x and r, as read from (7.14), in 
(7.15) gives an equation from which the two critical frequencies 
that determine the width of the absorption characteristic can 
be computed. As in § 2 the ratio's of these frequencies and 
the resonance frequency will be labelled and Thej^ can 
be called the relative critical frequencies. Bearing in mind that 
^res. ^0 VWv, one easily verifies that and /„ satisfy the 
equation 


Sound absorbing materials 


10 
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(716) 

TJiJs equation is the parallel of (711) of the single resonator 
theorj (716) is, hoi\e\er, much simpler All quantities except f 
are combined in one single variable the right hand side of (7 16) 
This is due to the simpler radiation impedance in this case, that 
IS independent of frequencj One simple feature of (7 16) is, that 
fJi — l, 1 e fi and fg are situated s^mmetncallj witli respect 
to 1 on a logarithmic frequency scale 
If the number of octaves 0 betueen and fg is small simple 
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To the selectivity of perforated panels 

■relations hold as in § 2 There is eten a close analogy in the 
formulae if analogous tanables are introduced 
Labelling 
E/nWo c= p.', 

(conductivity per m®), 

nY = (volume per m* « tndtii of the air gap between 
panel and tvall), 

G'A„. 

the following equations can be shown to be valid 
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§ 5 


ao,e,. ~ -i mV(1 + IJ-'y 

0 = 1.442 (1 + /) g' 

■ 0/X\,,. = 0.918 p.'g'/{l + pf) 

• 0/G' ^ = 0.918 p'/g' (1 + p') 

■ OA,es. 6^' = 0.918 /xV(l + p'), 

i. e., analogous formulae hold but for the factor \^res./2 that has 
disappeared from the equation of and now enters in 0. The 
equations in wliich the product • 0 is involved are the same, 
i. e., also in this case high values of p' and g' are advantageous. 
It should be borne in mind, however, that (7.17) only holds for 
rather small values of 0, thus bmiting the validity of (7.17) to 
certain maximum values of p' and o'. 

If 0 is no longer supposed to be small, graphs have to be 
used for exact design Avhieb this time are simpler however, since 
as Avell as 0 are functions of one independent variable only. 
The design is then based upon the formulae 



“o,„. =4/t'/(l 4- p'Y 

0 = a known function of the product Are8.G'(l + p') 


(7.18) 


The followng way of developing a construction seems suitable 


1. ohoose , 0, Arcs. > 

2. p' and Ares. G' then follow from (7.18) or from the graphs 
representing these formulae, 

3. p' yields R/n, 

4. Ares. G' and Ares, yield G' — nG and (7'A^reB./474, 
i. e., the gap between panel and wall, 

5. the amount nG can be achieved with manj^ small holes or 
with a smaller number of big holes. An approximate prac- 
tical solution may be obtained from nG ^ nD, whei’e D is 
the diameter of one hole. 


§ 5 RULES POE DESIGNING PEEPOBATBD PANEL ABSORBENTS 

In this section we shall give some graphs for the design of 
practical constructions, based upon the theoiy of § 4, but adapted 
to the use vuthout understanding this theorj-. 
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Fjg 82 shows what can be achieved with perforated panels 
If , 1 e , the absorption coefficient at resonance (maximuin 
^alue of Cg) and the number of octaves 0 in which the absorption 
coefficient is greater than /2 are given, Pig 82 immediately 
yields the necessary values of R/nW^ (right scale on the ordinate 
axis), the ratio WAre» of the thickness of the air space F' and 
the wave length at resonance Are, , and the product Ares nG = Ar„ G', 


^f.A 

1 2 4 6 8 1012H1&1820 30 40 50 



where n is the number of holes per m*, G the conductivity of one 
hole and, therefore, G' is the conductivity per m* 

Owing to the relation =CgV^'/T" the ratio F'/Ar«s is 

equal to Af^, G'/4 so that the curves of constant F'/Ares m 
Fig 82 are also curves of constant Ar«aG' (constant values approx 
40 times greater) 

In order to give another idea of the frequency band in which 
the absorption is reasonable also the abscissa axis has a double 
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scale. T.he main scale gives tlie number of octaves between the 
two critical frequencies and f^, the supplementary scale the 
ratio between these frequencies. 

As a rule the frequenej’ at which the absorption is maximum is 
also specified, so that from F'Ares. the air space V' between panel 
and wall can be computed. If the resonance frequency is 340 Hz 
Arcs, is one meter and V'/Xrea. is equal to the air space in metei’s. 
One sees that for 75 % maximum absorption in 3 octaves (150 — 
1200 Hz) V'/Xres. should be about 0.1, i. e., V' should be about 
5 to 10 cm. 

The numerical value of Xrcs.iit?, found from the requirements, 
ean be obtained in many ways. If the perforation consists of 
cylindrical holes of cross-sectional area S in a panel of thickness I 

G = S/(l+ 0.8 D), (7.19) 

where H — diameter of a hole. If the damping of the resonator 
is provided by a thin air resisting layer in or quite near to the 
channel this expression holds. If, however, the entire volume of 
the resonator is filled with porous material or the air resisting 
layer near the channel is of a thickness equal to the radius of 
the hole or more, G is decreased. We then have to replace 0.8 D 
in (7.19) by 0.4 (1 + k/k) D, Avhere k is the structure factor 
and 7i the porosity of the porous material. As a rule k/h may be 
assumed to amount to 2 for generally used filbng materials. If 
the panel is of small thickness in comparison to the diameter of 
the perforations, (7.19) reduces to 

G==D (7.20) 

or to 

G=2D/{1 + k/h) (7.21) 

if the entire resonator is filled %vith porous material. In practical 
eases (7.20) may be a rough approximation, verj* helpful in the 
design. If the thickness is not negligible and if k can be put equal 
to unity (air resistance concentrated in a thin layer or unfilled 
resonators) an adequate couple of I, 1) values can be read from 
Pig. 83, giving (? as a function of 1 and D. E. g., the conductivity 
G — 3 mm may be obtained with holes of 5 mm diameter in a 
panel of 2.5 mm thickness. 
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Example Dcsipi a construction -witli resonance frequency 300 Hz, alworbing' 
60 % at 300 Hz and 30 % at 100 Ha and 000 Hz 

From Fig 82 we find that, in order to fulfil the requirements (a„ =0 6, 

/1//2 = 9) we are to take 

^Vv. =0 03, 

a,„C' = 32, 

E/nW^ = iJ5 

Since =345/300=115 m we obtain 

= 0 00 m = 9 cm 
G' = «G = 2 8 

The air space la, therefore, fixed and ^ould be 9 cm A few solutions that 
fulfil G'=28 ia~* are 



G 

1 

1 ® 

I 

% perf 

250 m"® 

1 11 2 mm 

14 mm 

1 mm 

4% 



^ 20 

10 

S 

500 

50 

i 7 

1 

2 



12 

10 

55 

1000 

28 

35 

1 

1 



8 

10 

5 

2000 

14 

2 

1 

06 



5 

10 

4 


These rwults are directlj read from Fig 83 The last column gives the 
percentage perforation la order to Iccep this low many HnaJl holes in thin 
paneJa are required If aesthetic requirements do not plav a part a solution 
with a small number of big holes might bo much more economical and equally 
effective 

The real difficulty of the design is to adjust li to the required 
value by means of some practical construction This may be done 
by filbng the whole gap with some porous medium (e g , glass 
silk) of the required air resistance 

For thm panels (thickness small in comparison to the diameter 
of the holes) the resistance B can be expressed in the conduetiv itj 
G of one hole and the specific air resistance of the filling material 
This relation holds even for non circular perforation owing to the 
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close analogy of the two concepts. 1/G is, according to its definition, 

G~J S 

the integral being taken over the channel, end corrections included. 
One should bear in mind that the effective area -S' through which 
the flow of air takes place is fc//t times smaller if porous material 


-2 



Pig. S3 

The conductiyity G of a circular cylindrical channel 
of diameter P in a panel of thickness I 


is present at the site under consideration. R is, also according 
to its definition, 

R— f 

J vS 

ii v= specific velocity of volume displacement, the integral being 
taken over the same path. Now {dp/dl)/v is the specific air 
resistance v of the medium. If o- is great at the site of the internal 
end correction only, the B-integral can be restricted to that part 
of the channel; if the panel is verj' thin, the G-integral is taken 
over exactly the double lengtli, and owing to the symmetry 

K = a.l/Cf (1 4- /-r//i) (7.22) 
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Now G and R are found from the requirements, so a can be com 
puted if k/h IS put equal to, say, 2 

Therefore, if the whole gap is filled with porous medium, its 
specific air resistance should be 

a=(l + k/h)RG (7 23) 

Since resonators are generally used at low frequencies, the steady- 
flow -value of a can be taken as an approximation For higher 
precision the dynamical -value at the corresponding frequency must 
be used, which can be greater (see Chapter I § 9) 

With practical constructions the product RG = G’R/n only 
vanes in a Iinuted range (see Pig 82) from which values for <r 
can be computed lying in the neighbourhood of that of hair felt 
or lower (Chapter II § 5) If the panel is thick in comparison 
to the diameter of the holes c should be higher than according 
to (7 23) in the ratio 

{ i + 0 4 (1 + k/h) Z) } / 0 4 (1 + k/h) D 

A simpler way of adjusting R to the required value is the coating 
of the panel at the outside or the inside with a thin lajer of 
adequate air resistance, eg a high resistance textile coating, 
blotting paper, etc Now there is no objection against measur 
ing the resistance under steady flow, since such high resistance 
layers are certainly governed, even m the dynamical case, by 
the P 0 1 s e u 1 1 1 e law Moreover the quantities k and h do not 
enter into G in this case 

An advantage of filling the entire gap wth porous material is 
probably that we then obtain constructions that also absorb reason 
ably in the frequency range where the resonator conception fails, 

1 e , far above resonance 

J 6 RULES FOR DESIGNING PANEL A^OEBENTS WITH SLITS 

If the perforations are not circular in cross-section, but have 
the form of infinite or very long sbts in the panel the rules of 
^ 5 can he used as well We then have to consider V, G and R 
as the volume, conductivity and air resistance per m length of 
a slit A rough approximation for the conductivity G might 
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be computed from the formula for a single slit, according to 
Pedersen (see § 1 of this chapter) 

llG^d/b + 0.7 + (2/Tr) ln2 CoM, 

although there is some evidence that this expression lacks aceuraej" 
in the ease of panels with manj- slits. 

For thin panels and gaps completely filled with porous material 
we, again, have 

<r= (1 + k/h) RG. 

§ 7 EXPERIMENTAL RESULTS 

The literature does not provide us vdth many adequate data 
corroborating the theoiy. In this section results of a few funda- 
mental experiments will be described 

Impedance mea-surements have been carried out on simple reso- 
nators; they merely consisted of a circular metal plate of 5 mm 
thickness, provided with a circular hole of 4 or 8 mm diameter 
in the middle, which plate was fixed in the sample holder of the 
interferometer of Fig. 33 (diameter of the tube 90 mm) at an 
appropriate distance before the rigid back plate of the holder. In 
this way resonators of tw'o conductivities and adjustable volume 
could be built. 

Experiments of this kind yield results that should comply vrith 
the theory of perforated panels (§ § 4, 5), not ndth that of a single 
resonator (§ 3). As a matter of fact the reflected wave is plane 
at a reasonable distance from the resonator, as is the case with 
perforated panels. The behaviour is the same as that of a wall 
completely occupied by resonators of the same kind, i. e., with a 
number n per m" equal to 1 divided by the surface of the inter- 
ferometer tube in m^. It makes no difference that the interfero- 
meter tube is circular in cross-section. It is merely not feasible 
to cover a wall n-ith n round resonators of the described form. 

In Fig. 84 three impedance contours are given drami in full, 
the data of which are 

i? = 4mm, F' = 1.8mm, 8.1mm, 11.3 mm resp. 


1 Wc thank Mr. E. G. dc Brnyn, phys. cand., for carrying out the 
experiments. 
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that the agreement is very good except for the resonator of 2 mm 
depth {V' = 1.8 mm as a consequence of the thickness of the 
sample holder). Obviously, if the resonator has a depth smaller 
than the diameter of the hole the conductivity is lowered appre- 
ciably, as was to be expected^. Another feature of the impedance 
contours is readily explained. The resistance (distance from the 
imaginarj" axis) is materially increased by the vicinity of the 
bottom of the resonator if the depth is of the order of the dia- 



Absorption characteristics to Tig. 84 

meter D. The curvm width in octaves turns out to be in fair 
agreement with theory (Fig. 82). We may conclude from these 
results that practical non-filled resonators comply noth the theorj’. 

Similar results are obtained (see Pig. 84 and 85 dotted lines) 
with 

D — 8 mm, Y' — 1.8 mm, 7.2 mm, 16.2 mm resp. 

res. freq. measured 870 Hz, 515 Hz, 353 Hz resp. 

res. freq. from H and V 1063 Hz, 532 Hz, 354 Hz resp. 


' See also U. Ingard, J. Acousi. Soc. Am., 20 (1S4S) C63. 
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The «mdnetiTitv api>eais to be lowered <;ome 7 even if the 
depth li equal to the diameter For re^uatoi^ of normal depth 
the theorr hohL. verr well The number of octaves is amun ap- 
proAimatelr as theoretieaUr predicted from <i# and VVA,^ 

In. Fisr S6 results of more practical re^nator^ are jnvm* The 
resonator was a thm perforated metal plate phced in the inter 
le’^meter as prerjouslr The air -space behind the perforated plate 
was SO mm tie cap was filled with glass sitfe-. Lnfortnnateir no 
impedance^ were measured. The perforation cons.s’ed of circular 
hole^ in s “cubical lattice ,56 mm apart from each o*her, 4 mm 
m diameter, the whole plate had about 175 holes. Variations of 



ns 

Alwiyj*! a «iarsclenstic» o* the w'vaater ^tJi d.ffe’e- p^rfera i •- 

the leoaator were made bv clcsinsr a great number oi thc'C ho’es. 
The eurres m Pi<» S6 from left to ncht correspond to 1 J X 2 
3X3, 5X5, 7X7, OX'* and 11 X 11 hole^ opened m the middle 
of the plate 

The curres of 1 to o X 5 hoV s:how a -second re*-oaance at *^00 
to 900 Hr. This has nothinc to do however with the re<»nalor 
the^ irreculanties being probablv due to a vibrntion of the me'al 
covenne plate * Ii manv holes are opened the fo*w on the me al 
plate Is nechcible and the efiect therefore, absent The <hnting 


1 re«=l o' J ran dea Eui ^ 

f “ J^tEUsT’O*' of pi5blj«*lKl. 

5 See aho T I.. Jordan, Jtwf Z, 5 (1*^^) ” 
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of this resonance to higher freqnenev if more holes are opened is 
explained by the decrease of the mass of the plate by opening 
the holes. 

Qualitatively we see: 

a. increasing the number of holes leads to a higher resonance 
frequency caused by the inerea.sed eonduetix-ity G of the resonator, 

b. if the number of holes increases, the air resistance of the 

resonator must decrease and therefore should increa.se (see 

Fig. 11, 81 and 82), as is actually found, 

c. the curve for 9X9 holes almost reaches a,, = 100 %, the 

cuiv'e 11 X 11 holes shows a decreased maximum absorption. 
Obviously the resistance reaches the critical value TT’’,, at about 
9X9 holes. 

It is not very easy to estimate the conductivity' as a function of 
the number of holes. Firstly the resonator is filled, so the struc- 
ture factor might be of influence. Furthermore the holes are so 
close to each other that the conducti^-ity is probably not propor- 
tional to the number of holes. Moreover the depth of the resonator 
is .small, so if many holes are opened the “effective diameter” of 
the perforations will be greater than the depth, so that then com- 
putation is quite impossible. 

One easy verification can be made, however. and Are,, are 

known for every curve. Since the volume is known the cum'e %vidth 
can be read from Fig. 82. The measured curve width agrees fairly 
well with the computed ones. The discrepancies are systematical 
and show the trend of a somewhat greater curve width than com- 
puted. Almost certainly this is due to a curvature of the impedance 
contour. 

We may conclude from the experimental results of this 
section that 

a. the simple theory holds for practical constructions, 

b. the curve width actually found is in most eases greater than 
that theoretically predicted. 

{ S COilBEN'ATIOXS OF EESOXATORS 

A drawback of any resonator construction is always that it ab- 
sorhs more or less selectively in the neighbourhood of resonance, 
although the best constructions of this type are very satisfactory' 



158 


ABSORPTION BT BESONATOE8 


VII 


jjideed 'Espeeisiiy in thos,, cases -nbere absorption at high frcquen 
cies as not necessary because the room in question alreadj contains 
a considerable amount of porous material (upholsterj of seats, 
auditorium, wood fibre plates, carpets) excellent results can be 
obtained mth a simple perforated panel device 
If however a resonator construction is needed that absorbs non 
selectively m a very broad frequency range recourse should be 
had to a combination of resonators Senes connections (Fig 87o) 
and parallel connections (Pig 89o) seem equallj advantageous in 
this respect 



Fig 87 

Double resonator u» senes 


Whereas a single perforated panel has only one degree of free 
dom (one resonance frequency), the double construction of Fig 87a 
has two, which manifests itself in two peaks in the absorption 
curve (Fig 87b) ^ 

The theory of such double resonators has been gnen by 
Rayleigh* For our purpose it is, however, more convenient 
to make use of the impedance concept If namely the two per 
forated panels have a non identical perforation we cannot always 
divide the whole construction into double resonators, each having 

1 P V Bruel, Lydtsolaitiyn og rumattKtii, p 120, Goteborg 1946 

2 Rayleigh, Theory of sound II, 2iid ed, J 310, Cambridge 1890 
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two holes, but the construction must be treated per m-. The con- 
struction is governed by the eonducthdties and G, per m- (in 
the volumes V^ and Fj per m’-' (equal to the thickness of 
the respective layers of air apart from the volume of the porous 
filling material, if any) and the resistances and if, per m^ (see 
Pig. STffi). 

The electrical analogue of the specific impedance is given in 
Pig. 88. In view of the six variables it seems hardly possible to 
use this scheme direetlj* or the computed impedance therefrom for 
designing adequate constructions. It is, however, easy to give a 
rough idea of the frequency range in which the construction ab- 
sorbs, viz., by considering the two values of resonance. 

The resonance frequencies are approximated by putting R^ and 
ifj equal to zero and consequently computing the frequencies for 



Pig. SS 

The impedance of double resonator in series connection (see Pig. 87). 


which the (imaginary) impedance becomes zero. Labeling the in- 
ductances and capacities for brevity Lj , L, , Ci , C, (see Pig. 88), 
we have (ifj = 0, if, = 0) 


which becomes zero for 


" li:) " 1^)-^ 


i 'f' L2C2 P LiGi , 1 — n (1 94.1 

L^L.GX: ^ 


Labeling the resonance angular frequencies wres.j , urca., , we im- 
mediately read from (7.24) 
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li Tr»e IS the r esQnancfe ang nlar frequency of the sjngle reso- 
nator (=c„]/G7y = l/KLC) 

From the second oquatitm (725) we see that the loganthmie 
mean value of the resonance freqoenctes is the same as when the 
resonators were not coupled at all, from the first equation (7 25) 
we see that the resonance frequencies are more or less separated 
due to the coupling This effect maj be neglected for the first 



Ftg 89 

Double resonator in parallel 


calculation or estimated to some 20 % of the free frequencies If 
for example a construction is wanted absorbing from 50 to 3000 Hz 
the coupled resonance frequencies will have to be about 150 and 
I(M\? Hz, the free resonance frequencies are to fie cAosen, tAere/ofe, 
at about 180 and 800 Hz This can be done wath the aid of the 
rules gi%en in § 5 It is more difficult to find the best adjust 
ment of the resistances B, and As long as the results of 

extensne experiments are lacking a good approximation will be 
the values that make the non coupled resonators not too selectn e 
A parallel connection unit is shown in Fig 89a It consists of 
two resonators quite near to each other with different resonance 
frequencj If the absorbing wall is supposed to be complctelj 
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CHAPTER VITI 


OBLIQUE AND RANDOM INCIDENCE PRACTICAL ASPECTS 
j 1 INTKODOCTrON 

In tho previons chapters ue confined oursches to absorption at 
normal incidence The reasons arc that both theory and experiment 
become much more difficult when considering oblique incidence 
Something must be said, however, about the behaviour of matcrnls 
at oblique mcidonco with a \icw to the great importance of this 
case m practice 

Tho measurement of sound absorption at normal incidence is 
quite simple As a consequence the properties of materials for this 
COSO have boon studied extensively The reverse is true for oblique 
incidence, the measuring techniques arc more complicated more 
over a new variable is introduced with the angle of incidence thus 
materially increasing the number of measurements necessary for 
describing the properties of o given matcnal Little is knoun, 
therefore from experiment of the bcliavioiir of materials at oblique 
incidence 

Of the measuring dmiccs the following ncod mentioning 

a exploration of tho interference pattern^ Analogous to the 
interference pattern in norma’ incidence interferometers a 
pattern is found for oblique incidence, from which the ab 
sorption coefficient can be denved Plane waves and largo 
samples are required., moreover the cxpenment should be 
carried out in free space m order to avoid disturbances of 
the sound field 

b Measurement at 45® incidence* If two walls of a room are 
covered mth the material under tost and a plane waic is 
supposed to impinge upon the material in a direction nor- 


1 L. Crewer, riek Tech, 10 (1033-) OOC, 13 (1030) 30 

2 K ScJiustor, JUiMl Z 3 (1088) 157 
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mallj' to the line of intersection of the walls and symmetric- 
ally to both walls an interference pattern is obtained. The 
absorption coefficient can be computed from the ratio of the 
extreme values of the pressure-gradient or the pressure along 
the direction of incidence. In this case too, free space and 
large samples are required. 

c. Measurement at 55.7° incidenceh If a plane wave impinges 
sj'mmetrically upon a three-plane-comer consisting of three 
mutually perpendicular walls covered noth the material under 
test, a similar pattern as sub b is obtained. The measurement 
is as sub b. 

d. Pulse method 2 . If a sound source, situated at a certain dis- 
tance from the sample (large size) gives a short sound pulse, 
the microphone at the reeehdng place will receive two pulses, 
if the experiment is carried out in free space, viz., a pulse 
due to the direct sound from source to microphone, and a 
pulse reflected from the sample. The ratio of the intensities 
is governed by geometrical quantities and the absorption coef- 
ficient of the sample, that can be calculated by adequately 
taking into account the geometrical configuration or by cali- 
bration \vith a totally reflecting surface. 

e. Standing waves in a rectangular room This method is 
based upon the weU-kno^vn fact that the eigenfunctions of a 
rectangular room with homogeneously covered walls can be 
described as the superposition of plane waves. This method 
seems to be the most promising one since no free space nor 
extremely large samples are necessarj% Under certain circum- 
stances it is possible to excite onlj’ one single characteristic 
frequency. Prom the decay of the intensity after stopping 
the sound source the absorption coefficient at the correspond- 
ing angle of incidence of the wall under test (other walls 
being totally reflecting) can be calculated. 

f. Reverberation method. In this well-known method the samples 


1 K. Schuster, Aku^sl. Z., 3 (193S) 137. 

= F. Spandock, Ann. Phyiil:, 20 (1934) 32S. 

3 H. and L. C r e m e r, Akni.it. Z., 2 (1937) 296. 
-1 H. Neubert, Akust. Z., 5 (1940) 189. 
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are used under more or less random meidenee The inter- 
pretation of the results is, however, extremely difficult ‘ 

Since we know onlj little from experiment about the behaviour 
of materials under oblique incidence, it seems more adequate to 
study the problem theoretically Among all matenals those witii 
so called local reaction are tiie easiest to consider A laj er of 
matenal is called locallj reacting® if the \elocitj component per 
pendicular to the surface only depends upon tlie pressure at the 
point and not upon the angle of incidence Typical examples are 
the samples with glass tabes (VI § 2), the single resonator in an 
infimte wall (VII § 2), a wall eo\ered with resonators of this tjpe, 
1 e, with resonators the volumes of which are not interconnected 
in another waj than via the air m front of the wall In these 
cases the ratio p/v„ is independent of the angle of incidence and 
equal to the impedance used for normal incidence 

In the case of a normal porous lajcr or of perforated panels 
with an uninterrupted air gap between panel and wall the im 
pedanee p/un will depend upon the angle of incidence, the more 
so when the sound veiocitj m the absorbing medium is greater 
and the damping of the travelling waves in a sideways direction 
13 smaller This is elearlj shown bj Morse® with different 
examples 

Now owing to the structure factor k of porous materials alwajs 
being much greater than unitj (saj 4 to 9) the sound velocit> 
in these media is 2 to 3 times smaller than in air Mo^eo^ er there 
is, as a rule, ample damping so there is reason for neglecting the 
extended reaction also in these eases That a low sound T,clocitj 
makes a material locally reacting can also be elucidated by con 
sidering the direction of the refracted wa\e in the material® The 
incident wave is refracted towards the normal as a light beam 
impinging upon jj glass mirfacr E\en at glancing incidence the 
refracted wave only deviates by 20° to 30° from the normal so 
that also m porous media the interconnection of the pores m a 
sidewajs direction wall not materially affect the impedance p/in 


1 See e g-, E Meyer and A Bebocb Z, 4 (1939) 51 

2PM Morse Vibration and Sound, 2nd ed , p 3dl, New lork 1948 
3C W Kosten, Thesis, p 99, Ddft 1942 
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If perforated panels are used -vrith non-fiUed air gaps extended 
reaction is quite sure. This seems to be a disadvantage of resonators 
of this type as compared with resonators -with filled gaps. Extended 
reaction can be eliminated of coirrse by splitting up the gap in 
many non-communicating volumes ’ 

§ 2 OBLIQUE LN'CIDEYGE OX LOCALLY EEACTIXG SURFACES 

Suppose a plane sound wave to impinge on the surface with 
impedance z, in a direction making an angle a with the normal 
to the surface (see Pig. 90). The reflected wave will be as shown 
in the figure. Let p, and Pr be the pressures at the surface in 
any point, then 



Fig. 90 

Oblique incidence 


if Vin and tVn are the normal components of the respective particle 
velocities towards the surface. Now 

1 Pr 

a.n = cos f and Vm — ^ cos <p, 

»» 0 ” 0 

and therefore 

Z cos o Pi ~ Pr 
~ P,—Pr' 

The complex reflection coefficient pjpi, and, therefore, the ab- 
sorption coefficient 1 — | Pr/p. | " bears the same relation to z cos e 
as to z in the case of normal incidence. Therefore, the circle- 

1 Xote added in proof: The behaviour of perforated panels nith non-fillcd 
air gaps at oblique incidence has recently been studied by J. B r 1 1 1 o u i n, 
Cahiers du Centre sci. et tech, du Bailment, Cahier 31, Jannaiy 1949. 
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diagram Fig 11 also appbes to this case, provided it is supposed 
to gue the absorption coefficient and phase jump in the complex 
2 cos plane 

A direct conclusion is the foUoxring since the impedance of 
manj materials is larger than both as to their real and imagin 
arj part, an apparent decrease of the impedance as a consequence 
of the multiplication bj cos if -will generally increase the absorption 
coefficient Onlj for -very small impedances (near resonance) or 
for \er 5 small cosy (almost glancing incidence) the diminution 
of z to scosy uill be disad\antagcons It, therefore, is phiisiblc 
tliat the mean value of the absorption coefficient, if all directions 
are taken into account, is greater than that for normal incidence 

An adequate mean lalue is found if the energj is considered 
that IS absorbed by a surface element that is exposed to a com 
pletely diffuse sound field m terms of the incident energj The 
incident energy comprised m a small solid angle d fi in a direction 
with an angle f to the normal is proportional to cosjxin, since 
the apparent surface of the surface element under consideration 
IS proportional to cosy Let the fraction o be absorbed Then the 
mean ^aluc o is defined is 

_ absorbed energj f cos y d n 
incident energj ^ cos y d n 

Non fl IS only dependent upon y, therefore take ns elementarj 
solid angle d n the small solid angle between a circular cone around 
the normal with top angle 2y and a similar cone with top angle 
2 (y + dy) , then dn = 2 - sin ydy and 
»/» 

a=y OySm 2 ydy (SOI) 

u 

For loeallj reacting surfaces is known as a function of y, 
so the integral can be e\alualed This seems to haie been done 
independontlj bj 31 o r s e et al ^ and the authors * If the im 
pedance z is written as f? + jX, we obtain 


1 P M Jlorse, R H Bolt, and R I* Brown, J Jcoiwf Soc Ami 
12 (1040) 217 

t C v, Kosten and C Zwikker, Ptyica 8, (Febr 1941) 251 




A more adequate representation is, of course, giving contours of 
constant a in the complex 2 -plane. This is done in Fig. 91, which 
is the parallel of Fig. 11 for random incidence. The contours can 



Fig. 91 

The mean value of the absorption coefficient of a locally reacting 
surface as a function of the impedance 


hardly be distinguished from circles, but actually are not. Another 
way of representation is given by Morse^. The maximum value 
of a turns out to be 0.95. In Fig. 92 a is compared to on the 
supposition that the impedance is real. The dotted part of the 
curve refers to impedances less than IFo. One sees that for im- 
pedances greater than TT,, a is appreciably larger than Uo- 


1 P. II. Morse, Vihratton and sound, 2ntl ed , p. 458, Itew York 194S. 
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$ 3 OBLIQUE INCIDENCE ON A S^^GLE BESONATOB 

Suppose a plane sound wave impinges obliquely upon a surface 
containing one single resonator Since tbe aeoiistioal ‘emf” is 
unchanged equal to 2p, and the internal source impedance is the 
radiation impedance as previously and moreover the resonator im 
pedance is unchanged too, the amount of the absorbed power is 
independent of the angle of incidence Dividing this absorbed 
power fay the power impinging npon 1 m* Wall surface, i e, bj 



Tie 92 

The companson between mean and oormsl ealue of the absorption 
coefficient for locally reacting surfaces with real impedance 

(p,^/2 IFfl) Cos <p yields the surface in m* of the equivalent open 
window, and since the absorbed power is independent of <p, the 
open window surface is inversely proportional to cos <p The average 
surface of open window for landom. incidence of sound within a 
cone with top angle 2f is readily computed to be 

(8 03) 

A 1 — cosf ’ 

if A denotes the mean value and A the surface of open window 
at normal incidence 
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y eannot be extended to 90°, since then A increases to infinity, 
which seems meaningless. There seems to be no theoretical ob- 
jection against this extension. It is only of no practical value since 
the effect could only be verified with samples of unlimited di- 
mensions. 

One conclusion is however certain. The effective number of m- 
open -window at random incidence wiU always be greater than the 
computed number at normal incidence; experiments in a revei'- 
beration room are in agreement hermvith. 

In order to give some idea about the magnitude of this effect, 
we give the ratio A/ A for a few values of the maximum angle 
of incidence considered, ymm: 

y„,ax 30° 45° 60° 75° 

A/A 1.07 1.18 1.38 1.82 

Bearing in mind that a practical resonator is always somewhat 
less selective than according to the theory as a consequence of the 
slightly curved impedance contour, it seems justified to expect the 
number of m^ open window in practice to be about two times 
higher than the theoretical value at normal incideime. 
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